杭电ACM——1171，Big Event in HDU（母函数）

这道题虽然不是计数问题，但同样可以用母函数解决，好在这道题的时间要求较宽。

代码如下：

#include<cstdio>
#include<algorithm>
#include<iostream>
#include<cstring>
typedef long long ll;
using namespace std;
ll c1[250005],c2[250005];
struct fac
{
       ll v;
       ll n;
}facs[55];
bool cmp(fac x,fac y)
{
     return x.v<y.v;
}
int main()
{
    int N,ma,i,j,k;
    ll A,B;
    while((cin>>N)&&(N>0))
    {
        ma=0;
        for(i=1;i<=N;i++)
        {
            cin>>facs[i].v>>facs[i].n;
            ma+=facs[i].v*facs[i].n;     //ma记录所有facilities的总值
        }
        sort(facs+1,facs+N+1,cmp); //从小到大排序，按v
        for(i=0;i<=ma;i++)       
     		c1[i]=c2[i]=0;
        for(i=0;i<=facs[1].n;i++)  //母函数第一个多项式系数初始化为1
     	    c1[i*facs[1].v]=1;
        for(i=2;i<=N;i++)
        {
            for(j=0;j<=ma;j++)
                for(k=0;(k+j<=ma)&&(k<=facs[i].v*facs[i].n);k+=facs[i].v)
                    c2[j+k]=c2[j+k]+c1[j];
            for(j=0;j<=ma;j++)
            {
                 c1[j]=c2[j];
                 c2[j]=0;
            }
        }
        /*for(i=1;i<=N;i++)
        {
           for(j=1;j<=facs[i].n;j++)
            cout<<facs[i].v*j<<"  "<<c1[facs[i].v*j]<<endl;
        }*/
        for(i=0;i<=ma;i++)
        {
        	if(c1[i])   //如果c1[i]不为0，说明v为i的方案存在
        	{//cout<<c1[i]<<' ';
        		A=i;
        		B=ma-A;
        		if(A>=B)      //因为i是从小到大，所以如果发现A>=B，就可以输出A,B，直接结束循环了
        		{
        			cout<<A<<' '<<B<<endl;
        			break;
				}
			}
		}
    }
	return 0;
}