[NOI Online 3 提高组]魔法值

一、题目

点此看题

二、解法

首先要理解这个$a$天，一个$f_0$有贡献的话一定是$1$和这个点有一条路径，那么我们只需要考虑这个长度为$a$的路径个数奇偶就行了，奇数的话计入贡献。现在就不难想到矩阵快速幂，求出邻接矩阵的$a$次方即可。

这样做的时间复杂度是$O(n^{3}q\log a)$，就算用$bitset$优化也要$\text{T}$，考虑我们只需要知道$1$到各个点的情况，那么我们预处理邻接矩阵的$2$幂次次方，每次用第一行去乘，相当于一个倍增，不难得到时间复杂度$O(n^3\log a)+O(q{n}^2\log a)$

#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");
	}
}