AtCoder Product of Arithmetic Progression

/**
链接:https://atcoder.jp/contests/m-solutions2019/tasks/m_solutions2019_e
题意:给出等差数列的首项x,公差d,项数n。求该等差数列的前n项积;
分析:考虑从模数为1e6+3出发,很显然算法复杂度是一个多log的或者更优的写法。
    tag1: d=0时,ans = ksm(x,n);
    tag2: d!=0 时,则:
    ans = x*(x+d)*(x+2*d)*(x+3*d)*...*(x+(n-1)*d)
        = (x/d)*(x/d+1)*(x/d+2)*(x/d+3)*...*(x/d+(n-1))*ksm(d,n)
        由于mod==1e6+3,因此当 x/d + (n-1) >= mod时,ans=0;
        其它情况可以直接预处理出逆元和阶乘直接询问即可;
*/

#include
#define ll long long
using namespace std;

const ll mod=1000003;

ll ksm(ll a,ll b){
    ll ans=1ll;
    while(b){
        if(b&1) ans=ans*a%mod;
        b>>=1;
        a=a*a%mod;
    }
    return ans;
}


const int maxn=2e6+70;
ll fac[maxn],inv[maxn],finv[maxn];//阶乘,逆元,阶乘逆元
void init(){
    inv[0]=inv[1]=fac[0]=fac[1]=finv[0]=finv[1]=1;
    for(int i=2;i=mod) {
        cout<<0<=mod 的情况 
        return ;
    }
    ll ret=ksm(d*1ll,n);
    ret=ret*fac[zx+n-1]%mod;
    ret=ret*finv[zx-1]%mod;
    cout<

 

你可能感兴趣的:(At,Coder)