费马小定理之（阿里巴巴协助征战SARS（困难）

最新推荐文章于 2024-09-18 16:37:53 发布

浦柳人

最新推荐文章于 2024-09-18 16:37:53 发布

阅读量226

点赞数

分类专栏：数论黑洞&思维暴力

第二期水题同时被 2 个专栏收录

42 篇文章 0 订阅

订阅专栏

数论黑洞&思维暴力

36 篇文章 7 订阅

订阅专栏

费马小定理：

如果p是一个质数，而整数a不是p的倍数，则有a^（p-1）≡1（mod p）。
在这里插入图片描述
2与1000000007互质。2 ^ 0 = 2 ^ (p-1) 所以循环节为(p-1-0)=1000000006。

2 ^ a%1000000007 = 2 ^ （a%1000000006) %1000000007。

给出题目：阿里巴巴协助征战SARS（困难）

公式：sum = 2^(n-1) + 4 ^ (n-1);

#include<stdio.h>
#include<iostream>
#include<map>
#include<algorithm>
#include<cstring>
#include<string.h>
#include<string>
#include<math.h>
#include<vector>
#include<map>
using namespace std;
typedef long long ll;
#define MAXN 100005
#define INF 0x3f3f3f3f//将近int类型最大数的一半，而且乘2不会爆int
#define MOD 1000000007

char sh[MAXN];
ll ksm(ll a, ll b)
{
    ll ans=1;
    a%=MOD;
    while(b)
    {
        if(b&1)
            ans = ans*a%MOD;
        a = a*a%MOD;
        b>>=1;
    }
    return ans;
}

int main()
{
    ll n;
    while(scanf("%s", sh) && sh[0] != '0')
    {
        ll cnt=0;
        int len = strlen(sh);
        for(int i=0; i<len; ++i)
        {
            cnt = cnt*10 + (sh[i]-'0');
            if(cnt > MOD-1) cnt%=(MOD-1);	//mod 1000000006
        }
        cnt%=(MOD-1);
        ll sum = ksm(2, cnt-1) + ksm(4, cnt-1);
        cout << sum%MOD << '\n';
    }
	return 0;
}