結果

問題 No.503 配列コレクション
コンテスト
ユーザー latte0119
提出日時 2017-04-07 22:42:46
言語 C++11（廃止可能性あり）
(gcc 13.3.0)
結果
AC  
実行時間 65 ms / 2,000 ms
コード長 1,443 bytes
コンパイル時間 1,428 ms
コンパイル使用メモリ 160,288 KB
実行使用メモリ 55,552 KB
最終ジャッジ日時 2024-07-16 02:39:43
合計ジャッジ時間 3,844 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge4
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 25
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

#include<bits/stdc++.h>
using namespace std;

#define int long long

typedef vector<int>vint;
typedef pair<int,int>pint;
typedef vector<pint>vpint;
#define rep(i,n) for(int i=0;i<(n);i++)
#define reps(i,f,n) for(int i=(f);i<(n);i++)
#define all(v) (v).begin(),(v).end()
#define each(it,v) for(__typeof((v).begin()) it=(v).begin();it!=(v).end();it++)
#define pb push_back
#define fi first
#define se second
template<typename A,typename B>inline void chmin(A &a,B b){if(a>b)a=b;}
template<typename A,typename B>inline void chmax(A &a,B b){if(a<b)a=b;}

const int mod=1000000007;


int mpow(int n,int m){
    int ret=1;
    while(m){
        if(m&1)ret=ret*n%mod;
        n=n*n%mod;
        m>>=1;
    }
    return ret;
}
int fact[3333333];
int inv[3333333];

int C(int n,int k){
    return fact[n]*inv[k]%mod*inv[n-k]%mod;
}

int H(int n,int k){
    if(n==0&&k==0)return 1;
    return C(k+n-1,n-1);
}

signed main(){
    fact[0]=1;
    for(int i=1;i<3333333;i++)fact[i]=fact[i-1]*i%mod;
    inv[3333333-1]=mpow(fact[3333333-1],mod-2);
    for(int i=3333333-2;i>=0;i--)inv[i]=inv[i+1]*(i+1)%mod;


    int N,K,D;

    cin>>N>>K>>D;

    int tmp=0;
    int n=N;
    while(n>=K){
        n-=K-1;
        tmp++;
    }

    if(D==1){
        cout<<n<<endl;
        return 0;
    }

    int ans=0;
    int po=1;
    for(int i=0;i<=tmp;i++){
        ans=(ans+H(n-1,tmp-i)*po)%mod;
        po=po*D%mod;
    }

    cout<<ans*n%mod<<endl;
    return 0;
}
0