Another Easy Problem fzu1753

Another Easy Problem

Time Limit: 1000MS

Memory Limit: 32768KB

64bit IO Format: %I64d & %I64u

[Submit]   [Go Back]   [Status]  

Description

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

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

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

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

................

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

xtt希望你告诉他满足条件的最大的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

#include<stdio.h>
#include<string.h>
#include<stdlib.h>
#include<time.h>
#include<iostream>
#include<algorithm>
using namespace std;
int p[100000]= {0},t,pp[10000];
void init()
{
    long long i,j;
    t=0;
    for(i=2; i<1000; i++)
    {
        if(!p[i])
        {
            p[t++]=i;
            j=i*i;
            while(j<100000)
            {
                p[j]=1;
                j+=i;
            }
        }
    }
}
void fun(int n,int k)
{
    int i,j,m,sum;
    for(i=0; i<t; i++)
    {
        m=n;
        sum=0;
        while(m)
        {
            m/=p[i];
            sum+=m;
        }

        m=k;
        while(m)
        {
            m/=p[i];
            sum-=m;
        }

        m=n-k;
        while(m)
        {
            m/=p[i];
            sum-=m;
        }
        if(pp[i]!=-1)
        pp[i]=pp[i]>sum?sum:pp[i];
        else pp[i]=sum;
    }
}
int pow3(int a,int b)//求a^b
{
    int r = 1,base = a;
    while(b!= 0 )
    {
        if(b&1)
            r*= base;
        base*= base;
        b>>= 1;
    }
    return r;
}

int main()
{
    init();
    int tt,m,n,i;
    while(~scanf("%d",&tt))
    {
        memset(pp,-1,sizeof(pp));
        while(tt--)
        {
            scanf("%d%d",&m,&n);
            fun(m,n);
        }
        long long sum=1;
        for(i=0; i<t; i++)
        {
            if(pp[i])
                sum*=pow3(p[i],pp[i]);
        }
        printf("%I64d
",sum);
    }
}