poj 3260 双端队列优化多重背包+完全背包

#include<cstdio>
#include<cstring>
#define MAX(x,y) ((x)>(y)?(x):(y))
#define MIN(x,y) ((x)>(y)?(y):(x))
#define MAX_N 100200
using namespace std;
int dp1[MAX_N],dp2[MAX_N];
int n,m,num;
int c[MAX_N],v[MAX_N];
void solve()
{
	memset(dp2,0x3f,sizeof(dp2));
	dp2[0]=0;
	for(int i=0;i<n;i++)
	{
		for(int j=v[i];j<=m+num;j++)
		{
			dp2[j]=MIN(dp2[j],dp2[j-v[i]]+1);
		}
	}
}
void solve1()
{
	int deqv[100000],deq[100000];
	memset(dp1,0x3f,sizeof(dp1));
	dp1[0]=0;
	for(int i=0;i<n;i++)
	{	 
		for(int a=0;a<v[i];a++)
		{
			int s=0,t=0;
			for(int j=0;j*v[i]+a<=m+num;j++)
			{
				int val=dp1[j*v[i]+a]-j;
				while(s<t&&deqv[t-1]>=val)
				t--;
				deq[t]=j;
				deqv[t++]=val;
				dp1[j*v[i]+a]=deqv[s]+j;
				if(deq[s]==j-c[i])
				s++;
			}
		}
	}
}
int main()
{
	scanf("%d%d",&n,&m);
	num=0;
	for(int i=0;i<n;i++)
	{
	   scanf("%d",&v[i]);
	   num=MAX(num,v[i]);
    }
    num*=num;
	for(int i=0;i<n;i++)
	scanf("%d",&c[i]);
	solve1();
	solve(); 
	int res=0x3f3f3f3f;
	for(int i=num+m;i>=0;i--)
	{
		res=MIN(res,dp1[i+m]+dp2[i]);
	}
	if(res==0x3f3f3f3f)
	printf("-1\n");
	else
	printf("%d\n",res);
	return 0; 
}