费马小定理

upupup

Time Limit : 3000/1000ms (Java/Other)   Memory Limit : 65535/32768K (Java/Other)

Total Submission(s) : 239   Accepted Submission(s) : 14

Font: Times New Roman | Verdana | Georgia

Font Size: ← →

Problem Description

最近侯ry感觉自己在数学方面的造诣不忍直视；他发现他的学习速率呈一个指数函数递增，疯狂的陷入学习的泥潭，无法自拔；他的队友发现了他的学习速率y=e^（b*lna+lnc);
e是科学界非常重要而常见的常数，e=2.718281828……。
侯ry由于数学很差不会算学习数率y，现求助于学弟，感激不尽；

Input

多组数据，每组数据输入三个整数a，b，c（a，c<=10^12,b<=10^100000)

Output

一个整数y，对10^9+7取模

Sample Input

2 3 10

Sample Output

80

　　　　费马小定理(Fermat Theory)是数论中的一个重要定理，其内容为： 假如p是质数，且gcd(a,p)=1，那么 a^(p-1)≡1（mod p）。即：假如a是整数，p是质数，且a,p互质(即两者只有一个公约数1)，那么a的(p-1)次方除以p的余数恒等于1。

　　　周末打的场比赛，里面的一道题目，用费马小定理可以分分钟秒杀hhh，这道题可以直接转化成(a^b)*c的问题，但由于b是个大数，如果看懂了费
马小定理，这时候应该就知道怎么解决了，直接用一个高精度取余单精度就可以把这个棘手的问题解决了，将b直接转化成了一个比模数10^9+7还小的数，之
后就一个快速幂就可以得到结果了，是不是很腻害的样子>.<

#include<iostream>
#define LL long long
#define N 1000000007
using namespace std;

struct num
{
    static const int L=100005;
    int a[L];
    int len;
    void ini(string s)//字符串输入
    {
        fill(a,a+L,0);
        for(int i=s.length()-1;i>=0;i--)
            a[s.length()-i-1]=s[i]-'0';
        len=s.length();
    }
};
LL mod(num a)//高精度取余单精度 输出余数
{
    LL temp=0;
    for(int i=a.len-1; i>=0; i--)
        temp=(temp*10+a.a[i])%(N-1);
    return temp;
}
LL quickpow(LL m,LL n)
{
    LL b = 1;
    while (n > 0)
    {
          if (n & 1)
             b = (b*m)%N;
          n = n >> 1 ;
          m = (m*m)%N;
    }
    return b;
}
int main()
{
    LL a,c;
    num str;
    string ss;
    while(cin>>a>>ss>>c)
    {
        str.ini(ss);
        LL temp=mod(str);
        c=c%N;
        a%=N;
        LL ans=(quickpow(a,temp)*c)%N;
        cout<<ans<<endl;
    }
    return 0;
}