HDU 小明系列故事——师兄帮帮忙 高速幂

小明系列故事——师兄帮帮忙

Time Limit: 3000/1000 MS (Java/Others)    Memory Limit: 65535/32768 K (Java/Others)

Total Submission(s): 4850    Accepted Submission(s): 1275

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取余以后的结果。

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

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

高速幂问题，发现最后结果是  a[i]*k^t 然后依据 t 和 n的关系把序列又一次排一下就OK啦

#include <iostream>
#include <stdio.h>
#include <cmath>
#include <cstring>
#include <string>
#include <algorithm>
#define N 10009
#define ll __int64
#define mod 1000000007

using namespace std;

ll n,t,k,T;
ll a[N];

ll pow_mod(ll m,ll n)
{
    ll b=1;
    while(n)
    {
        if(n&1)
        b=(b*m)%mod;
        n=n>>1;
        m=m*m%mod;
    }
    return b;
}

ll b[N];

int main()
{
   while(~scanf("%I64d",&T))
   {
        while(T--)
        {
            scanf("%I64d %I64d %I64d",&n,&t,&k);
            for(ll i=0;i<n;i++)
            scanf("%I64d",&a[i]);

            ll mm=pow_mod(k,t);
            for(ll i=0;i<n;i++)
            {
                a[i]=((mm*a[i])%mod);
            }

            ll f=t%n;

            for(ll i=0;i<n;i++)
            {
                if(i+f>=n)
                b[(i+f)%n]=a[i];
                else
                b[i+f]=a[i];
            }

            for(ll i=0;i<n-1;i++)
            printf("%I64d ",b[i]);
            printf("%I64d\n",b[n-1]);

        }
   }

    return 0;
}