bzoj2287【POJ Challenge】消失之物 缺一01背包

bzoj2287【POJ Challenge】消失之物 缺一01背包

链接

bzoj

思路

分治solve(l,r,arr)表示缺少物品\([l,r]\)的dp数组arr。

然后solve(l,mid,arr)用右边的物品更新，solve(mid+1,r,arr)同理。

\(f(n)=2*f(\frac{n}{2})+(r-l+1)*m\)

复杂度为\(O(nmlog{n})\)
缺点最短路也是这样，用\(floyd\)

代码

#include <bits/stdc++.h>
using namespace std;
const int _=2e3+7;
int n,m,w[_],tmp[_<<2][_],ans[_],cnt;
void solve(int l,int r,int *f) {if(l==r) {for(int i=1;i<=m;++i) printf("%d",f[i]);puts("");return;}int mid=(l+r)>>1,p=++cnt;for(int i=0;i<=m;++i) tmp[p][i]=f[i];for(int i=mid+1;i<=r;++i)for(int j=m;j>=w[i];--j)f[j]+=f[j-w[i]],f[j]%=10;solve(l,mid,f);for(int i=0;i<=m;++i) f[i]=tmp[p][i];for(int i=l;i<=mid;++i)for(int j=m;j>=w[i];--j)f[j]+=f[j-w[i]],f[j]%=10;solve(mid+1,r,f);   
}
int main() {scanf("%d%d",&n,&m);for(int i=1;i<=n;++i) scanf("%d",&w[i]);tmp[0][0]=1,solve(1,n,tmp[0]);return 0;
}