首页 > 代码库 > hdu4506小明系列故事——师兄帮帮忙 (用二进制,大数高速取余)
hdu4506小明系列故事——师兄帮帮忙 (用二进制,大数高速取余)
Problem Description
小明自从告别了ACM/ICPC之后,就開始潜心研究数学问题了,一则能够为接下来的考研做准备,再者能够借此机会帮助一些同学,尤其是美丽的师妹。这不,班里唯一的女生又拿一道数学题来请教小明,小明当然非常高兴的就接受了。只是等他细致读题以后,发现自己也不会做,这下小明囧了:假设回复说自己不懂,岂不是非常没面子?
所以,他如今私下求你帮忙解决这道题目,题目是这种:
给你n个数字,各自是a1,a2,a3,a4,a5……an,这些数字每过一个单位时间就会改变,如果上一个单位时间的数字为a1’,a2’,a3’……an’,那么这个单位时间的数字a[i] = a[i - 1]’ * K(i == 1的时候a[1] = a[n]’ * K),当中K为给定的系数。
如今的问题就是求第t单位时间的时候这n个数字变成了什么了?因为数字可能会非常大,所以仅仅要你输出数字对10^9 + 7取余以后的结果。
所以,他如今私下求你帮忙解决这道题目,题目是这种:
给你n个数字,各自是a1,a2,a3,a4,a5……an,这些数字每过一个单位时间就会改变,如果上一个单位时间的数字为a1’,a2’,a3’……an’,那么这个单位时间的数字a[i] = a[i - 1]’ * K(i == 1的时候a[1] = a[n]’ * K),当中K为给定的系数。
如今的问题就是求第t单位时间的时候这n个数字变成了什么了?因为数字可能会非常大,所以仅仅要你输出数字对10^9 + 7取余以后的结果。
Input
输入数据第一行是一个正整数T,表示有T组測试数据;
每组数据有两行,第一行包括输入三个整数n, t, k,当中n代表数字个数,t代表第t个单位时间,k代表系数;第二行输入n个数字ai,代表每一个数字開始的时候是多少。
[Technical Specification]
T <= 100
1 <= n <= 10 ^ 4
0 <= t <= 10 ^ 9 当中 t = 0 表示初始状态
1 <= k <= 10 ^ 9
1 <= ai<= 10 ^ 9
每组数据有两行,第一行包括输入三个整数n, t, k,当中n代表数字个数,t代表第t个单位时间,k代表系数;第二行输入n个数字ai,代表每一个数字開始的时候是多少。
[Technical Specification]
T <= 100
1 <= n <= 10 ^ 4
0 <= t <= 10 ^ 9 当中 t = 0 表示初始状态
1 <= k <= 10 ^ 9
1 <= ai<= 10 ^ 9
Output
对于每组数据请输出第t单位时间后这n个数字变成了什么,输出的时候每两个数字之间输出一个空格,行末不要输出多余的空格,详细见例子。
Sample Input
2 3 2 5 1 2 3 3 0 5 1 2 3
Sample Output
50 75 25 1 2 3#include<stdio.h> #define mod 1000000007 int main() { int n,m,k,i,t; __int64 aa[60],ans[100005],sum; scanf("%d",&t); while(t--) { scanf("%d%d%d",&n,&m,&k); for(i=1;i<=n;i++) { scanf("%I64d",&ans[i]);ans[i]%=mod; } int tk=m,ti=0,a[60]; while(tk) { a[++ti]=tk%2; tk/=2; } aa[1]=k%mod; for(i=2; i<=ti;i++) aa[i]=(aa[i-1]*aa[i-1])%mod; sum=1; for(i=1; i<=ti; i++) if(a[i]) sum=(sum*aa[i])%mod; tk=m%n; if(tk)printf("%I64d",(ans[n-tk+1]*sum)%mod),ti=n-tk+1; else printf("%I64d",(ans[1]*sum)%mod),ti=1; if(ti==1) for(i=2;i<=n;i++) printf(" %I64d",(ans[i]*sum)%mod); else { i=ti-1; for(ti++; ti<=n; ti++) printf(" %I64d",(ans[ti]*sum)%mod); for(ti=1;ti<=i; ti++) printf(" %I64d",(ans[ti]*sum)%mod); } printf("\n"); } return 0; }
声明:以上内容来自用户投稿及互联网公开渠道收集整理发布,本网站不拥有所有权,未作人工编辑处理,也不承担相关法律责任,若内容有误或涉及侵权可进行投诉: 投诉/举报 工作人员会在5个工作日内联系你,一经查实,本站将立刻删除涉嫌侵权内容。