首页 > 代码库 > 【数论】【中国剩余定理】poj1006 生理周期
【数论】【中国剩余定理】poj1006 生理周期
CRT用于求解一元线性同余方程组(模数互质),实际上模数不互质我们也可以解决,在之前的某篇文章里提过。如下
http://www.cnblogs.com/autsky-jadek/p/6596010.html
#include<cstdio> using namespace std; typedef long long ll; ll m[4],a[4],D,ans; void exgcd(ll a,ll b,ll &d,ll &x,ll &y){ if(!b){ d=a; x=1; y=0; } else{ exgcd(b,a%b,d,y,x); y-=x*(a/b); } } void CRT(int r){ ll M=1,Mi,x0,y0,d; ans=0; for(int i=1;i<=r;++i){ M*=m[i]; } for(int i=1;i<=r;++i){ Mi=M/m[i]; exgcd(Mi,m[i],d,x0,y0); ans=(ans+Mi*x0*a[i])%M; } if(ans<=D){ ans+=M; } } int main(){ // freopen("poj1006.in","r",stdin); int zu=0; m[1]=23; m[2]=28; m[3]=33; while(1){ ++zu; scanf("%lld%lld%lld%lld",&a[1],&a[2],&a[3],&D); if(a[1]==-1 && a[2]==-1 && a[3]==-1 && D==-1){ break; } CRT(3); printf("Case %d: the next triple peak occurs in %lld days.\n",zu,ans-D); } return 0; }
【数论】【中国剩余定理】poj1006 生理周期
声明:以上内容来自用户投稿及互联网公开渠道收集整理发布,本网站不拥有所有权,未作人工编辑处理,也不承担相关法律责任,若内容有误或涉及侵权可进行投诉: 投诉/举报 工作人员会在5个工作日内联系你,一经查实,本站将立刻删除涉嫌侵权内容。