ACM-ICPC 2018 焦作赛区网络预赛 G. Give Candies 打表+指数循环节 or欧拉降幂（费马小定理） 一题多解

博客目录

原题

传送门

•  26.61%
•  1000ms
•  65536K

There are NN children in kindergarten. Miss Li bought them NN candies. To make the process more interesting, Miss Li comes up with the rule: All the children line up according to their student number (1...N)(1...N), and each time a child is invited, Miss Li randomly gives him some candies (at least one). The process goes on until there is no candy. Miss Li wants to know how many possible different distribution results are there.

Input

The first line contains an integer TT, the number of test case.

The next TT lines, each contains an integer NN.

1 \le T \le 1001≤T≤100

1 \le N \le 10^{100000}1≤N≤10100000

Output

For each test case output the number of possible results (mod 1000000007).

样例输入复制

1
4

样例输出复制

8

题目来源

ACM-ICPC 2018 焦作赛区网络预赛

思路

通过打表我们发现，结果是pow(2,n-1)

然后本题可以用两种解法，第一个是欧拉降幂，第二个是指数循环节。我用的第二种。

欧拉降幂只能用于对质数取模的情况（不然可能要用到中国剩余定理反正各种麻烦而且复杂度就飙升了），质数循环节则无所谓质数合数

指数循环节：

寻找循环节长度的代码：

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
ll const mod=1e9+7;
int main(){// 1 2 4 8 16
    int i;
    f[1]=1;
    ll w=1;
    for(i=1;i<=mod;i++){
        w*=2;
        w%=mod;
        if(w==1){
            cout<<i<<endl;
            break;
        }
    }
}

输出是500000003

也就是每500000003循环一次，所以有：

2^n%mod=2^（n%500000003）

然后我们用大数运算计算下式：

 

fpow(2,n%500000004)

其中n%500000004是大数取模运算，fpow是快速幂

欧拉降幂(费马小定理)：

由欧拉公式有：phi(p)=p-1,(p为质数)    pow(a,phi(p))%p==1

所以有：

pow(a,b)%p=pow(a,b%(p-1)+p-1)%p

AC代码（指数循环节）

#include<cstring>
#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
char s[110000];
ll divMod(char* ch,ll num)
{
    ll s = 0;
    for(int i=0; ch[i]!='\0'; i++)
        s = (s*10+ch[i]-'0')%num;
        s=((s-1+num)%num+num)%num;
    return s;
}
ll PowerMod(ll a, ll b, ll c)

{

    ll ans = 1;
    a = a % c;
    while(b>0)
    {
        if(b % 2 == 1)
            ans = (ans * a) % c;
        b = b/2;
        a = (a * a) % c;
    }
    return ans;

}
int main()
{
    int t;
    ll n;
    cin>>t;
    while(t--)
    {
        cin>>s;
        ll d=divMod(s,500000003);
        cout<<PowerMod(2,d,1000000007)<<endl;
    }
    return 0;
}

欧拉降幂（费马小定理）

#include<cstring>
#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
char s[110000];
ll divMod(char* ch,ll num)
{
    ll s = 0;
    for(int i=0; ch[i]!='\0'; i++)
        s = (s*10+ch[i]-'0')%num;
        s=((s-1+num)%num+num)%num;
    return s;
}
ll PowerMod(ll a, ll b, ll c)

{

    ll ans = 1;
    a = a % c;
    while(b>0)
    {
        if(b % 2 == 1)
            ans = (ans * a) % c;
        b = b/2;
        a = (a * a) % c;
    }
    return ans;

}
int main()
{
    int t;
    ll n;
    cin>>t;
    while(t--)
    {
        cin>>s;
        ll d=divMod(s,1000000007-1);
        cout<<PowerMod(2,d,1000000007)<<endl;
    }
    return 0;
}