BZOJ 1044: [HAOI2008]木棍分割

最新推荐文章于 2014-12-11 22:51:02 发布

Killerfear

最新推荐文章于 2014-12-11 22:51:02 发布

阅读量735

点赞数

分类专栏：动态规划

本文链接：https://blog.csdn.net/wsx1754175/article/details/38980211

版权

动态规划专栏收录该内容

74 篇文章 0 订阅

订阅专栏

题目：http://www.lydsy.com/JudgeOnline/problem.php?id=1044

分析：第一问直接用二分贪心就行了，不说了。第二问有dp[i][j]=sum{ dp[k][j-1] | 1<=k<j && sum[i]-sum[k] <= maxlen } 首先注意到，i在往右移的同时，min(k) 只能往右移或者不动，我们只需要维护dp[min(k)][j-1] ~ dp[i-1]j-1]的和就行了。

我的方法比较麻烦。。。经常会有小错误，临界想得我头都大了。。。

代码：

#include <iostream>
#include <cstdio>
#include <cstring>
#include <stack>
#include <string.h>
#include <queue>
#include <vector>
#include <algorithm>
#include <cassert>
#include <set>
#include <map>
#include <cmath>
#include <ctime>
using namespace std;
#define rep(i,a,b) for(int i=(a);i<(b);++i)
#define rrep(i,a,b) for(int i=(a);i>=(b);--i)
#define clr(a,x) memset(a,(x),sizeof(a))
#define eps 1e-8
#define LL long long
#define mp make_pair
const int maxn=50000+5;
const int mod=10007;
int L[maxn],n,m;
int f[1005];
int sum[1005];
queue<int> tmp[1005];
//dp[i][j] = sum(dp[k][j-1]); k<j && sum[k+1][j]<=maxsum
int maxlen=0;

void init()
{
    rep(i,0,1005) {
        f[i]=sum[i]=0;
        while(tmp[i].size()) tmp[i].pop();
    }
}

bool ok(int up)
{
    int k=0,S=0;
    rep(i,1,n+1) {
        S += L[i];
        if(S>up) S=L[i],++k;
        if(k>m) return false;
    }
    return true;
}

void up_len()
{
    int l=0,r=0;
    rep(i,1,n+1) r += L[i] , l = max(l,L[i]);
    int m;
    while(l<=r) {
        m=(l+r)>>1;
        if(ok(m)) { maxlen=m; r=m-1; }
        else l=m+1;
    }
}


int main()
{
    //freopen("in.txt","r",stdin);
    while(scanf("%d%d",&n,&m)==2) {
        rep(i,1,n+1) scanf("%d",L+i);
        up_len();
        printf("%d ",maxlen);
        init();
        int h=1,S=L[1];
        sum[0]=1; f[0]=1;
        rep(j,0,m) tmp[j].push(f[j]);
        rep(i,2,n+1) {
            S+=L[i];
            while(S>maxlen) {
                f[0]=0;
                if(S-L[h]<=maxlen) break;
                S -= L[h++];
                rep(j,0,m) {
                    sum[j] -= tmp[j].front();
                    tmp[j].pop();
                    if(sum[j]<0) sum[j]+=mod;
                }
            }
            rrep(j,m,0) {
                if(j-1>=0&&j<=i-1) f[j]=sum[j-1];
                tmp[j].push(f[j]);
                sum[j] += f[j];
                if(sum[j]>=mod) sum[j] -= mod;
            }
        }
        int ans=0;
        rep(i,0,m+1) ans = (ans+f[i])%mod;
        ans = (ans+mod)%mod;
        printf("%d\n",ans);
    }
}