码迷,mamicode.com
首页 > 其他好文 > 详细

数论之中国剩余定理

时间:2019-01-24 00:18:56      阅读:182      评论:0      收藏:0      [点我收藏+]

标签:print   12px   res   space   amp   bsp   不定方程   同余   splay   

中国剩余定理是中国古代求解一次同余式组(见同余)的方法,是数论中一个重要定理,又称孙子定理。

孙子定理

  孙子定理主要用来求解模线性不定方程组。

  对于方程组:

         x0≡a0(mod b0)

       x1≡a1(mod b1)

       x2≡a2(mod b2) 

       .  .  .  

       xn≡an(mod bn)

  设m=b0*b1*b2*...*bn,mi=m/bi

  得到方程组

      mixi+biyi=1

  扩欧求解后,原模线性方程组的解为:

      ∑mi*ai*xi  mod m

看代码吧:

技术分享图片
 1 #include <cstdio>
 2 using namespace std;
 3 const int MAXN=10050;
 4 int n;
 5 int a[MAXN],b[MAXN];
 6 
 7 void exgcd(long long a,long long b,long long &x,long long &y)
 8 {
 9     if(a%b==0)
10     {
11         x=0;y=1;return ;
12     }
13     exgcd(b,a%b,x,y);
14     long long t=x;
15     x=y;y=t-a/b*y;
16 }
17 
18 long long China_Reminder(int n)
19 {
20     long long m=1,mi,res=0,x,y;
21     for(int i=1;i<=n;i++)m*=b[i];
22     for(int i=1;i<=n;i++)
23     {
24         mi=m/b[i];
25         exgcd(mi,b[i],x,y);
26         x=(x%b[i]+b[i])%b[i];
27         res=(res+mi*a[i]*x%m)%m;
28     }
29     return (res+m)%m;
30 }
31 
32 int main()
33 {
34     scanf("%d",&n);
35     for(int i=1;i<=n;i++)scanf("%d %d",&a[i],&b[i]);
36     printf("%lld\n",China_Reminder(n));
37     return 0;
38 }
View Code

例题:fzu1402 猪的安家

数论之中国剩余定理

标签:print   12px   res   space   amp   bsp   不定方程   同余   splay   

原文地址:https://www.cnblogs.com/lzxzy-blog/p/10311968.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!