結果

問題 No.562 超高速一人かるた small
ユーザー ei1333333ei1333333
提出日時 2017-08-25 23:51:56
言語 C++11
(gcc 11.4.0)
結果
AC  
実行時間 1,094 ms / 3,000 ms
コード長 1,830 bytes
コンパイル時間 1,511 ms
コンパイル使用メモリ 162,428 KB
実行使用メモリ 5,248 KB
最終ジャッジ日時 2024-10-15 16:47:28
合計ジャッジ時間 15,247 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge3
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
5,248 KB
testcase_01 AC 2 ms
5,248 KB
testcase_02 AC 2 ms
5,248 KB
testcase_03 AC 2 ms
5,248 KB
testcase_04 AC 2 ms
5,248 KB
testcase_05 AC 3 ms
5,248 KB
testcase_06 AC 2 ms
5,248 KB
testcase_07 AC 7 ms
5,248 KB
testcase_08 AC 111 ms
5,248 KB
testcase_09 AC 1,090 ms
5,248 KB
testcase_10 AC 236 ms
5,248 KB
testcase_11 AC 509 ms
5,248 KB
testcase_12 AC 112 ms
5,248 KB
testcase_13 AC 5 ms
5,248 KB
testcase_14 AC 1,088 ms
5,248 KB
testcase_15 AC 1,088 ms
5,248 KB
testcase_16 AC 1,094 ms
5,248 KB
testcase_17 AC 1,087 ms
5,248 KB
testcase_18 AC 1,089 ms
5,248 KB
testcase_19 AC 1,088 ms
5,248 KB
testcase_20 AC 1,090 ms
5,248 KB
testcase_21 AC 1,087 ms
5,248 KB
testcase_22 AC 1,090 ms
5,248 KB
testcase_23 AC 1,091 ms
5,248 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include<bits/stdc++.h>

using namespace std;

typedef long long int64;
const int mod = 1e9 + 7;

inline int64 modPow(int64 x, int64 n)
{
  if(n == 0) return (1);
  int64 ret = modPow(x, n / 2);
  (ret *= ret) %= mod;
  if(n & 1) (ret *= x) %= mod;
  return (ret);
}

inline int64 modInv(int64 a)
{
  return (modPow(a, mod - 2));
}

int N;
string S[2000];
int dat[2000][2000];
int64 fact[2001], inv[2001];

int main()
{
  fact[0] = 1;
  for(int i = 1; i < 2001; i++) {
    (fact[i] = fact[i - 1] * i) %= mod;
  }
  for(int i = 0; i < 2001; i++) inv[i] = modInv(fact[i]);


  cin >> N;
  for(int i = 0; i < N; i++) {
    cin >> S[i];
    S[i] += "*";
  }

  for(int i = 0; i < N; i++) {
    dat[i][i] = 1;
    for(int j = 0; j < N; j++) {
      if(i == j) continue;
      int diff = 0;
      while(S[i][diff] == S[j][diff]) ++diff;
      dat[i][j] = diff + 1;
    }
  }


  int dp[1 << 20] = {};
  dp[0] = 0;

  vector< int > ans(N + 1, 0);

  for(int i = 0; i < (1 << N); i++) {

    vector< int > yuki(N, 1);
    for(int j = 0; j < N; j++) {
      if(!((i >> j) & 1)) {
        for(int k = 0; k < N; k++) {
          if(!((i >> k) & 1)) yuki[j] = max(yuki[j], dat[j][k]);
        }
      }
    }

    int sz = __builtin_popcount(i);

    for(int j = 0; j < N; j++) {
      if((i >> j) & 1) continue;
      // (dp[i | (1 << j)] += (dp[i] + 1LL * yuki[j] * fact[sz] % mod) % mod) %= mod;
      // (アが及ぼすアへの道のりは 残りbit数!)
      for(int k = sz + 1; k <= N; k++) {
        int rest = N - (sz + 1); // 現時点での残りbit
        (ans[k] += yuki[j] * fact[sz] % mod * fact[rest] % mod * inv[N - k] % mod) %= mod;
      }

    }

  }

/*
  for(int i = 0; i < (1 << N); i++) {
    (ans[__builtin_popcount(i)] += dp[i]) %= mod;
  }
  */
  for(int i = 1; i <= N; i++) {
    cout << ans[i] << endl;
  }

}
0