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;
}