Period of an Infinite Binary Expansion-----欧拉函数，欧拉定理

来源：互联网发布：linux coredump 编辑：程序博客网时间：2024/05/17 04:05

题目：http://poj.org/problem?id=3358

这一题很经典。

题解：源于http://www.cnblogs.com/ACKOKO/articles/2119216.html

这里讲一下2^i*(2^(j-i)-1)%q==0，为什么将q中的2的个数是i，因为2^i中的i要尽量小(i是循环的前一位，题目要求i要小)，并且要满足去除q中所有的2，所以i要等于q中的2的个数。

源代码：

#include <stdio.h>#include <algorithm>using namespace std;typedef long long ll;int cnt;int m,n,d;int factor[1000000];int  gcd(int a,int b){    if(b==0)  return a;    return gcd(b,a%b);}int euler(int m){    int s=m;    for(int i=2;i*i<=m;i++)     {         if(m%i==0)         {             while(m%i==0) {m=m/i;}             s=s/i*(i-1);         }     }     if(m>1)     s=s/m*(m-1);     return s;}ll quickPow(int b,int m){   ll  s=1,a=2;   while(b>0)   {     if(b%2==1) s=s*a%m;     a=a*a%m;     b=b/2;   }   return s;}int solve(int  m){    int flag=0,y,k=0;    int sum,x=euler(m);     if(x==1)  return 1;     for(int i=2;i*i<=x;i++)      if(x%i==0)      {          factor[k++]=i;          factor[k++]=x/i;      }      sort(factor,factor+k);      for(int i=0;i<k;i++)          {              sum=quickPow(factor[i],m);              if(sum==1)              {                  y=factor[i];                 flag=1;                break;              }           }       if(!flag) y=x;      return y; }int   main(){   int  cas=1;   while(scanf("%d/%d",&n,&m)!=EOF)  {      cnt=1;      d=gcd(n,m);      n=n/d;  m=m/d;      while(m%2==0)      {          m=m/2;          cnt++;      }      printf("Case #%d: %d,%d\n",cas++,cnt,solve(m));  }   return 0;}