FZU1753:组合数因数分解

Description
小TT最近学习了高斯消元法解方程组,现在他的问题来了，如果是以下的方程，那么应该如何解呢?

C(n1,m1)==0 (mod M)

C(n2,m2)==0 (mod M)

C(n3,m3)==0 (mod M)

…………….

C(nk,mk)==0 (mod M)

小TT希望你告诉他满足条件的最大的M

其中C(i,j)表示组合数,例如C(5,2)=10,C(4,2)=6…

Input
输入数据包括多组，每组数据的第一行是一个正整数T(1<=T<=150)表示接下来描述的T个方程

接下来T行，每行包括2个正整数ni,mi (1<=mi<=ni<=100000)

Output
输出一行答案，表示满足方程组的最大M。

Sample Input
3
100 1
50 1
60 1
Sample Output
10

思路：
看似高斯消元，实则因数分解，因为方程数小于未知数的个数，必定解不出来；
刚开始TLE，后来看了人家的博客才发现了一个可以优化的地方；
只需要考虑比最小的n小的因子！！！！！！（打表不可能）

#include <cstdio>
#include <cstring>
#include <algorithm>
#include <iostream>
#define inf 0x3f3f3f3f
using namespace std;
const int maxn=100001;
int prime[maxn+1];
int num[maxn+1],ans[maxn+1],Min;
void getprime()
{
    memset(prime,0,sizeof prime);
    for(int i=2; i<=maxn; i++)
    {
        if(!prime[i]) prime[++prime[0]]=i;
        for(int j=1; j<=prime[0]&&(long long)i*prime[j]<=maxn; j++)
        {
            prime[i*prime[j]]=1;
            if(i%prime[j]==0) break;
        }
    }
}
void cal(int n,int op)
{
    int x=1;
    for(int i=prime[1]; i<=Min; i=prime[++x])//只需要考虑比Min小的因子！！！！！
    {
        int ans=0;
        int t=n;
        while(t)
        {
            ans+=t/i;
            t/=i;
        }
        if(op) num[i]+=ans;
        else num[i]-=ans;
    }
}
int main()
{
    getprime();
    int t;
    int a[159],b[159];
    while(~scanf("%d",&t))
    {
        memset(ans,inf,sizeof ans);
        for(int i=0; i<t; i++)
            cin>>a[i]>>b[i];
        Min=inf;
        for(int i=0; i<t; i++)
            Min=min(Min,a[i]);
        for(int i=0; i<t; i++)
        {
            memset(num,0,sizeof num);
            cal(a[i],1);
            cal(b[i],0);
            cal(a[i]-b[i],0);
            for(int i=prime[1],x=1; i<=Min; i=prime[++x])
                ans[i]=min(ans[i],num[i]);
        }
        long long r=1;
        for(int i=prime[1],x=1; i<=Min; i=prime[++x])
            if(ans[i])
            {
                for(int j=1; j<=ans[i]; j++)
                    r*=i;
            }
        cout<<r<<endl;
    }
    return 0;
}