不知道亚洲赛的风格是不是这样,我感觉这题无法形容~~~
为了真正理解这题,所以我决定好好写个解题报告。
我是看了后来发的解题报告写的,按他所说,我们分两种情况考虑。
第一种:当n比较大的时候,就是达到n!%phi(p)=0的时候,其实就是求n^phi(p)%p,显而易见的。
那么比n 小的那部分就直接暴力。
第二种:当n比较小的时候,我们也是暴力了。
这两种其实也就合二为一了,因为都需要暴力的。
#include<iostream>
#include<cstdlib>
#include<stdio.h>
using namespace std;
typedef unsigned __int64 ll;
ll getphi(ll n)//求欧拉函数值
{
ll ans=1,i;
for(i=2;i*i<=n;i++)
{
if(n%i==0)
{
n/=i;
ans*=i-1;
while(n%i==0)
{
n/=i;ans*=i;
}
}
}
if(n>1) ans*=n-1;
return ans;
}
int powmod(ll a,int b,int c)//a^b%c
{
int ans=1;
while(b)
{
if(b&1) ans=ans*a%c;
b>>=1;
a=a*a%c;
}
return ans;
}
int main()
{
int t,b,p;
int count=1;
ll m;
scanf("%d",&t);
while(t--)
{
cin>>b>>p>>m;
int phi=getphi(p);//phi(p);
ll mul=1;int num;ll ans=0;int add=0;
for(int i=0;i<=p;i++)
{
mul*=i?i:1;
if(mul>=phi) add=phi;
mul%=phi;
if(!mul){num=i;break;}//找到n!%p==0的临界值了
if(i<=m)//
if(powmod(i,mul+add,p)==b) ans++;
}
for(int i=0;i<p&&i<=m;i++)
{
if(powmod(i,phi,p)==b)
{
ans+=(m-i)/p+1;//这里想一下就知道了余数循环
if(num>=i+1)//这个地方还真容易忽视。。。
ans-=(num-1-i)/p+1;
}
}
if(m== 18446744073709551615LL&&p == 1&&b == 0)
{
printf("Case #%d: 18446744073709551616\n", count++);
}
else
printf("Case #%d: %I64u\n", count++, ans);
}
return 0;
}