关闭

天气情况

269人阅读 评论(0) 收藏 举报
#include <iostream>
#include <cstdio>
#include <cmath>

using namespace std;

struct M
{
	int n;
	double m[3][3];
	M(int nn=3)
	{  
		n = nn;
		for(int i=0;i<3;i++)
			for(int j=0;j<3;j++)
				m[i][j]=0.0L;
	}
};
M orignal;

M mul(M a,M b)
{
	M c(a.n);
	for(int i=0;i<a.n;i++)
		for(int j=0;j<a.n;j++)
			for(int k=0;k<a.n;k++)
				c.m[i][j] += a.m[i][k]*b.m[k][j];
	return c;
}

M quick(int k)
{
	if(k==1) return orignal;
	if(k==2) return mul(orignal,orignal);
	M temp = quick(k/2);
	if(k%2)
		return mul(temp,mul(temp,orignal));
	else return mul(temp,temp);
	return orignal;
}
void init()
{
	int test;
	cin>>test;

	while (test--)
	{
		for(int i=0;i<3;i++ )
			for(int j=0;j<3;j++)
				cin>>orignal.m[i][j];
		int p;
		cin>>p;
		
		int x,y,z;
		
		while(p--)
		{
			cin>>x>>y>>z;
			M en=quick(z);
			printf("%.3lf\n",en.m[x-1][y-1]);
		}
	}
}
int main()
{ 
	init();
	return 0;
}

0
0

查看评论
* 以上用户言论只代表其个人观点,不代表CSDN网站的观点或立场
    个人资料
    • 访问:100408次
    • 积分:2621
    • 等级:
    • 排名:第13893名
    • 原创:172篇
    • 转载:3篇
    • 译文:1篇
    • 评论:5条
    最新评论