bzoj2287. 【POJ Challenge】消失之物

方法一:

$\mathrm{f[x]}$ 表示所有物品凑成体积为 $\mathrm{x}$ 的方案数.  

$\mathrm{g[x][j]}$ 表示不用 $\mathrm{x}$ 物品组成体积为 $\mathrm{j}$ 的方案数.  

然后 $\mathrm{g}$ 数组可以用 $\mathrm{f,g}$ 容斥搞一波.   

时间复杂度为 $\mathrm{O(nm)}$.  

方法二:   

主要介绍一下分治解决 01 背包的方法.   

不难发现暴力算的话很多部分是重复的,而用线段树可以简化操作.   

令 $\mathrm{solve(l,r,dep)}$ 表示当前 $[l,r]$ 的物品没有加入,其他物品都加入的背包.    

每次进入左边分治前把右子树算进去.

每次进入右子树前撤销原本的右子树并加入左子树.   

时间复杂度: $O(nm \log n)$.   

#include <cstdio>
#include <vector>
#include <cstring>
#include <algorithm>
#define N  2005   
#define ll long long 
#define pb push_back 
#define setIO(s) freopen(s".in","r",stdin) 
using namespace std; 
const int mod=10;  
int n,m, v[N], f[20][N];  
void solve(int l, int r, int dep) {
    if(l == r) {
        for(int i=1;i<=m;++i) {
            printf("%d",f[dep-1][i]);  
        }
        printf("\n");  
        return ; 
    }
    int mid=(l+r)>>1;   
    // 复制.  
    for(int i=0;i<=m;++i) f[dep][i]=f[dep-1][i];       
    for(int i=mid+1;i<=r;++i) {
        for(int j=m;j>=v[i];--j) {
            (f[dep][j]+=f[dep][j-v[i]])%=mod;    
        }
    }        
    solve(l, mid, dep+1);   
    for(int i=0;i<=m;++i) f[dep][i]=f[dep-1][i];  
    for(int i=l;i<=mid;++i) {
        for(int j=m;j>=v[i];--j) 
            (f[dep][j]+=f[dep][j-v[i]])%=mod;   
    }
    solve(mid+1,r,dep+1);   
}
int main() {  
    // setIO("input"); 
    scanf("%d%d",&n,&m);  
    for(int i=1;i<=n;++i) {
        scanf("%d",&v[i]); 
    }
    f[0][0]=1;            
    solve(1, n, 1);   
    return 0; 
}

  

上一篇:P5200 [USACO19JAN]Sleepy Cow Sorting


下一篇:CC攻击