HDU 5171

这道题本来很水,以前做过一样的,斐波那契数列,用矩阵快速幂的方法求,本来很水,以前做过很多次,为毛做的时候没想到T_T

#include <iostream>
#include <cstdio>
#include <algorithm>
#include <cmath>
#define LL __int64
using namespace std;

const LL MOD=10000007;
int a[100005];

struct Matrix{
	LL p[2][2];
};
Matrix per,s;

LL tmp[2];

Matrix operator *(Matrix a,Matrix b){
	Matrix c;
	for(int i=0;i<2;i++){
		for(int j=0;j<2;j++){
			c.p[i][j]=0;
			for(int k=0;k<2;k++)
			c.p[i][j]=(c.p[i][j]+a.p[i][k]*b.p[k][j])%MOD;
		}
	}
	return c;
}

Matrix cal_quick(int k){
	Matrix ans=per,p=s;
	while(k){
		if(k&1)
		ans=ans*p;
		k>>=1;
		p=p*p;
	}
	return ans;
}

int main(){
	per.p[0][0]=per.p[1][1]=1;
	per.p[0][1]=per.p[1][0]=0;
	s.p[0][0]=s.p[0][1]=s.p[1][0]=1;
	s.p[1][1]=0;
	int n,k;
	LL first,second,pos;
	LL ans;
	while(scanf("%d%d",&n,&k)!=EOF){
		ans=0; first=second=pos=-1;
		for(int i=0;i<n;i++){
			scanf("%I64d",&a[i]);
			ans=(ans+a[i])%MOD;
			if(a[i]>first){
				first=a[i]; pos=i;
			}
		}
		ans-=first;
		for(int i=0;i<n;i++){
			if(pos!=i&&second<a[i])
			second=a[i];
		}
		tmp[0]=1;tmp[1]=0;
		Matrix one=cal_quick(k+1);
		LL ans_a=(one.p[0][0]*tmp[0]+one.p[0][1]*tmp[1])%MOD;
		ans=((ans+((ans_a-1)*second)%MOD)%MOD+MOD)%MOD;
		one=one*s;
		ans_a=(one.p[0][0]*tmp[0]+one.p[0][1]*tmp[1])%MOD;
		ans=((ans+((ans_a-1)*first)%MOD)%MOD+MOD)%MOD;
		printf("%I64d\n",ans);
	}
	return 0;
}

  

时间: 2024-11-26 01:33:18

HDU 5171的相关文章

hdu 5171 GTY&#39;s birthday gift(数学,矩阵快速幂)

题意: 开始时集合中有n个数. 现在要进行k次操作. 每次操作:从集合中挑最大的两个数a,b进行相加,得到的数添加进集合中. 以此反复k次. 问最后集合中所有数的和是多少. (2≤n≤100000,1≤k≤1000000000) 思路: 写出来发现是要求Fibonaci的前n个数的和. Fibonaci是用矩阵快速幂求的,这个也可以. [Sn,Fn,Fn-1]=[某个矩阵]*[Sn-1,Fn-1,Fn-2] [S2,F2,F1]=[2,1,1] 然后写,,, 这个代码有些繁琐,应该把矩阵操作单独

HDU 5171 GTY&#39;s birthday gift 矩阵快速幂

点击打开链接 GTY's birthday gift Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 65536/65536 K (Java/Others) Total Submission(s): 225    Accepted Submission(s): 78 Problem Description FFZ's birthday is coming. GTY wants to give a gift to ZZF. He as

hdu 5171 GTY&#39;s birthday gift (BestCoder Round #29)

GTY's birthday gift                                                                       Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 65536/65536 K (Java/Others) Total Submission(s): 209    Accepted Submission(s): 71 Problem Description F

HDU 5171 GTY&#39;s birthday gift(矩阵快速幂 )

HDU 5171 GTY's birthday gift ( 矩阵快速幂裸题目 ) #include <cstdio> #include <cstring> #include <algorithm> using namespace std; #define MAX_SIZE 3 #define MOD 10000007 #define clr( a, b ) memset( a, b, sizeof(a) ) typedef long long LL; struct M

hdu 5170 5171

hdu 5170题意: 给定四个整数 a,b,c,d; 要比较a^b 与c^d的大小, 如果数据小的话直接搞就行,现在数据比较大,可以两边同取对数比较: 但是要注意精度问题! #include <iostream> #include <stdio.h> #include <string> #include <string.h> #include <algorithm> #include <cmath> using namespace

hdu 1869 六度分离(最短路floyd)

题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=1869 六度分离 Time Limit: 5000/1000 MS (Java/Others)    Memory Limit: 32768/32768 K (Java/Others)Total Submission(s): 5171    Accepted Submission(s): 2089 Problem Description 1967年,美国著名的社会学家斯坦利·米尔格兰姆提出了一个名为“

HDU 6203 ping ping ping [LCA,贪心,DFS序,BIT(树状数组)]

题目链接:[http://acm.hdu.edu.cn/showproblem.php?pid=6203] 题意 :给出一棵树,如果(a,b)路径上有坏点,那么(a,b)之间不联通,给出一些不联通的点对,然后判断最少有多少个坏点. 题解 :求每个点对的LCA,然后根据LCA的深度排序.从LCA最深的点对开始,如果a或者b点已经有点被标记了,那么continue,否者标记(a,b)LCA的子树每个顶点加1. #include<Bits/stdc++.h> using namespace std;

HDU 5542 The Battle of Chibi dp+树状数组

题目:http://acm.hdu.edu.cn/showproblem.php?pid=5542 题意:给你n个数,求其中上升子序列长度为m的个数 可以考虑用dp[i][j]表示以a[i]结尾的长度为j的上升子序列有多少 裸的dp是o(n2m) 所以需要优化 我们可以发现dp的第3维是找比它小的数,那么就可以用树状数组来找 这样就可以降低复杂度 #include<iostream> #include<cstdio> #include<cstring> #include

hdu 1207 汉诺塔II (DP+递推)

汉诺塔II Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 65536/32768 K (Java/Others)Total Submission(s): 4529    Accepted Submission(s): 2231 Problem Description 经典的汉诺塔问题经常作为一个递归的经典例题存在.可能有人并不知道汉诺塔问题的典故.汉诺塔来源于印度传说的一个故事,上帝创造世界时作了三根金刚石柱子,在一根柱子上从下往