BZOJ1044/HAOI2008木棍分割

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<algorithm>
#include<cmath>
#include<cstring>
#include<string>
using namespace std;
const int imax=50000+229;
const int mmax=1000+229;
const int mod=10007;
int n,m;
int a[imax],sum[imax];

void iread()
{
scanf("%d%d",&n,&m);
for(int i=1;i<=n;i++) scanf("%d",&a[i]);
for(int i=1;i<=n;i++) sum[i]=sum[i-1]+a[i];
}

bool pd(int delta)
{
int last=0; int need=0;
if(a[1]>delta) return false;
for(int i=2;i<=n;i++)
{
if(a[i]>delta) return false;
if(sum[i]-sum[last]>delta) { last=i-1; need++;}
}
if(need>m) return false; return true;
}

int f[2][imax],q[imax];
int get_ans(int delta)
{
int ans=0;
int pre,now,head,tail;
f[0][0]=1;
for(int i=1;i<=m;i++)
{
int now=i&1;
int pre=now^1;
memset(f[now],0,sizeof(f[now]));
head=tail=1; q[head]=0;
int tot=f[pre][0];
for(int j=1;j<=n;j++)
{
while(head<=tail && sum[j]-sum[q[head]]>delta)
{
tot=(tot-f[pre][q[head]]+mod)%mod;
head++;
}
q[++tail]=j; f[now][q[tail]]=tot;
tot=(tot+f[pre][q[tail]])%mod;
}
for(int j=n-1;j>=1;j--)
{
if(sum
-sum[j]<=delta) ans=(ans+f[now][j])%mod;
else break;
}
}
return ans;
}

void iwork()
{
int l=1; int r=1000*imax;
while(l<=r)
{
if(l+1==r || l==r)
{
if(!pd(l)) l=r;
break;
}
int Mid=(l+r)>>1;
if(pd(Mid)) r=Mid;
else l=Mid+1;
}
int ans=get_ans(l);
printf("%d %d\n",l,ans);
}

int main()
{
iread();
iwork();
return 0;
}