一、题目
二、解法
首先要理解这个 a a a天,一个 f 0 f_0 f0有贡献的话一定是 1 1 1和这个点有一条路径,那么我们只需要考虑这个长度为 a a a的路径个数奇偶就行了,奇数的话计入贡献。现在就不难想到矩阵快速幂,求出邻接矩阵的 a a a次方即可。
这样做的时间复杂度是 O ( n 3 q log a ) O(n^3q\log a) O(n3qloga),就算用 b i t s e t bitset bitset优化也要 T \text{T} T,考虑我们只需要知道 1 1 1到各个点的情况,那么我们预处理邻接矩阵的 2 2 2幂次次方,每次用第一行去乘,相当于一个倍增,不难得到时间复杂度 O ( n 3 log a ) + O ( q n 2 log a ) O(n^3\log a)+O(qn^2\log a) O(n3loga)+O(qn2loga)
#include <cstdio>
#include <cmath>
#include <map>
using namespace std;
#define int long long
int read()
{
int num=0,flag=1;char c;
while((c=getchar())<'0'||c>'9')if(c=='-')flag=-1;
while(c>='0'&&c<='9')num=(num<<3)+(num<<1)+(c^48),c=getchar();
return num*flag;
}
int n,m,k;map<int,int> mp;
int qkpow(int a,int b,int p)
{
int r=1;
while(b>0)
{
if(b&1) r=r*a%p;
a=a*a%p;
b>>=1;
}
return r;
}
int gcd(int a,int b)
{
return !b?a:gcd(b,a%b);
}
bool bsgs(int y,int z,int p)
{
int k=0,a=1;
while(1)
{
int d=gcd(y,p);if(d==1) break;
if(z%d) return 0;
z/=d;p/=d;++k;a=a*y/d%p;
if(z==a) {printf("%lld\n",k);return 1;}
}
mp.clear();int m=sqrt(p)+1;
for(int i=0,t=z%p;i<m;i++,t=t*y%p) mp[t]=i;
for(int i=1,tt=qkpow(y,m,p),t=a*tt%p;i<=m;i++,t=t*tt%p)
{
if(!mp.count(t)) continue;
printf("%lld\n",i*m-mp[t]+k);return 1;
}
return 0;
}
signed main()
{
while(~scanf("%lld %lld %lld",&n,&m,&k) && n+m+k)
{
if(!bsgs(n,k,m)) puts("No Solution");
}
}