POJ 1742 Coins(多重背包|单调队列优化)

用二分拆包就TLE了，学了个单调队列优化，具体东西写在代码的注释中。

奋战LTC的男人八题，两天两道题，成为1/4男人，3/4不明生物，所以。。。。至今情况不知道算是阴阳人还是人妖。。。。

底下的如果没解完八题的仔细琢磨下再喷我。。。。

/*
    这是我至今遇到的最恶心的多重背包，居然可以把二分拆包卡掉，ORZ LTC。。。
    不过确实很有意思，题意为有n种钱币，每种钱币面值a[i],数量为c[i],问组成的总钱数不超过m的组合数
    首先第一反应，二进制优化，交上去TLE。。。上网搜来搜去，得到结果——>单调队列优化
    优化的思想如下：
    如果c[i]=1 01背包直接解
    若a[i]*c[i]>=m 等价于完全背包
    剩下的用单调队列优化：
        若当前在尝试第i种硬币，那么dp[i]可以出现的情况就是存在一个j
        dp[j]=1&&(i-j)%a[i]==0&&(i-j)/a[i]<=c[i]
        说的直白点，就是在算步长，在一定距离内等步长的几个点如果存在真，那么这个情况也为真
        把原来的情况存在队列que中，在更新的时候就不会影响到原来情况
        而用sum统计当前队列种的值得和
        如果sum>0就可以知道存在这种情况
        从左向右扫，维护que和sum
*/
#include<iostream>
#include<cstdio>
#include<cstring>
using namespace std;

#define N 100
#define M 100010
int a[N];
int c[N];
bool dp[M],que[M];

int main(){
    int i,j,k,n,m,tot;
    while(scanf("%d %d",&n,&m)!=EOF&&n+m){
        for(i=0;i<n;i++){
            scanf("%d",&a[i]);
        }
        for(i=0;i<n;i++){
            scanf("%d",&c[i]);
        }
        memset(dp,0,sizeof(dp));
        dp[0]=1;
        int ret=0;
        for(i=0;i<n;i++){
            if(c[i]==1){    //01背包
                for(j=m;j>=a[i];j--){
                    if(!dp[j]&&dp[j-a[i]]){
                        dp[j]=1;
                        ret++;
                    }
                }
            }
            else if(c[i]*a[i]>=m){   //完全背包
                for(j=a[i];j<=m;j++){
                    if(!dp[j]&&dp[j-a[i]]){
                        dp[j]=1;
                        ret++;
                    }
                }
            }
            else{
                for(j=0;j<a[i];j++){
                    int s=0,e=-1,sum=0;
                    for(k=j;k<=m;k+=a[i]){
                        if(s+c[i]==e){
                            sum-=que[s++];
                        }
                        sum+=dp[k];
                        que[++e]=dp[k];
                        if(!dp[k]&&sum){
                            dp[k]=1;
                            ret++;
                        }
                    }
                }
            }
        }
        printf("%d\n",ret);
    }
    return 0;
}