Crisis of HDU（母函数）

Crisis of HDU

Time Limit: 3000/1000 MS (Java/Others)    Memory Limit: 32768/32768 K (Java/Others)
Total Submission(s): 4258    Accepted Submission(s): 1197

Problem Description

话说上回讲到HDU大战东洋小苟，结果自然是中方大胜，这一战也使得海东集团在全球同行业中的地位更加巩固。随着集团的发展，很多创业时期的元老逐步功成身退，先是8600移民海外，然后是linle夫妇退隐山林，逐渐的，最初众多的元老只剩下XHD夫妇和Wiskey三人了。
到了2020年，因为扩张过度加上老鼠数量逐年减少，公司的发展遇到了前所未有的危机，此时集团已经没有任何流动资金，更可怕的是，这个时候，wiskey也决定退出了！
退出本身并不麻烦，麻烦的是，退出的人需要取走相应比例（1/3）金额的资产。
假设公司此时一共有n种价值的资产，每种价值的资产数量已知，请帮助心烦意乱的XHD夫妇计算一共有多少种分割资产的方法。

 

Input

输入包含多个测试实例，每个实例的第一行是一个整数n(n<100)，表示一共有n种价值的资产，接着的n行每行包含两个整数pi和mi(0<pi,mi<10)，分别表示某种价值和对应的数量，n为0的时候结束输入。

 

Output

对于每个测试实例，请输出分割资产的方案数%10000，如果不能分割，请输出“sorry”，每个实例的输出占一行。

 

Sample Input

2 1 1 2 1 0

 

Sample Output

1

代码：

 #include<stdio.h>
 #include<string.h>
 #define mem(a) memset(a,0,sizeof(a))
 const int MAXN=;
 int main(){
     int n,p[MAXN],m[MAXN],a[MAXN],b[MAXN];
     while(~scanf("%d",&n),n){
         int sum=;
         mem(a);mem(b);
         for(int i=;i<=n;i++)
             scanf("%d%d",&p[i],&m[i]),sum+=m[i]*p[i];
         for(int i=;i<=m[];i++)
             a[p[]*i]=,b[i]=;//p[1]
             if(sum%!=)puts("sorry");
     else{
         sum/=;
         for(int i=;i<=n;i++){
             for(int j=;j<=sum;j++)
                 for(int k=,t=;k<=sum&&t<=m[i];k+=p[i],t++)
                     b[j+k]+=a[j];
             for(int j=;j<=sum;j++)
                 a[j]=b[j]%,b[j]=;
         }
     //    for(int i=0;i<=sum;i++)printf("%d ",a[i]);
     //    puts("");
         if(a[sum]==)puts("sorry");
         else printf("%d\n",a[sum]);
         }
     }
     return ;
 }

题解：取走总资产的1/3；

该题的母函数为：

extern "C++"{
#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<queue>
using namespace std;
const int INF = 0x3f3f3f3f;
#define mem(x,y) memset(x,y,sizeof(x))
typedef long long LL;
typedef unsigned long long ULL;

void SI(int &x){scanf("%d",&x);}
void SI(double &x){scanf("%lf",&x);}
void SI(LL &x){scanf("%lld",&x);}
void SI(char *x){scanf("%s",x);}

}
const int MAXN = 10010;
int a[MAXN],b[MAXN];
int p[MAXN],m[MAXN];
int main(){
	int n;
	while(scanf("%d",&n),n){
		int sum = 0;
		for(int i = 0;i < n;i++)scanf("%d%d",&p[i],&m[i]),sum += p[i] * m[i];
		mem(a,0);mem(b,0);
		if(sum % 3 != 0){
			puts("sorry");continue;
		}
		for(int i = 0;i <= m[0];i++){
			a[p[0] * i] = 1;
		}
		for(int i = 1;i < n;i++){
			for(int j = 0;j <= sum ;j++){
				for(int k = 0;k <= m[i];k++){
					b[k * p[i] + j] += a[j];
				}
			}
			for(int j = 0;j <= sum;j++){
				a[j] = b[j] % 10000; b[j] = 0;
			}
		}
		if(a[sum/3] == 0)puts("sorry");
		else printf("%d\n",a[sum/3] % 10000);
	}
	return 0;
}