P6569 [NOI Online #3 提高组]魔法值--倍增

n<=100,很明显矩阵存储图,设原图为A(n*n)，每次fi都是前i-1转移而来，f[i]=W(1*n)*A^i也就是长度为i的边所影响的节点。

图上可达矩阵普通乘法是 意义:（aij）^k代表ij路径长度为k的条数。

本题为异或W(1*n)^(A运算A.....),这个运算是什么运算呢？和普通矩阵乘法有什么区别？

ai*j代表有aij个wi相异或。根据异或运算的规则，偶数个形同的数异或相当于没有参与运算。因此可以普通矩阵运算再算奇偶性，也可以如果进行普通矩阵运算，注意数据范围，数据值将会很大。或者直接进行矩阵上的异或运算。

将每次的矩阵进行倍增后存储，备用。

对每个询问，进行矩阵快速幂运算，时间复杂度为n*n*n*logn,总时间复杂是为q*n*n*n*logn，只能得40分，还是超时。

因为W为1*n的矩阵，先进行w*A的运算，时间复杂度为n*n,总复杂度变为q*n*n*logn

参考代码：

#include<bits/stdc++.h>
using namespace std;
const int N=110;
typedef long long ll;
ll w[N], A[35][N][N],dp[N][N];
ll n,m,q,u,v;
ll work(ll x) {
	ll ans=0;
	memset(dp,0,sizeof(dp));
	for(int i=1; i<=n; i++)dp[i][i]=1;
	for(int l=0; l<=31&&x; l++) {
		if(x&(1ll<<l)) {
			ll tmp[N][N]= {0};
			for(int i=1; i<=n; i++)
				for(int j=1; j<=n; j++)
					for(int k=1; k<=n; k++) {
						tmp[i][j]^=dp[i][k]*A[l][k][j];
					}
			memcpy(dp,tmp,sizeof(tmp));
			x-=(1<<l);
		}
	}
	for(int i=1; i<=n; i++)if(dp[i][1])ans=ans^w[i];
	return ans;
}
int main() {

	scanf("%lld%lld%lld",&n,&m,&q);
	for(int i=1; i<=n; i++)scanf("%lld",&w[i]);
	for(int i=1; i<=m; i++) {
		scanf("%lld%lld",&u,&v);
		A[0][u][v]=A[0][v][u]=1;
	}
	for(int l=1; l<=31; l++) { //倍增2^l,求出路径长度为2^l的路径条数
		for(int i=1; i<=n; i++)
			for(int j=1; j<=n; j++)
				for(int k=1; k<=n; k++)
					A[l][i][j]^=A[l-1][i][k]*A[l-1][k][j];

	}
	for(int i=1; i<=q; i++) {
		ll x;
		scanf("%lld",&x);
		printf("%lld\n",work(x));
	}
	return 0;
}