HDU 2110——Crisis of HDU

话说上回讲到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

这道题是一个母函数的题，利用母函数模板就可以解决问题，注意在j循环的时候不是每次都会达到最大值，而是相应的最大值是相应的cmax+=p[i]*m[i]，

#include <iostream>
#include <cstring>
using namespace std;
int n,p[105],m[105];
int ans[100005],tmp[100005];
int sum;
void mu()
{
    memset(ans,0,sizeof(ans));
    memset(tmp,0,sizeof(tmp));
    ans[0]=1;
    int cmax=0;
    for(int i=1;i<=n;i++)
    {
        cmax+=p[i]*m[i];
        for(int j=0;j<=cmax;j++)
        for(int k=0;k<=m[i]&&k*p[i]+j<=cmax;k++)
        {
            tmp[j+k*p[i]]+=ans[j];
            tmp[j+k*p[i]]%=10000;
        }
        memcpy(ans,tmp,sizeof(tmp));
        memset(tmp,0,sizeof(tmp));
    }


}
int main()
{


    while(cin>>n)
    {
        if(!n)
        break;
        sum=0;
        for(int i=1;i<=n;i++)
        {
            cin>>p[i]>>m[i];
            sum+=p[i]*m[i];
        }
        if(sum%3)
        {cout<<"sorry"<<endl;continue;}
        else
        {
            sum/=3;
            mu();
            if(ans[sum]==0)
            cout<<"sorry"<<endl;
            else
            cout<<ans[sum]<<endl;


        }
    }
    //cout << "Hello world!" << endl;
    return 0;
}