关闭

hdu4565 So Easy!(矩阵快速幂)

标签: hdu
1360人阅读 评论(1) 收藏 举报
分类:
题意:
告诉你a,b,n,m,求

0<a,m<2^15 , (a-1)^2<b<a^2 , 0<b,n<2^31
解题思路:
    我们首先对(a+sqrt(b))^n进行处理,a+sqrt(b))^n的展开式我们可以知道
(a+sqrt(b))^n = Xn + Yn * sqrt(b);
那么
( a + sqrt( b ) )^( n + 1 ) = ( a + sqrt( b ) ) * ( Xn + Yn * sqrt(b) ) = ( a * Xn + b * Yn ) + ( a * Yn + Xn ) * sqrt( b );
Xn+1 = ( a * Xn + b * Yn )
Yn+1 = ( a * Yn + Xn )
将其转化为矩阵,则

递推下去可以得到

所以
Xn = a * X0 + b * Y0;
Yn = X0 + a * Y0;

由于 (a-1)^2<b<a^2 ,
那么 0 < ( a - sqrt( b ) )^ n < 1
 又  ( a + sqrt( b ) )^ n + ( a - sqrt( b ) )^ n = 2 * Xn
那么 ( a + sqrt( b ) )^ n 向上取整的值就是 2 * Xn
最终
    Sn = ( 2 * Xn )% m;

参考代码:
#include <iostream>
#include <math.h>
#include <stdio.h>
using namespace std;
typedef long long ll;
struct Matrix{
	ll mat[2][2];
};
Matrix mul(Matrix a,Matrix b,ll mod){
	Matrix ans;
	for (int i=0;i<2;i++){
		for (int j=0;j<2;j++){
			ans.mat[i][j]=0;
			for (int k=0;k<2;k++){
				ans.mat[i][j]=(ans.mat[i][j]+a.mat[i][k]*b.mat[k][j]);
				ans.mat[i][j]%=mod;
			}
		}
	}
	return ans;
}
Matrix Init(){
	Matrix ans;
	for (int i=0;i<2;i++){
		for (int j=0;j<2;j++){
			if (i==j)
				ans.mat[i][j]=1;
			else
				ans.mat[i][j]=0;
		}
	}
	return ans;
}
Matrix exp(Matrix a,ll k,ll m){
	Matrix ans=Init();
	while (k){
		if (k&1)
			ans=mul(ans,a,m);
		a=mul(a,a,m);
		k>>=1;
	}
	return ans;
}
int main(){
	Matrix M;
	ll a,b,n;
	ll m;
	while (~scanf("%lld%lld%lld%lld",&a,&b,&n,&m)){
		M.mat[0][0]=M.mat[1][1]=a;
		M.mat[0][1]=b;
		M.mat[1][0]=1;
		Matrix ans=exp(M,n,m);
		/*
		for (int i=0;i<2;i++){
			for (int j=0;j<2;j++)
				cout<<ans.mat[i][j]<<" ";
			cout<<endl;
		}
		*/
		ll x=ans.mat[0][0];

		cout<<(2*x)%m<<endl;
	}
	return 0;
}

1
0
查看评论
发表评论
* 以上用户言论只代表其个人观点,不代表CSDN网站的观点或立场

HDU4565 So Easy! (矩阵)

HDU4565 So easy (矩阵+推公式)
  • u013790563
  • u013790563
  • 2015-05-19 11:22
  • 773

hdu 4565 So Easy! /2013 ACM-ICPC 长沙赛区全国邀请赛A题 矩阵乘法

#include #include #include #include #include #include using namespace std; #define LL __int64 ...
  • a601025382s
  • a601025382s
  • 2013-08-18 14:05
  • 2086

hdu4565---So Easy!(矩阵)

Problem Description   A sequence Sn is defined as:Where a, b, n, m are positive integers.┌x┐is the ...
  • Guard_Mine
  • Guard_Mine
  • 2015-03-14 11:32
  • 657

HDU 4565 So Easy! (矩阵快速幂)

题目地址:http://acm.hdu.edu.cn/showproblem.php?pid=4565 题意:给你一个a,b,n,m。其中(a-1)^2 Sn=⌈(a+b√)n⌉%m,(a−1)2...
  • ilovexiaohao
  • ilovexiaohao
  • 2013-08-14 19:42
  • 1511

hdu4565之矩阵快速幂

矩阵快速幂
  • xingyeyongheng
  • xingyeyongheng
  • 2013-08-05 16:13
  • 899

HDU 4565 So Easy!

题解: 因为b有一个范围,即a-1 Cn*((a+sqrt(b))+(a-sqrt(b))),可得递推公式Cn+1=2a*Cn-(a^2-b)*Cn-1,而Pn+1=aPn+bPn-1,这种形式可...
  • zamac1993
  • zamac1993
  • 2013-09-07 01:52
  • 375

HDU4565/2013年长沙赛区So easy

Description A sequence Sn is defined as: Where a,b,n,m are positiveintegers.┌x┐is the ceil ofx. ...
  • lphy2352286B
  • lphy2352286B
  • 2013-06-03 21:16
  • 1657

ZOJ 3690 & HDU 3658 (矩阵快速幂+公式递推)

ZOJ 3690题意:有n个人和m个数和一个k,现在每个人可以选择一个数,如果相邻的两个人选择相同的数,那么这个数要大于k 求选择方案数。思路:打表推了很久的公式都没推出来什么可行解,好不容易有了想...
  • qq_15714857
  • qq_15714857
  • 2015-08-16 23:36
  • 2443

hdu4565 so easy 矩阵

的确和题目描述一样,见过这个类型就so easy了 然而无脑的我还是调了一个晚上,就因为写乘法的时候忘记return了 然后就写了两个版本的,其实并没有多大不同 #include #include #...
  • BPM136
  • BPM136
  • 2016-03-09 21:17
  • 167

HDU 4565 So Easy!

题目分析Sn=⌈(a+b√)n⌉%m S_n = \lceil(a+\sqrt{b})^n\rceil \% m 因为(a−1)2<b<a2(a-1)^2 < b < a^2 ,所以0<(a−...
  • chen_ze_hua
  • chen_ze_hua
  • 2016-07-30 15:37
  • 123
    个人资料
    • 访问:144984次
    • 积分:3016
    • 等级:
    • 排名:第13397名
    • 原创:157篇
    • 转载:0篇
    • 译文:0篇
    • 评论:29条
    最新评论