xdu_1077:循环节长度

题意很简单，就是给出p,q，求p/q的循环节长度。

由循环小数的循环部分的值等于等比数列求和的值S，列公式得到最简分数分母的值。最终得10^x%q==1（其中q为经过modify之后的值）。搞清这些之后，套kuangbin大神的BSGS模板就可以了。

题目链接： http://202.117.120.247/problem.php?id=1077

 #include<bits/stdc++.h>
 using namespace std;
 //baby_step giant_step
 // a^x = b (mod n) n为素数，a,b < n
 // 求解上式 0<=x < n的解
 #define MOD 76543
 #define next nxt
 int hs[MOD],head[MOD],next[MOD],id[MOD],top;
 void insert(int x,int y)
 {
     int k = x%MOD;
     hs[top] = x, id[top] = y, next[top] = head[k], head[k] = top++;
 }
 int find(int x)
 {
     int k = x%MOD;
     ; i = next[i])
         if(hs[i] == x)
             return id[i];
     ;
 }
 int gcd(int a,int b)
 {
     return b? gcd(b,a%b):a;
 }
 int BSGS(int a,int b,int n)
 {
     memset(head,-,sizeof(head));
     top = ;
     ) ;
     int m = sqrt(n*1.0), j;
     , p = ;
     ; i < m; ++i, p = p*a%n)insert(p*b%n,i);
     for(long long i = m; ; i += m)
     {
          && i != j )
         {
 //            cout<<i<<"=="<<j<<endl;
             return i-j;
         }
         if(i > n)break;
     }
     ;
 }
 int main()
 {
     //freopen("in.txt","r",stdin);
     //freopen("out.txt","w",stdout);
     ,b=,n;
     int p,q;
     while(~scanf("%d%d",&p,&q))
     {
         n=q/gcd(p,q);
 //        while(n%2==0) n/=2;
 //        while(n%5==0) n/=5;
         int ans = BSGS(a,b,n);
         printf("%d\n",ans);
     }
 }