广东工业大学第十四届程序设计竞赛 1010(最后一题,矩阵快速幂)

	题目链接:http://acm.hdu.edu.cn/contests/contest_showproblem.php?pid=1010&cid=843

是真的菜啊,有点难受.

(假)手抖,(真) 菜.jpg. 继续努力了.

题目意思 F(n+1)=F(n)+2*F(n-1)+(n+1)^3. ( ^为乘方运算 );

由题意来看,(显然)这是矩阵快速幂.
1.F(2),不适用公式,前期花了很长时间没推争取.后面直接从F(4)开始加.
2.初始化错了(见图).
3.初始化问题,在自己电脑上没有问题,编译器默认成了0;在队友的电脑上居然就是随机值(现场debug,学会了debug不亏).
4.头脑不清容易混乱,最后要做出来的时候,有的虚了,不够稳定.
5.群消息疯狂助攻下,写了三题
贴一下垃圾代码(硬核初始化(真 · 菜))

#include<bits/stdc++.h>
#define ll long long
using namespace std;
ll T,m,tot,cnt;
ll n;
ll p=123456789;
struct Mat{
	ll M[6][6];
};
struct mat{
	ll m[6];//1*6的矩阵(为啥一开始我写的是  m[1][6],直接一维就可以)
};
Mat Mul1(Mat x){
	Mat c;
	for(int i=0;i<6;i++){
		for(int j=0;j<6;j++){
			c.M[i][j]=0;
		}
	}
	for(int i=0;i<6;i++){
		for(int j=0;j<6;j++){
			for(int k=0;k<6;k++){
				c.M[i][j]=(c.M[i][j]+x.M[i][k]*x.M[k][j])%p;
			}
		}
	}//  矩阵快速里的  x*x  
	return c;
}
mat Mul2(mat x, Mat y){
	mat c;
	
	for(int i=0;i<6;i++)	c.m[i]=0;
	
	for(int i=0;i<6;i++){
			for(int k=0;k<6;k++){
				c.m[i]=(c.m[i]+x.m[k]*y.M[i][k])%p;
			}
	}//矩阵快速幂里的	ans *  x  返回一个mat(1*6)
	return c;
}

mat pow_mod(Mat x,ll y){
		mat ans;
		ans.m[0]=99;
		ans.m[1]=31;
		ans.m[2]=64;
		ans.m[3]=16;
		ans.m[4]=4;
		ans.m[5]=1;//初始化为F(4),手算的
	while(y){
		if(y&1)	ans=Mul2(ans,x);
		x=Mul1(x);
		y>>=1;
	}
	return ans;
}

int main(){
	Mat x;
	x.M[0][0]=1;x.M[0][1]=2;x.M[0][2]=1;x.M[0][3]=3;x.M[0][4]=3;x.M[0][5]=1;
	x.M[1][0]=1;x.M[1][1]=0;x.M[1][2]=0;x.M[1][3]=0;x.M[1][4]=0;x.M[1][5]=0;
	x.M[2][0]=0;x.M[2][1]=0;x.M[2][2]=1;x.M[2][3]=3;x.M[2][4]=3;x.M[2][5]=1;
	x.M[3][0]=0;x.M[3][1]=0;x.M[3][2]=0;x.M[3][3]=1;x.M[3][4]=2;x.M[3][5]=1;
	x.M[4][0]=0;x.M[4][1]=0;x.M[4][2]=0;x.M[4][3]=0;x.M[4][4]=1;x.M[4][5]=1;
	x.M[5][0]=0;x.M[5][1]=0;x.M[5][2]=0;x.M[5][3]=0;x.M[5][4]=0;x.M[5][5]=1;//失智型初始化
	scanf("%lld",&T);
	while(T--){
	scanf("%lld",&n);
	if(n==1) puts("1");
	else if(n==2) puts("2");//前两行不需要,题目给定n>=3;
	else if(n==3)  puts("31");
	else{
			mat y=(pow_mod(x,n-4));
			printf("%lld\n",y.m[0]%p);
 }
	}
	return 0;
}

由于结束就不能提交了,所以这不一定是AC代码(逃
(明天还会更菜么
Dmymax_ly
2019 / 3 / 16