CF_603B(欧几里得算法_规律_数论)

打了一下欧几里得算法,又wrong了一次才找到规律;

注意k==1时特判;

找到规律后又错了好多次,好郁闷(;′⌒`)(;′⌒`)(;′⌒`)

#include<stdio.h>
#include<math.h>
#include<stdlib.h>
#include<algorithm>
#include<string.h>
using namespace std;
long long x,y;
long long res[1000000+7];
long long vis[1000000+7];
long long pow4(long long a,long long b)
{
	long long  r=1,base=a;
	while(b!=0)
	{
		if(b&1)
		r*=base;
		r=r%(1000000000+7);
		base*=base;
		base=base%(1000000000+7);
		b>>=1;
	}
	return r;
}
long long  extgcd(long long a,long long b)
{
    if(b==0)
    {
        x=1;
        y=0;
        return a;
    }
    long long m=extgcd(b,a%b);
    long long tmp=x;
    x=y;
    y=tmp-a/b*y;
    return m;
}
int main()
{
   long long a,b;
   while(scanf("%I64d%I64d",&a,&b)!=EOF)
   {
   	memset(res,0,sizeof(res));
   	memset(vis,0,sizeof(vis));
    int zige=1;
    long long n;
    long long aa=a,bb=b;
    zan=a;
    a=b;
    b=-zan;
    for(n=0;n<=aa-1;n++)
    {
        long long e=extgcd(a,b);
        {
        	res[n]=(x*n/e%(b/e)+(b/e))%(b/e);
        	while(res[n]<0)
            {
                res[n]+=(b/e)<0?-(b/e):(b/e);
            }
        }
    }
   else
   {
   	  long long t=1;
   	  for(long long i=1;i<=aa-1;i++)
   	  {
        if(vis[i])
        {
            continue;
        }
   	  	if(res[i]==i)
   	  	{
   	  	    t+=1;
   	  		continue;
   	  	}
   	  	else
   	  	{
   	  		long long q=i;
   	  		for(;;)
   	  		{
   	  			vis[q]=t;
   	  			q=res[q];
   	  			if(vis[q]!=0)
   	  			{
   	  				t+=1;
					 break;
				}
   	  		}
   	  	}
   	  }
   	  if(bb!=1)
   	  printf("%I64d\n",pow4(aa,t-1));
   	  else printf("%I64d\n",pow4(aa,aa));
   	  continue;
   }
   }
}

附上百度百科模运算性质

基本概念

给定一个 正整数
   
,任意一个整数
   
,一定存在等式
   
其中
   
   
是整数,且
   
,称
   
   
除以
   
的商,
   
   
除以
   
余数
 
对于 正整数和整数
   
,
   
,定义如下运算:
取模运算:a % p(或a mod p),表示a除以p的 余数
模p 加法:(a + b) % p ,其结果是a+b算术和除以p的余数,也就是说,(a+b) = kp +r,则(a + b) % p = r。
模p减法:(a-b) % p ,其结果是a-b算术差除以p的余数。
模p 乘法:(a * b) % p,其结果是 a * b算术乘法除以p的余数。
说明:
1. 同余式:正整数a,b对p取模,它们的余数相同,记做 a ≡ b % p或者a ≡ b (mod p)。
2. n % p得到结果的正负由 被除数n决定,与p无关。例如:7%4 = 3, -7%4 = -3, 7%-4 = 3, -7%-4 = -3。

基本性质

(1)若p|(a-b),则a≡b (% p)。例如 11 ≡ 4 (% 7), 18 ≡ 4(% 7)
(2)(a % p)=(b % p)意味a≡b (% p)
(3) 对称性:a≡b (% p)等价于b≡a (% p)
(4)传递性:若a≡b (% p)且b≡c (% p) ,则a≡c (% p)

运算规则

模运算与基本四则运算有些相似,但是除法例外。其规则如下:
(a + b) % p = (a % p + b % p) % p (1)
(a - b) % p = (a % p - b % p) % p (2)
(a * b) % p = (a % p * b % p) % p (3)
(a^b) % p = ((a % p)^b) % p (4)
结合律:
((a+b) % p + c) % p = (a + (b+c) % p) % p (5)
((a*b) % p * c)% p = (a * b*c) % p (6)// (a%p*b)%p=(a*b)%p
交换律:
(a + b) % p = (b+a) % p (7)
(a * b) % p = (b * a) % p (8)
分配律:
((a +b)% p * c) % p = ((a * c) % p + (b * c) % p) % p (9)
重要定理:
若a≡b (% p),则对于任意的c,都有(a + c) ≡ (b + c) (%p);(10)
若a≡b (% p),则对于任意的c,都有(a * c) ≡ (b * c) (%p);(11)
若a≡b (% p),c≡d (% p),则 (a + c) ≡ (b + d) (%p),(a - c) ≡ (b - d) (%p),
(a * c) ≡ (b * d) (%p),(a / c) ≡ (b / d) (%p); (12)

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值