BZOJ 1898: [Zjoi2005]Swamp 沼泽鳄鱼

邻接矩阵求路线数

特殊性在于周期，不过求下LCM还是可以搞的。

可以发现图是周期变化的，以12为一周期，于是建13个邻接矩阵，第13个为前12个的乘积，于是有G(K)=(G(13))^(k/12)*(G(1)*G(2)*……*G(k%12))，O(n^3*logk)

话说没注意邻接矩阵的含义WA了一次是闹哪样。。。。。

#include<iostream>
#include<cstdio>
#include<cstring>
using namespace std;
typedef long long ll;
const int N=50+5;
const int p=1e4;
int n;
struct matrix{
	ll a[N][N];
	matrix(){memset(a,0,sizeof(a));}
	void clear(){memset(a,0,sizeof(a));}
	matrix operator * (const matrix &b)const{
		static matrix ans;
		ans.clear();
		for(int i=0;i<n;i++)
		for(int j=0;j<n;j++)
		for(int k=0;k<n;k++)
		ans.a[i][j]+=a[i][k]*b.a[k][j];
		for(int i=0;i<n;i++)
		for(int j=0;j<n;j++)
		ans.a[i][j]%=p;
		return ans;
	}
	matrix operator ^ (int k)const{
		static matrix ans,b;
		ans.clear();
		for(int i=0;i<n;i++)ans.a[i][i]=1;
		memcpy(b.a,a,sizeof(a));
		for(;k;k>>=1,b=b*b)if(k&1)ans=ans*b;
		return ans;
	}
}g[13];
bool can[13][N];
int u[N*N],v[N*N],pi[10];
int main(){
	//freopen("a.in","r",stdin);
	int m,s,t,k;
	scanf("%d%d%d%d%d",&n,&m,&s,&t,&k);
	for(int i=1;i<=m;i++)
	scanf("%d%d",&u[i],&v[i]);
	int nfish;
	scanf("%d",&nfish);
	while(nfish--){
		int T;scanf("%d",&T);
		for(int i=0;i<T;i++)
		scanf("%d",&pi[i]);
		for(int i=0;i<12;i++) 
		can[i][pi[i%T]]=1;
	}
	for(int i=0;i<n;i++)can[12][i]=can[0][i];
	for(int i=0;i<n;i++)g[0].a[i][i]=1;
	for(int i=1;i<=12;i++){
		for(int j=1;j<=m;j++){
			if(!can[i-1][u[j]]&&!can[i][v[j]])
			g[i].a[u[j]][v[j]]++;
			if(!can[i-1][v[j]]&&!can[i][u[j]])
			g[i].a[v[j]][u[j]]++;
		}
		g[0]=g[0]*g[i];
	}
	g[0]=g[0]^(k/12);
	for(int i=1;i<=k%12;i++)
	g[0]=g[0]*g[i];
	printf("%lld\n",g[0].a[s][t]);
	return 0;
}