HDU 2256 Problem of Precision （矩阵快速幂）

最新推荐文章于 2019-05-14 15:56:22 发布

hcbbt

最新推荐文章于 2019-05-14 15:56:22 发布

阅读量2.3k

点赞数 1

CC 4.0 BY-SA版权

分类专栏： =====算法相关===== +基础算法 +组合数学文章标签： ACM

本文链接：https://blog.youkuaiyun.com/hcbbt/article/details/38363431

=====算法相关===== 同时被 3 个专栏收录

319 篇文章

订阅专栏

+基础算法

143 篇文章

订阅专栏

+组合数学

19 篇文章

订阅专栏

本文探讨使用矩阵快速幂解决HDU2256ProblemofPrecision中的精度问题，包括题意解析、推导过程和代码实现，重点在于解决高精度计算的难点。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

HDU 2256 Problem of Precision （矩阵快速幂）

ACM

题目地址：HDU 2256 Problem of Precision

题意：
给出一个式子，求值。

分析：
推起来最后那步会比较难想。
具体过程见：

表示共轭只听说过复数的和图的...
这构题痕迹好明显...
跟基友开玩笑说：如果遇到这种题，推到Xn+Yn*sqrt(6)这步时，打表最多只能打到10就爆int了，这是输出正解和Xn，说不定眼神好能发现ans = Xn * 2 - 1呢。= =...

代码：

/*
*  Author:      illuz <iilluzen[at]gmail.com>
*  Blog:        http://blog.youkuaiyun.com/hcbbt
*  File:        2256.cpp
*  Create Date: 2014-08-03 14:27:15
*  Descripton:   
*/

#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
#include <cmath>
using namespace std;
#define repf(i,a,b) for(int i=(a);i<=(b);i++)
typedef long long ll;

const int SIZE = 3;		// max size of the matrix
const int MOD = 1024;

struct Mat{
	int n;
	ll v[SIZE][SIZE];	// value of matrix

	Mat(int _n = SIZE) {
		n = _n;
		memset(v, 0, sizeof(v));
	}

	void init(ll _v) {
		repf (i, 0, n - 1)
			v[i][i] = _v;
	}

	void output() {
		repf (i, 0, n - 1) {
			repf (j, 0, n - 1)
				printf("%lld ", v[i][j]);
			puts("");
		}
		puts("");
	}
} a, b;

Mat operator * (Mat a, Mat b) {
	Mat c(a.n);
	repf (i, 0, a.n - 1) {
		repf (j, 0, a.n - 1) {
			c.v[i][j] = 0;
			repf (k, 0, a.n - 1) {
				c.v[i][j] += (a.v[i][k] * b.v[k][j]) % MOD;
				c.v[i][j] %= MOD;
			}
		}
	}
	return c;
}

Mat operator ^ (Mat a, ll k) {
	Mat c(a.n);
	c.init(1);
	while (k) {
		if (k&1) c = c * a;
		a = a * a;
		k >>= 1;
	}
	return c;
}

void solve(int n) {
	Mat a(2);
	a.v[0][0] = 5;
	a.v[0][1] = 12;
	a.v[1][0] = 2;
	a.v[1][1] = 5;

	a = a ^ (n - 1);
	int xn = 5 * a.v[0][0] + 2 * a.v[0][1];
	printf("%d\n", (xn * 2 - 1) % MOD);
}

int t, n;

int main() {
	scanf("%d", &t);
	while (t--) {
		scanf("%d", &n);
		solve(n);
	}
	return 0;
}