BZOJ 1044 HAOI2008 木棍切割 二分答案+动态规划

Posted mfmdaoyou

tags:

篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了BZOJ 1044 HAOI2008 木棍切割 二分答案+动态规划相关的知识,希望对你有一定的参考价值。

题目大意:给定n个连在一起的木棍。分成m+1段。使每段最大值最小,求最大值的最小值及最大值最小时切割的方案数

第一问水爆了……二分答案妥妥秒过

第二问就有些难度了 首先我们令f[i][j]表示用前j个棒♂子得到i段的方案数

诶我没打什么奇怪的符号吧

于是我们有动规方程

f[i][j]=Σf[i-1][k] (sum[j]-sum[k]<=ans,k<j)

这个最坏情况下是O(m*n^2)的,肯定挂

我们发现k的下界是单调上升的 于是我们直接令k为当前j时k的下界。开一个变量记录k~j的f值之和

每次j++时将k向后调整 此外f数组要开滚动数组不然MLE

写完交上去各种神慢……搞不懂其它人都写了啥

#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
#define M 50500
#define Mo 10007
using namespace std;
int n,m,limit,ans;
int maxlen=0,a[M],sum[M];
int f[2][M];
inline bool Judge(int x)
{
	int i,re=0,now=x;
	for(i=1;i<=n;i++)
	{
		if(now+a[i]>x)
			++re,now=a[i];
		else
			now+=a[i];
		if(re>m)
			return false;
	}
	return true;
}
inline int Bisection()
{
	int l=maxlen,r=sum[n];
	while(l+1<r)
	{
		int mid=l+r>>1;
		if( Judge(mid) )
			r=mid;
		else
			l=mid;
	}
	if( Judge(l) )
		return l;
	return r;
}
int main()
{
	int i,j,k;
	cin>>n>>m;++m;
	for(i=1;i<=n;i++)
		scanf("%d",&a[i]),sum[i]=sum[i-1]+a[i],maxlen=max(maxlen,a[i]);
	limit=Bisection();
	f[0][0]=1;
	for(i=1;i<=m;i++)
	{
		int temp=0;k=i-2;
		for(j=i;j<=n;j++)
		{
			temp+=f[~i&1][j-1],temp%=Mo;
			while(sum[j]-sum[k+1]>limit)
				temp+=Mo-f[~i&1][++k],temp%=Mo;
			f[i&1][j]=temp;
		}
		ans+=f[i&1][n],ans%=Mo;
	}
	cout<<limit<<' '<<ans<<endl;
}


以上是关于BZOJ 1044 HAOI2008 木棍切割 二分答案+动态规划的主要内容,如果未能解决你的问题,请参考以下文章

bzoj1044[HAOI2008]木棍分割 单调队列优化dp

BZOJ1044[HAOI2008]木棍分割 二分+DP

[BZOJ1044][HAOI2008]木棍分割 二分+贪心+dp+前缀和优化

BZOJ1044: [HAOI2008]木棍分割 二分+区间DP

bzoj 1044 [HAOI2008]木棍分割(二分+贪心,DP+优化)

BZOJ1044: [HAOI2008]木棍分割