第六届蓝桥杯省赛C++B组 垒骰子

垒骰子
赌圣atm晚年迷恋上了垒骰子，就是把骰子一个垒在另一个上边，不能歪歪扭扭，要垒成方柱体。
经过长期观察，atm 发现了稳定骰子的奥秘：有些数字的面贴着会互相排斥！
我们先来规范一下骰子：1 的对面是 4，2 的对面是 5，3 的对面是 6。
假设有 m 组互斥现象，每组中的那两个数字的面紧贴在一起，骰子就不能稳定的垒起来。 
atm想计算一下有多少种不同的可能的垒骰子方式。
两种垒骰子方式相同，当且仅当这两种方式中对应高度的骰子的对应数字的朝向都相同。
由于方案数可能过多，请输出模 10^9 + 7 的结果。
不要小看了 atm 的骰子数量哦～

「输入格式」
第一行两个整数 n m
n表示骰子数目
接下来 m 行，每行两个整数 a b ，表示 a 和 b 数字不能紧贴在一起。

「输出格式」
一行一个数，表示答案模 10^9 + 7 的结果。

「样例输入」
2 1
1 2

「样例输出」
544

「数据范围」
对于 30% 的数据：n <= 5
对于 60% 的数据：n <= 100
对于 100% 的数据：0 < n <= 10^9, m <= 36

资源约定：
峰值内存消耗 < 256M
CPU消耗  < 2000ms

请严格按要求输出，不要画蛇添足地打印类似：“请您输入...” 的多余内容。
所有代码放在同一个源文件中，调试通过后，拷贝提交该源码。
注意: main函数需要返回0
注意: 只使用ANSI C/ANSI C++ 标准，不要调用依赖于编译环境或操作系统的特殊函数。
注意: 所有依赖的函数必须明确地在源文件中 #include <xxx>， 不能通过工程设置而省略常用头文件。

提交时，注意选择所期望的编译器类型。

 

思路：矩阵快速幂，去年刷题刷到这个的时候完全不懂，后来听老师上课讲了觉得自己懂了，不过回去自己写还是写不出来，今年又研究了一下，总算是懂了……网上这道题的题解真的非常的多，觉得每个大佬都讲得很好，我就不重复啦，放一个我做时参考的链接~https://blog.csdn.net/hemk340200600/article/details/61618906

代码：

#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<cmath>
#include<iostream>
#include<algorithm>
#include<string>
#include<vector>
#include<queue>
#include<map>  
#include<set>
using namespace std;

#define MOD 1000000007 
struct matrix
{
    int n,m;
	long long a[7][7];
};
int up[7]={0,4,5,6,1,2,3};

matrix matrix_mul(matrix A,matrix B)
{
	matrix C;
	C.n=A.n;
	C.m=B.m;
	for(int i=1;i<=A.n;i++)
	{
		for(int j=1;j<=B.m;j++)
		{
			C.a[i][j]=0;
			for(int k=1;k<=A.m;k++)
			{
				C.a[i][j]+=(A.a[i][k]*B.a[k][j]%MOD);
				C.a[i][j]%=MOD;
			}
		}
	}
	return C;
}
matrix matrix_pow(matrix A,int n)
{
	matrix res,temp=A;
	res.m=6;
	res.n=6;
	for(int i=1;i<=6;i++)
	{
		for(int j=1;j<=6;j++)
		{
			if(i==j)
				res.a[i][j]=1;
			else
				res.a[i][j]=0;
		}
	}
	while(n)
	{
		if(n&1)
			res=matrix_mul(res,temp);
		temp=matrix_mul(temp,temp);
		n>>=1;
	}
	return res;
}
long long num_pow(long long num,int n)
{
	long long res=1,temp=num;
	while(n)
	{
		if(n&1)
			res=res*temp%MOD;
		temp=temp*temp%MOD;
		n>>=1;
	}
	return res;
}
int main()
{
	int m,i,j,a,b;
	long long n,ans=0;
	matrix coe,base,res;
	cin>>n>>m;
	coe.n=6;
	coe.m=6;
	for(i=1;i<=6;i++)
	{
		for(j=1;j<=6;j++)
			coe.a[i][j]=1;
	}
	for(i=0;i<m;i++) 
	{
		cin>>a>>b;
		coe.a[up[a]][b]=0;
		coe.a[up[b]][a]=0;
	}
	base.n=1;
	base.m=6;
	for(i=1;i<=6;i++)
		base.a[1][i]=1;
	res=matrix_pow(coe,n-1);
	res=matrix_mul(base,res);
	for(i=1;i<=6;i++)
		ans=(ans+res.a[1][i])%MOD;
	ans=ans*num_pow(4,n)%MOD;
	cout<<ans<<endl;
	return 0;
}