fzu1753 Another Easy Problem（多组C(n,m)求解gcd）

传送门：http://acm.fzu.edu.cn/problem.php?pid=1753

Problem 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

思路

  多组C(n,m)求解gcd模板https://blog.csdn.net/qq_43724031/article/details/100901070
唯一注意的是，输出，debug了好久……OJ编译器问题，%lld不被认可，可以%I64d或者直接cout

#include<cstdio>
#include<iostream>
#include<cmath>
#include<algorithm>
#include<cstring>
#include<cstdlib>
#include<cctype>
#include<vector>
#include<stack>
#include<queue>
#include<map>
#include<ctime>
#define ll long long
#define ld long double
#define ull unsigned long long
using namespace std;
const int inf = 0x3f3f3f3f;
const int maxn = 100100;
struct In{
    int n,m;
}a[maxn];
int pri[maxn],tot;
bool v[maxn];
int p[maxn];
void prime()
{
    memset(v,true,sizeof(v));
    tot = 0;
    for(int i=2;i<maxn;i++){
        if(v[i])    pri[++tot] = i;
        for(int j=1;j<=tot&&i*pri[j]<maxn;j++){
            v[i*pri[j]] = false;
            if(i%pri[j]==0) break;
        }
    }
}
int cal(int n,int p)
{
    int ans = 0;
    while(n){
        ans += n/p;
        n /= p;
    }
    return ans;
}
int getgcd(int n,int m,int p)
{
    int q = cal(n,p) - cal(m,p) - cal(n-m,p);
    return q;
}
int main(void)
{
    int t,m,n;
    prime();
    while(~scanf("%d",&t)){
        int _min = inf;
        for(int i=0;i<t;i++){
            scanf("%d%d",&a[i].n,&a[i].m);
            _min = min(_min,a[i].n);
        }

        ll ans = 1;
        for(int i=1;pri[i]<=_min;i++){
            p[i] = inf;
            for(int j=0;j<t;j++){
                int r = getgcd(a[j].n,a[j].m,pri[i]);
                p[i] = min(p[i],r);
            }
            for(int j=0;j<p[i];j++)
                ans *= (ll)pri[i];
        }
        printf("%I64d\n",ans);	
        //cout<<ans<<'\n';	
    }
    return 0;
}