69久久精品无码一区二区,91精品国产高清久久久久久91,久久综合久久自在自线精品自

原语饿狼 — Sat, 05 Sep 2009 11:49:00 GMT

求解模线性方�E?nbsp; axºb(mod n)

1.必备知识�Q�扩展欧几里得算法的知识�Q�可查看我的数论(3)------扩展�Ƨ几里得��法

2.基本思�\�Q?br>
   设d=gcd(a,n)�Q�用扩展�Ƨ几里得��法解线性方�E?ax'+ny'=d.
     如果d|b,则方�E�axºb(mod n)有一个解的值x0=x'(b/d)mod n

   ��法��D��里说�Q�（�q�没理解�Q?br>   方程axºb(mod n)有解�Q�即存在d|b,其中d=gcd(a,n)),x0是该方程的�Q意一个解�Q�则该方�E�对模n恰有d个不同的
   解，分别�?x(i)=x(0)+i(n/d)(i=1,2,...d).
   特别的设e=x0+n�Q�方�E�ax=b(mod n)的最��整数解x1=e mod (n/d),最大整数解x2=x1+(d-1)*(n/d)�?br>

   所以实际上用欧几里得算法记得x'��可以知道结果了�?br>
3.源代码模�?br>

1//扩展�Ƨ几里得��法
2int Extended_Euclid(int a,int b,int& x,int &y)
3{
4    if(b==0){
5        x=1;
6        y=0;
7        return a;
8    }
9    int d=Extended_Euclid(b,a%b,x,y);
10    int temp=x;x=y;y=temp-a/b*y;
11    return d;
12}

1//用扩展欧几里得解模线性方�E�ax=b (mod n)
2bool modularLinearEquation(int a,int b,int n)
3{
4    int x,y,x0,i;
5    int d=Extended_Euclid(a,n,x,y);
6    if(b%d)
7           return false;
8        x0=x*(b/d)%n;
9    for(i=1;i<=d;i++)
10       printf("%d\n",(x0+i*(n/d))%n);
11    return true;
12}
13

原语饿狼 2009-09-05 19:49 发表评论

数论(3)-------�Ƨ拉phi函数

原语饿狼 — Sat, 05 Sep 2009 10:51:00 GMT

�Ƨ拉phi函数

1.应用�Q?br> 对一个正整数n,求小于n且与n互质�Q�包�?�Q�的个数�?br>
2.公式�Q?br>

         I
Φ(n) =n ∏ (1 - 1 / pi)�Q�其中pi表示n的质因子�Q?br>        i=1

    I
n = ∏ (Pi)^ki(I �?n 的素因子的个�?
   i=1

如：Φ�Q?0�Q?10�Q?-1/2)(1-1/5)=4,其中2,5�?0的质因子.

3.证明�Q?br>
   要证明这个式子，我们先来看几个基本的公式�?br>
   �Q?�Q?nbsp; Φ(p)=p-1,p是质�?br>   �Q?�Q?nbsp; Φ(p*q)=Φ(p)Φ(q) p,q是质�?br>
         Φ(p*q)=p*q-1- (q-1)(注：【p,2p,...(q-1))】个数q-1�Q?-(p-1)(注：【q,2q,...(p-1)q】个数p-1)
              =(p-1)(q-1)=Φ(p)Φ(q)

    (3) 对于整数n,n=p^k
               φ(n) = p^k - p^k-1
         ��于 pk 的正整数个数�?pk - 1个，其中
         �?pk 不互质的正整数有{p * 1,p * 2,...,p * (p^k-1-1)} ��p�� p^k-1- 1 �?br>         所�?φ(n) = p^k -1 ( p^k-1- 1) = p^k-p^k-1 �?br>

     接下来要证明上面那个�Ƨ拉函数��是轻而易丄��事情了�?br>

            I
Φ(n) = Φ( ∏ (Pi)^ki₎
           i=1
          I
     = ∏[(Pi)^ki- (Pi)^ki-1 ]
      i=1
再除以n��可以求�?font face=黑体>Φ(n)�?/font>

4.源代码模�?/span>

1int phi(int n)
2{
3    int ans,i,k;
4    if(n==1)
5        ans=0;
6
7    else{
8      ans=n;
9      k=1;
10      for(i=2;n!=1;i+=k){
11        if(n%i==0){
12            ans*=(i-1);ans/=i;
13            while(n%i==0) n/=i;
14            i=k;
15        }
16
17      }
18    }
19    return ans;
20}

原语饿狼 2009-09-05 18:51 发表评论

数论(2)-------扩展�Ƨ几里得��法

原语饿狼 — Sat, 05 Sep 2009 08:51:00 GMT

扩展�Ƨ几里得��法------求解�U�性方�E�ax+by=c

1.应用�Q?br>      �U�性方�E�ax+by=c ,已知a,b,c�Q�求解x,y.

2.基本思�\�Q?br>
      ax+by=c有解 => c=k*gcd(a,b)=kd�Q�因为d=gcd(a,b)=>d|(ax+by))

      我们先考虑求解 ax+by=d
        由欧几里得算�?d=bx'+(a mod b)y'=bx'+(a-[a/b]b)y'=ay'+b(x'-[a/b])y'
        则由上述两式子，我们可以得出 x=y' ,y=x'-[a/b]y'
        �q�样子，在欧几里得算法添加x,y变量�Q�最后得到解。（可结合下面代码源代码�q�行理解�Q?br>
      接下来我们来看看ax'+by'=d和ax+by=c之间的关�p?br>
        �Q�c/d)ax'+(c/d)by'=(c/d)d �?nbsp;可以得到 x=(c/d)x',y=(c/d)y'

         所以可以得到ax+by的一�l�解
         那么ax+by=c所有解的�Ş式是什么呢�Q?br>             a(x+qb)+b(y-qa)=c; q��Z�Q意整�?br>(注意�Q�当要求y-qa的最��正整数min�?由y-qa>=0�Q?q取[y/a]最��，min=y-[y/a]y,但是,[y/a]可能�?�Q�如果y是负敎ͼ�min此时也�ؓ负数�Q�不好，此时令min+=a��可以取得最��正整数��g��Q�[y/a]=0所以|y|
3.源代码模�?br>

1int Extended_Euclid(int a,int b,int& x,int &y)
2{
3    if(b==0){
4        x=1;
5        y=0;
6        return a;
7    }
8    int d=Extended_Euclid(b,a%b,x,y);
9    int temp=x;x=y;y=temp-a/b*y;
10    return d;
11}
12//用扩展欧几里得算法解�U�性方�E�ax+by=c;
13bool linearEquation(int a,int b,int c,int& x,int &y)
14{
15    int d=Extended_Euclid(a,b,x,y);
16    if(c%d) return false;
17
18    int k=c/d;
19    x*=k;y*=k;//求的只是其中一个解
20    return true;
21}

�?学习心得
扩展�Ƨ几里得��法的妙用，值得好好体会�Q�对整个�q�程自己可以找个例子好好印证一下，肯定印象深刻。去Google一下POJ上的�q�方面的题目好好实践一下吧�Q�！�Q?/span>

原语饿狼 2009-09-05 16:51 发表评论

数论(1)-----�Ƨ几里得��法

原语饿狼 — Sat, 05 Sep 2009 07:48:00 GMT

一. �Ƨ几里得��法------求最大公�U�数

1.公式�Q?br>     gcd(a,b)=gcd(b,a mod b)    (a为非负整敎ͼ�b为正整数)

2.证明�Q?br>      思�\�Q?br>           两个整数a和b�Q�如果a|b&&b|a�Q�即a,b能互相整除）�Q�那么a=b.

      基础知识准备�Q?br>           A. (a mod b)=a-qb , q=(int)a/b;
           B. d|a&&d|b => d|(xa+yb) x,y��Z�Q意整�?br>           C. d|a&&d|b => d|gcd(a,b)

      �q�程�Q�（1�Q�证�Q�gcd(a,b)|gcd(b,a mod b)

                  设d=gcd(a,b)=>d|a&&d|b,
                  由A和B知道�Q�d|a&&d|b=>d|(xa+yb)=>d|(a mod b)
                  由C知道�Q�d|b&&d|(a mod b)=>d|gcd(b,a mod b)=>gcd(a,b)|gcd(b,a mod b);

            (2) 证：gcd(b,a mod b)|gcd(a,b)

                  设d=gcd(b,a mod b)=>d|b&&d|(a mod b),
                  由A和B知道, d|b&&d|(a mod b)=>d|(xb+y(a mod b)=>d|a(由A�Q�a=qb+(a mod b))
                  由C知道�Q�d|a&&d|b=>d|gcd(a,b)=>gcd(b,a mod b)|gcd(a,b)

             由（1�Q�和�Q?�Q�，可以知道我们得证了�?br>
3.�E�序模板�Q?

//递归版本
int gcd(int a,int b)
{
    return b?gcd(b,a%b):a;
}

//循环版本
int gcd1(int a,int b)
{
  for(int c=a%b;c;a=b,b=c,c=a%b);
  return b;
}

4.学习心得

    �Ƨ几里得��法是之后很多数论算法的基础�Q�了解它的原理是很有必要的�?br>    自己要�D几个例子来熟悉一下算法的执行�q�程中的每一步骤�Q�这��h��能记忆深刅R�?br>    上述的基��知识也是很有用的�Q��^时注意积累。不懂的地方��几个实例看一下�?nbsp;

原语饿狼 2009-09-05 15:48 发表评论

69久久精品无码一区二区,91精品国产高清久久久久久91,久久综合久久自在自线精品自

数论(3)-------�Ƨ拉phi函数

数论(2)-------扩展�Ƨ几里得���法

数论(1)-----�Ƨ几里得���法

数论(2)-------扩展�Ƨ几里得��法

数论(1)-----�Ƨ几里得��法