51nod 1035 最长的循环节 数学

1035 最长的循环节

题目连接：

https://www.51nod.com/onlineJudge/questionCode.html#!problemId=1035

Description

正整数k的倒数1/k，写为10进制的小数如果为无限循环小数，则存在一个循环节，求<=n的数中，倒数循环节长度最长的那个数。

1/6= 0.1(6) 循环节长度为1

1/7= 0.(142857) 循环节长度为6

1/9= 0.(1) 循环节长度为1

Input

输入n（10 <= n <= 1000)

Output

输出<=n的数中倒数循环节长度最长的那个数

Sample Input

10

Sample Output

7

Hint

题意

题解：

推荐一篇科普性的论文：http://w3.math.sinica.edu.tw/math_media/d253/25311.pdf

补充一下，这道题最后要求解一个 10x≡1(modC) 这个玩意儿

如果gcd(10,C)!=1的话，显然无解。

如果存在解的话，根据欧拉公式，那么这个解 x|phi(C)，所以直接暴力枚举x就好了。

代码

#include<bits/stdc++.h>
using namespace std;
const int maxn = 1005;
long long phi[maxn];
int n;
int gcd(int a,int b)
{
if(b==0)return a;
return gcd(b,a%b);
}
long long quickpow(long long  m,long long n,long long k)//返回m^n%k
{
long long b = 1;
while (n > 0)
{
if (n & 1)
b = (b*m)%k;
n = n >> 1 ;
m = (m*m)%k;
}
return b;
}
void phi1()
{
memset(phi,0,sizeof(phi));
phi[1]=1;
for(long long i=2;i<=n;i++)
{
if(!phi[i])
{
for(long long j=i;j<=n;j+=i)
{
if(!phi[j]) phi[j]=j;
phi[j]=phi[j]/i*(i-1);
}
}
}
}
int solve(int x)
{
for(int i=1;i*i<=phi[x];i++)
{
if(phi[x]%i==0)
{
if(quickpow(10,i,x)==1)
return i;
}
}
for(int i=sqrt(phi[x]);i>=1;i--)
{
if(phi[x]%i==0)
{
if(quickpow(10,phi[x]/i,x)==1)
return phi[x]/i;
}
}
return -1;
}
int main()
{
scanf("%d",&n);
phi1();
int ans1 = 0,ans2 = 0;
for(int i=1;i<=n;i++)
{
int C = solve(i);
if(ans2<C)
ans2 = C,ans1 = i;
}
cout<<ans1<<endl;
}