1. 程式人生 > >[bzoj2287][poj Challenge]消失之物_背包dp_容斥原理

[bzoj2287][poj Challenge]消失之物_背包dp_容斥原理

turn read zoj inline urn dot 允許 方案 efi

消失之物 bzoj-2287 Poj Challenge

題目大意:給定$n$個物品,第$i$個物品的權值為$W_i$。記$Count(x,i)$為第$i$個物品不允許使用的情況下拿到重量為$x$的方案數。

註釋:$1\le n,val_i\le 2\cdot 10^3$。


想法:只需要用取模瞎**容斥一下就行了。

Code:

#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
#define N 2010 
using namespace std;
int a[N],f[N],g[N];
inline char nc() {static char *p1,*p2,buf[100000]; return (p1==p2)&&(p2=(p1=buf)+fread(buf,1,100000,stdin),p1==p2)?EOF:*p1++;}
int rd() {int x=0; char c=nc(); while(!isdigit(c)) c=nc(); while(isdigit(c)) x=(x<<3)+(x<<1)+(c^48),c=nc(); return x;}
int main()
{
	// freopen("thing.in","r",stdin);
	// freopen("thing.out","w",stdout);
	int n=rd(),m=rd(); for(int i=1;i<=n;i++) a[i]=rd();
	f[0]=1;
	for(int i=1;i<=n;i++)
	{
		for(int j=m;j>=a[i];j--) (f[j]+=f[j-a[i]])%=10;
	}
	for(int i=1;i<=n;i++)
	{
		memset(g,0,sizeof g);
		g[0]=1;
		for(int j=1;j<=m;j++)
		{
			printf("%d",((f[j]-g[j%a[i]])%10+10)%10);
			g[j%a[i]]=((f[j]-g[j%a[i]])%10+10)%10;
		}
		puts("");
	}
	// fclose(stdin); fclose(stdout);
	return 0;
}

小結:好題。

[bzoj2287][poj Challenge]消失之物_背包dp_容斥原理