算法竞赛进阶指南,150页,欧拉函数,快速幂
题目意思:
题目给定一个正整数 L, 问有至少多少个连续的8组成的正整数是 L 的倍数。
书上的思路:
1、 888888 这个式子来表示:k * (10^x - 1) / 9
进而简化:8 * (10^x - 1) / 9 = L * k (k是一个整数)
8 * (10^x - 1) = 9L * k, 假设 d = gcd(8, L), 那么 (10^x - 1) = k * (9L / d), 从而
10^x = 1(mod 9L / d)
书上的引理:
正整数 a 和 n, 如果a,n互质, 则 a^x = 1(mod n) 的最小正整数x0 是 n的欧拉函数 phi(n) 的约数。
通过反证法证明。
2、 因此,求出 9L / d 的欧拉函数 ph, 扫描ph的所有约数x,看看 式子 10^x = 1(mod 9L / d) 是否满足即可。
选出满足条件的最小约数。
3、无解的时候就是 9L / d 与10不互质的时候;
4、 写快速幂的时候,不要直接相乘,会溢出。
#include <cstdio>
#include <cstring>
#include <iostream>
#include <cmath>
using namespace std;
long long L;
long long gcd(long long a, long long b)
{
return b? gcd(b, a % b) : a;
}
long long phi(long long n) //欧拉函数
{
long long ans = n;
for (long long i = 2; i * i <= n; i++)
if (n % i == 0)
{
ans = ans / i * (i - 1);
while (n % i == 0)
n /= i;
}
if (n > 1)
ans = ans / n * (n - 1);
return ans;
}
long long multi(long long a, long long b, long long mod)
{
long long ret = 0;
while(b)
{
if(b & 1)
{
ret = (ret + a) % mod;
}
a = (a << 1) % mod;
b >>= 1;
}
return ret;
}
long long quick_pow(long long a, long long b, long long mod)
{
long long ans = 1;
while(b > 0)
{
if(b & 1)
{
ans = multi(ans, a, mod); //这里直接相乘会溢出
}
b >>= 1;
a = multi(a, a, mod);
}
return ans;
}
void solve()
{
long long a = 10 , d = gcd(L, 8), n = 9 * L / d;
if(gcd(a, n) != 1)
{
printf("0\n");
return;
}
long long ph = phi(n), sq = sqrt(ph);
// 枚举 n的欧拉函数 phi(n) 的所有约数
for(long long i = 1; i <= sq; ++i) //求约数的时候,先扫描到 sqrt(ph), 全扫的话超时
{
if(ph % i == 0 && quick_pow(a, i, n) == 1)
{
printf("%lld\n", i);
return;
}
}
for(long long i = sq - 1; i; --i)
{
if(ph % i == 0 && quick_pow(a, ph / i, n) == 1)
{
printf("%lld\n", ph / i);
return;
}
}
printf("0\n");
}
int main()
{
int cnt = 0;
while(scanf("%lld", &L) != EOF && L)
{
printf("Case %d: ", ++cnt);
solve();
}
return 0;
}
/*
8
11
16
0
*/
/*
Case 1: 1
Case 2: 2
Case 3: 0
*/