PJU3696-The Luckiest number

#include<stdio.h>
#include<string.h>

#define LL long long

LL gcd(LL a,LL b)
{
	return b==0?a:gcd(b,a%b);
}

//φ(n)
LL eular(LL n)
{
	LL ans = n;
	for (LL i=2; i*i<=n; i++)
		if (n%i==0)
		{
			ans -= ans/i;
			while(n%i==0) n/=i;
		}
	if (n>1) ans -= ans/n;
	return ans;
}

//a*b%p
LL p;
LL multi(LL a,LL b)
{
	LL ans = 0;
	while(b)
	{
		if (b&1) ans = (ans+a)%p;
		a = a*2%p;
		b >>= 1;
	}
	return ans;
}

//a^b%p
LL power_mod(LL a,LL b)
{
	LL ans=1;
	while(b)
	{
		if (b&1) ans = multi(ans,a);
		a = multi(a,a);
		b >>= 1;
	}
	return ans;
}

int main()
{
	int iCase = 1;
	int L;
	while (scanf("%d",&L)!=EOF && L)
	{
		p = L/gcd(L,8)*9;
		if (gcd(10,p)!=1 || L%5==0)
		{
			printf("Case %d: 0\n",iCase++);
			continue;
		}
		LL n = eular(p);
		LL number = n;
		for (LL i=2; i*i<n; i++)
			if (n%i==0)
			{
				int cnt = 0;
				while (n%i==0)
				{
					cnt++;
					n /= i;
				}
				while (cnt)
				{
					if (power_mod(10,number/i)%p==1)
					{
						number /= i;
						cnt--;
					}
					else break;
				}
			}
		if(n>1 && power_mod(10,number/n)%p==1) number /= n;
		printf("Case %d: %lld\n",iCase++,number);
	}
	return 0;
}
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值