文件详细信息

下载本文件

本文件的大小为 1120 字节。

#include<cstdio>
#include<cstdlib>
#include<cstring>
typedef long long ll;
const int maxn=100;
const ll P=1000000007;
int popcnt[65536];
void init(){
	for(int i=0;i<65536;i++)popcnt[i]=popcnt[i>>1]+(i&1);
}
int popcount(ll x){
	return popcnt[x>>48]+popcnt[(x>>32)&0xffff]+popcnt[(x>>16)&0xffff]+popcnt[x&0xffff];
}
int n;
ll k,a[maxn];
struct Matrix{
	ll a[maxn][maxn];
	Matrix(){
		memset(a,0,sizeof(a));
	}
}I,start,trans;
Matrix operator*(const Matrix a,const Matrix b){
	Matrix res;
	for(int i=0;i<n;i++)
		for(int j=0;j<n;j++)
			for(int k=0;k<n;k++)
				res.a[i][k]=(res.a[i][k]+a.a[i][j]*b.a[j][k])%P;
	return res;
}
Matrix qpow(Matrix a,ll x){
	Matrix res=I;
	while(x){
		if(x&1)res=res*a;
		a=a*a;
		x>>=1;
	}
	return res;
}
int main(){
	init();
	scanf("%d%lld",&n,&k);
	for(int i=0;i<n;i++)scanf("%lld",a+i);
	for(int i=0;i<n;i++)I.a[i][i]=1;
	for(int i=0;i<n;i++)start.a[0][i]=1;
	for(int i=0;i<n;i++)
		for(int j=0;j<n;j++)
			if(popcount(a[i]^a[j])%3==0)
				trans.a[i][j]=1;
	Matrix ans=start*qpow(trans,k-1);
	ll output=0;
	for(int i=0;i<n;i++)output+=ans.a[0][i];
	printf("%lld\n",output%P);
	return 0;
}