A.chinese remainder theorem
套刘汝佳书上的中国剩余定理模板,注意无解的检验方法是将模板的结果依次代入数组各元素试验是否正确,还有套模板需确保各元素不重复。以前没用过这个定理,这个是看别人的blog了解的注意事项。
#include<cstdio>
#include<algorithm>
using namespace std;
long long t,a[105],m[105],n;
bool v[20][20];
void gcd(long long a,long long b,long long &d,long long&x,long long&y)
{
if(!b)
{
d=a;
x=1;y=0;
}
else
{
gcd(b,a%b,d,y,x);
y-=x*(a/b);
}
}
long long china()
{
long long M=1,d,y,x=0;
for(long long i=0;i<n;i++)M*=m[i];
for(long long i=0;i<n;i++)
{
long long w=M/m[i];
gcd(m[i],w,d,d,y);
x=(x+y*w*a[i])%M;
}
return (x+M)%M;
}
int main()
{
scanf("%d",&t);
while(t--)
{
for(int i=0;i<20;i++)for(int j=0;j<20;j++)v[i][j]=0;
scanf("%lld",&n);
for(long long i=0;i<n;i++){scanf("%lld%lld",&a[i],&m[i]);if(v[a[i]][m[i]]==1)i--,n--;v[a[i]][m[i]]=1;}
long long ans=china();
bool ok=1;
for(int i=0;i<n;i++)if(ans%m[i]!=a[i])ok=0;
if(ok)printf("%lld\n",ans);
else printf("IMPOSSIBLE\n");
}
return 0;
}
第二次做,其实暴力就可以过。枚举最大范围是18以内的素数乘积,只有不到60w,纯暴力就可以过,自己一开始枚举的最大范围是17!显然不行。
#include<iostream>
using namespace std;
int t,n,a[20],p[20];
int x,y;
bool ok;
int main()
{
cin>>t;
while(t--)
{
cin>>n;
for(int i=0;i<n;i++)cin>>a[i]>>p[i];
y=2*3*5*7*11*13*17;
for(x=0;x<y;x++)
{
ok=1;
for(int i=0;i<n;i++)if(x%p[i]!=a[i]){ok=0;break;}
if(ok)break;
}
if(x<y)cout<<x<<"\n";
else cout<<"IMPOSSIBLE\n";
}
return 0;
}
B.AC的概率
记录当前遇到的数中取模为i的数的个数,扫一遍数组。
#include<cstdio>
#include<cstring