ACM-ICPC 2018 焦作赛区网络预赛 Give Candies 题解

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

n个糖果分给n个小朋友
从1到n个小朋友依次给,每次随机给个数,至少一个,知道没有糖果为止。
问糖果的分布情况方案数。
输出方案数mod 109+710^9+7109+7

考虑只有前i个小朋友得到糖的情况,于是等价于将n个糖果分为i堆,插板法易得方案数是(n−1i−1)\binom{n-1}{i-1}(i1n1)
总方案数∑i=1n(n−1i−1)=2n−1\sum_{i=1}^{n}\binom{n-1}{i-1}=2^{n-1}i=1n(i1n1)=2n1
2n−1mod  10000000072^{n-1} \mod 10000000072n1mod1000000007
anmod  pa^n \mod panmodp
p是质数,只是n很大
an≡anmod  ϕ(p)(modp)a^n \equiv a^{n \mod \phi(p)} \pmod{p}ananmodϕ(p)(modp)
依据是费马-欧拉定理

更一般的情况简记

事实上,更为一般的是:
gcd(a,c)=1⇒ab≡abmod  ϕ(c)(modc)gcd(a,c)=1 \Rightarrow a^b \equiv a^{b \mod \phi(c)} \pmod{c}gcd(a,c)=1ababmodϕ(c)(modc)
如果a,c不互素呢?
b>ϕ(n)⇒ab≡abmod  ϕ(c)  +  ϕ(c)(modc)b \gt \phi(n) \Rightarrow a^b \equiv a^{b \mod \phi(c)\; + \;\phi(c)} \pmod{c}b>ϕ(n)ababmodϕ(c)+ϕ(c)(modc)
如果b≤ϕ(n)b \leq \phi(n)bϕ(n);那就不用换了。

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
ll n,m;
const ll mod = 1e+9+7; // is prime
const ll phi_mod = mod-1;

// pre: mod != 0, <a,n>!=<0,0> n>=0
ll mlt(ll a, ll n, ll mod) {
    if (n == 0)
        return 1;
    ll t = 1;
    a %= mod;
    while (n > 1) {
        if (n&1)
            t = (t*a)%mod;
        a = (a*a)%mod;
        n >>= 1;
    }
    return (t*a)%mod;
}
string s;
int main() {
    ios::sync_with_stdio(false);
    cin.tie(0);
    int t;
    cin >> t;
    while (t--) {
        cin>>s;
        n = 0;
        for (auto x : s)
            n = ((n*10)+x-'0')%phi_mod;
        n = (n-1+phi_mod)%phi_mod;
        cout<<mlt(2ll,n,mod)<<endl;
    }
    return 0;
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值