这题拿来以后以为是个神题.....T T
后来看了结题报告,说这是简单题,报告上说啥我没有理解。后来又看了看其他的东西,发现我有一个公式不知道,这个公式详见:http://hi.baidu.com/aekdycoin/item/e493adc9a7c0870bad092fd9
A^x=A^(x%phi(p)+phi(p)(mod p)...其中,x>=phi(m)
这样就可以通过公式求解,当x<phi(p)时,暴力求解。当x>=phi(p)时,把A用A%p替换,这样还是枚举暴力。。。
当p=1时,结果为m+1,由于m=2^64-1,所以当m为2^64-1时,要打印个字符串表示2^64,否则溢出,坑爹。。。。
另外,我当时做的时候WA了n多次,一个很2B的错误就是要用unsigned Long long......long long 的范围是2^63-1....
还有,我当时在处理p=2时搞了好半天。。。因为这个时候phi(p)=1,代码写不好的话,很容易杯具,样例给的很好,要把内个样例真正想明白才行,附代码:
#include <iostream>
#include <vector>
#ifdef _WIN32
#define LL unsigned __int64
#define out64 "%I64u"
#define in64 "%I64u"
#else
#define LL unsigned long long
#define out64 "%llu"
#define in64 "%llu"
#endif
#define cl(a) memset(a,0,sizeof(a))
#define ss(a) scanf("%d",&a)
#define sl(a) scanf(in64,&a)
#define pb push_back
using namespace std;
const int N=100100;
int a[N],c[N],e[N],p,b;
LL m;
void prime()
{
int i,j;
for (i=2;i<N;i++)
if (!a[i])
{
for (j=i*2;j<N;j+=i)
{
if (!a[j]) c[j]=i;
a[j]=1;
}
}
}
void euler()
{
int i,k;
for (i=2;i<N;i++)
if (!a[i]) e[i]=i-1;
else
{
k=i/c[i];
if (k%c[i]==0) e[i]=c[i]*e[k];
else e[i]=(c[i]-1)*e[k];
}
}
LL js(int x,int y)
{
LL k=y,r=x,res=1;
while (k>0)
{
if (k&1) res=(res*r)%p;
r=(r*r)%p;
k>>=1;
}
return res;
}
void print(int cas,LL z)
{
printf("Case #%d: "out64"\n",cas,z);
}
int main()
{
int T,i,cas;
ss(T);
prime();
euler();
for (cas=1;cas<=T;cas++)
{
ss(b);ss(p);sl(m);
if (p==1)
{
if (m==18446744073709551615ULL)
printf("Case #%d: 18446744073709551616\n",cas);
else print(cas,m+1);
continue;
}
LL n=0,r=1,z=0;
while (r<e[p]&&n<=m)
{
if (js(n,r)==b) z++;
n++;
r*=n;
}
r=r%e[p];
while (r!=0&&n<=m)
{
if (js(n,e[p]+r)==b) z++;
n++;
r=(r*n)%e[p];
}
if (n>m)
{
print(cas,z);
continue;
}
if (n>0) n--;
else if (b==0) z++;
for (i=0;i<p;i++)
if (js(i,e[p])==b)
{
LL k1=0,k2=0;
if (m>=i) k1=(m-i)/p+1;
if (n>=i) k2=(n-i)/p+1;
z=z+(k1-k2);
}
print(cas,z);
}
return 0;
}