HDU2110：Crisis of HDU（母函数）

Crisis of HDU
Time Limit: 3000/1000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others)
Total Submission(s): 6934 Accepted Submission(s): 2145

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

面对朋友的退出和资金短缺，海东集团能度过难关吗？作为老总的XHD如何来化解这场危机？欲知后事如何，且听下回分解——

这题的关键是要解决系数不超过1e4的问题，跟写快速幂的思想有点类似，就是每一步都要除以1e4，来确保得到的数一定不比1e4大。

还有一点就是如果总价值除以三不是整数的话就直接在一开始排除了啊，减少运行时间（虽然这题我试验了一下发现不一开始排除的话时间也不会增多：D），但还是养成好习惯吧。

下面附上ac代码：

#include <bits/stdc++.h>
#include <iostream>
#include <string.h>
#include <algorithm>
#include <ctype.h>
#include <cstdlib>
#include <cmath>
#include <queue>
#include <stdlib.h>
#include <time.h>
using namespace std;
#define INF 0x7ffffff
#define mod 10000019
typedef long long int ll;
ll sum,n;
ll c1[10100],c2[10100];
struct node {
    ll vlv,num;
} record[200];
int main() {
    while(cin >> n) {
        if(n == 0) {
            break;
        }
        sum = 0;
        for(int i = 1; i <= n; i++) {
            cin >> record[i].vlv >> record[i].num;
            sum += record[i].vlv * record[i].num;
        }
        if((sum / 3) * 3 != sum) { //无法平分为等量的三份
            cout << "sorry" << endl;
        } else {
            memset(c1, 0, sizeof(c1));
            memset(c2, 0, sizeof(c2));
            for(int i = 0; i <= record[1].vlv * record[1].num; i += record[1].vlv) { //解决第一个括号
                c1[i] = 1;
            }
            for(int i = 2; i <= n; i++) { //从第二个括号开始解决
                for(int j = 0; j <= sum; j++) {
                    for(int k = 0; k + j <= sum && k <= record[i].vlv * record[i].num; k += record[i].vlv) {
                        c2[j + k] += c1[j];
                        c2[j + k] = (c2[j + k] % 10000);
                    }
                }
                for(int j = 0; j <= sum; j++) {
                    c1[j] = c2[j];
                    c2[j] = 0;
                }
            }
            if(c1[sum/3] == 0) {
                cout << "sorry" << endl;
            } else {
                cout << c1[sum/3] << endl;
            }
        }
    }
    return 0;
}