单调队列优化多重背包

 

#include <iostream>
#include <cstring>
#include<algorithm>
#include <cmath>
using namespace std;
int dp[1000];
int main()
{
    int n,m,s,v,w;
    cin>>n>>m;
    for(int i=1;i<=n;i++)
    {
        cin>>v>>w>>s;
        for(int j=m;j>=v;j--)
        {
            for(int k=0;k<=s&&j-k*v>=0;k++)
            {
                dp[j]=max(dp[j],dp[j-k*v]+k*w);
//每一个dp[j]都是由v的倍数转移过来的因此我们可以通过余数把它们分成相同的组然后组之内进行状态转移。组与组之间没有转移。
            }
        }
    }
    cout<<dp[m]<<endl;
}

 

//板子ac码
#include <iostream>
#include <cmath>
#include <string.h>
#include <algorithm>
#include <map>
#include <queue>
typedef  long long ll;
using namespace std;
int n,m;
int g[200000],dp[200000],q[200000];
int main()
{
    cin>>n>>m;
    for(int i=1;i<=n;i++)
    {
        int s,v,w;
        cin>>v>>w>>s;
        memcpy(g,dp,sizeof(dp));
        for(int j=0;j<v;j++)
        {
            int ll=1,rr=0;//每个组都新建一个单调队列维护最优的决策的下标
            for(int k=j;k<=m;k+=v)
            {
                dp[k]=g[k];
                if(ll<=rr&&k-s*v>q[ll])
                    ll++;//取队头操作如果长度超了的话就说明不符合了
                //这时候就应该把队头剔除
                if(ll<=rr)
                    dp[k]=max(dp[k],g[q[ll]]+(k-q[ll])/v*w);//用最优的决策来维护
                while(ll<=rr&&g[q[rr]]-(q[rr]-j)/v*w<=g[k]-(k-j)/v*w)rr--;//这里的减号 
//即"还原".(个数计算,因为k是下标而且枚举的都是余数为j的下标那么他们肯定在减去j之后都能被v整除也 
//就是个数。
//新决策入队把那些不可能用到的决策都删了
                q[++rr]=k;
            }
        }
    }
    cout<<dp[m]<<endl;
}

看了一天yxc的讲解和课本总算是能看懂了觉得最关键的还是看出来同余

代码可能有些许错误 如果发现还望指出。

顺便附上二进制优化码

#include <iostream>
#include <cmath>
#include <string.h>
#include <algorithm>
#include <map>
#include <queue>
typedef  long long ll;
using namespace std;
int cnt=0;
int w[1000040],v[1000040],dp[3000];
int main()
{
    int x,y,n,c,m;
    cin>>n>>m;
    for(int i=1;i<=n;i++)
    {
        cin>>x>>y>>c;
        for(int k=1;k<=c;k*=2)
        {
            c-=k;
            w[++cnt]=y*k;
            v[cnt]=x*k;
        }
        if(c>0)
        {
            w[++cnt]=y*c;
            v[cnt]=x*c;
        }
    }
    memset(dp,0,sizeof(dp));
    dp[0]=0;
    for(int i=1;i<=cnt;i++)
    {
        for(int j=m;j>=v[i];j--)
        {
            dp[j]=max(dp[j],dp[j-v[i]]+w[i]);
        }
    }
    cout<<dp[m]<<endl;
}