結果

問題 No.563 超高速一人かるた large
ユーザー latte0119latte0119
提出日時 2017-08-26 00:00:56
言語 C++11
(gcc 11.4.0)
結果
AC  
実行時間 1,831 ms / 3,000 ms
コード長 2,327 bytes
コンパイル時間 1,647 ms
コンパイル使用メモリ 169,096 KB
実行使用メモリ 9,672 KB
最終ジャッジ日時 2024-04-23 15:51:46
合計ジャッジ時間 10,830 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge1
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 6 ms
6,912 KB
testcase_01 AC 6 ms
6,912 KB
testcase_02 AC 7 ms
6,912 KB
testcase_03 AC 9 ms
7,108 KB
testcase_04 AC 18 ms
7,112 KB
testcase_05 AC 27 ms
7,240 KB
testcase_06 AC 266 ms
7,880 KB
testcase_07 AC 402 ms
8,264 KB
testcase_08 AC 920 ms
8,776 KB
testcase_09 AC 1,831 ms
9,416 KB
testcase_10 AC 412 ms
7,112 KB
testcase_11 AC 400 ms
7,240 KB
testcase_12 AC 420 ms
7,108 KB
testcase_13 AC 662 ms
7,368 KB
testcase_14 AC 477 ms
7,116 KB
testcase_15 AC 232 ms
6,980 KB
testcase_16 AC 700 ms
7,168 KB
testcase_17 AC 682 ms
9,672 KB
testcase_18 AC 477 ms
9,288 KB
testcase_19 AC 13 ms
9,256 KB
testcase_20 AC 635 ms
8,448 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

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

#define int long long

#define rep(i,n) for(int i=0;i<(n);i++)
#define pb push_back
#define all(v) (v).begin(),(v).end()
#define fi first
#define se second
typedef vector<int>vint;
typedef pair<int,int>pint;
typedef vector<pint>vpint;

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[222222];
int inv[222222];

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

int N;
string S[2222];

typedef unsigned long long ull;
const ull p=1000000009;
vector<ull>H[2222];

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

    cin>>N;
    rep(i,N)cin>>S[i];

    rep(i,N){
        H[i].pb(0);
        rep(j,S[i].size()){
            H[i].pb(H[i].back()*p+S[i][j]);
        }
    }

    int ans[2222]={};
    rep(i,N){
        vint lis;
        rep(j,N){
            if(i==j)continue;
            int lb=0,ub=min(S[i].size(),S[j].size())+1;
            while(ub-lb>1){
                int mid=(ub+lb)/2;
                if(H[i][mid]==H[j][mid])lb=mid;
                else ub=mid;
            }
            lis.pb(ub);
        }
        sort(all(lis));
        reverse(all(lis));

        vint cost,cnt;
        rep(j,lis.size()){
            if(j==0||lis[j-1]!=lis[j]){
                cost.pb(lis[j]);
                cnt.pb(N-1-j);
            }
        }

        for(int k=1;k<=N;k++){
            if(k==N){
                ans[k]=(ans[k]+fact[k-1])%mod;
                continue;
            }
            for(int j=0;j<cost.size();j++){
                if(cnt[j]<N-k)break;
                int tmp=C(cnt[j],N-k);
                if(j+1!=cost.size()&&cnt[j+1]>=N-k)tmp=(tmp-C(cnt[j+1],N-k)+mod)%mod;
                ans[k]=(ans[k]+tmp*cost[j]%mod*fact[k-1])%mod;
            }
        }
    }

    for(int k=1;k<=N;k++){
        ans[k]=(ans[k]+ans[k-1]*(N-k+1))%mod;
        cout<<ans[k]<<endl;
    }
}
0