hdu 2110 基础母函数

题意：退出本身并不麻烦，麻烦的是，退出的人需要取走相应比例（1/3）金额的资产。
假设公司此时一共有n种价值的资产，每种价值的资产数量已知，请帮助心烦意乱的XHD夫妇计算一共有多少种分割资产的方法。

 

现在我们引用《组合数学》上最经典的一个例题:

我们要从苹果、香蕉、橘子和梨中拿一些水果出来，要求苹果只能拿偶数个，香蕉的个数要是5的倍数，橘子最多拿4个，梨要么不拿，要么只能拿一个。问按这样的要求拿n个水果的方案数。

g(x)=(1+x^2+x^4+...)(1+x^5+x^10+..)(1+x+x^2+x^3+x^4)(1+x)  就是最终的生成函数

 

 这个母函数让我想起了去年暑假集训，那个时候母函数搞了几天都不懂，bfs完全不会，dijkstra搞不懂就反复裸拍代码，acm入门确实很难，但是一旦入了门就能感受到其中的乐趣，学的越多，感觉自己的知识就越匮乏，也就想学的更多。正是这样的循环让我的代码能力得以提升，孜孜不倦的去追求更多的知识。

 1 #include<cstdio>

 2 #include<iostream>

 3 #include<algorithm>

 4 #include<cstring>

 5 #include<cmath>

 6 #include<queue>

 7 #include<map>

 8 using namespace std;

 9 #define MOD 1000000007

10 const int INF=0x3f3f3f3f;

11 const double eps=1e-5;

12 #define cl(a) memset(a,0,sizeof(a))

13 #define ts printf("*****\n");

14 const int MAXN=105;

15 int n,tt;

16 int c1[10010],c2[10010],m[MAXN],p[MAXN];

17 int main()

18 {

19     int i,j,k;

20     #ifndef ONLINE_JUDGE

21     freopen("1.in","r",stdin);

22     #endif

23     while(scanf("%d",&n)!=EOF)

24     {

25         if(n==0)    break;

26         int sum=0;

27         for(i=1;i<=n;i++)

28         {

29             scanf("%d%d",&p[i],&m[i]);

30             sum+=p[i]*m[i];

31         }

32         if(sum%3!=0)

33         {

34             printf("sorry\n");

35             continue;

36         }

37         sum/=3;

38         cl(c1),cl(c2);

39         for(i=0;i<=p[1]*m[1];i+=p[1])    c1[i]=1;   //初始化第一个函数

40         for(i=2;i<=n;i++)   //需要乘的次数

41         {

42             for(j=0;j<=sum;j++) //生成函数

43             {

44                 for(k=0;k<=p[i]*m[i]&&(k+j)<=sum;k+=p[i])   //相乘的函数

45                 {

46                     c2[k+j]+=c1[j];

47                 }

48             }

49             for(j=0;j<=sum;j++)

50             {

51                 c1[j]=c2[j]%10000;

52                 c2[j]=0;

53             }

54         }

55         if(c1[sum]==0)

56             printf("sorry\n");

57         else printf("%d\n",c1[sum]);

58     }

59 }