結果

問題 No.562 超高速一人かるた small
ユーザー ikd
提出日時 2017-08-26 03:44:20
言語 D
(dmd 2.109.1)
結果
AC  
実行時間 958 ms / 3,000 ms
コード長 1,070 bytes
コンパイル時間 692 ms
コンパイル使用メモリ 104,692 KB
実行使用メモリ 12,992 KB
最終ジャッジ日時 2024-06-12 21:47:46
合計ジャッジ時間 12,732 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge2
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 21
権限があれば一括ダウンロードができます

ソースコード

diff # 

void main(){
  import std.stdio, std.string, std.conv, std.algorithm;
  import std.range;

  int n; rd(n);
  auto ss=new string[](n);
  foreach(i; 0..n) rd(ss[i]);

  sort(ss);
  auto mm=new int[][](n, n);
  foreach(i; 0..n)foreach(j; i..n){
    int k=0;
    while(k<ss[i].length && k<ss[j].length
          && ss[i][k]==ss[j][k]) k++;
    mm[i][j]=mm[j][i]=k;
  }

  long mod=10^^9+7;
  auto f=new long[](n+1); f[0]=1;
  foreach(i; 1..f.length) f[i]=i*f[i-1]%mod;

  auto ex=new long[](n+1);
  auto dp=new long[](1<<n); // 0->remaining, 1->used
  foreach(i; 0..(1<<n)){
    int[] v;
    foreach(j; 0..n){
      if(!(i>>j&1)) v~=j;
    }
    foreach(k; 0..v.length){
      auto c=0;
      if(k>0) c=max(c, mm[v[k-1]][v[k]]);
      if(k+1<v.length) c=max(c, mm[v[k]][v[k+1]]);
      (dp[i|(1<<v[k])]+=dp[i]+(c+1)*f[n-v.length])%=mod;
    }
    (ex[n-v.length]+=dp[i])%=mod;
  }

  foreach(k; 1..ex.length) writeln(ex[k]);

}

void rd(T...)(ref T x){
  import std.stdio, std.string, std.conv;
  auto l=readln.split;
  foreach(i, ref e; x){
    e=l[i].to!(typeof(e));
  }
}
0