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

最新推荐文章于 2021-08-19 19:06:55 发布

iewlhf

最新推荐文章于 2021-08-19 19:06:55 发布

阅读量361

点赞数

分类专栏：数学

本文链接：https://blog.csdn.net/iewlhf/article/details/51628995

版权

数学专栏收录该内容

1 篇文章 0 订阅

订阅专栏

打了一下欧几里得算法，又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;
   }
   }
}

附上百度百科模运算性质