POJ 3696 The Luckiest number（算法竞赛进阶指南，欧拉函数，快速幂）

算法竞赛进阶指南，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
*/