結果

問題 No.140 みんなで旅行
コンテスト
ユーザー moti
提出日時 2015-04-16 15:56:10
言語 C++11（廃止可能性あり）
(gcc 13.3.0)
結果
AC  
実行時間 12 ms / 5,000 ms
コード長 948 bytes
コンパイル時間 545 ms
コンパイル使用メモリ 65,716 KB
実行使用メモリ 11,904 KB
最終ジャッジ日時 2024-07-04 14:57:21
合計ジャッジ時間 1,665 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge2
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 19
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <iostream>
#include <algorithm>
#include <vector>

using namespace std;

#define REP(i,a,b) for(int i=a;i<(int)b;i++)
#define rep(i,n) REP(i,0,n)

typedef long long ll;

int const MOD = 1e9+7;

int main() {

  ll N; cin >> N;
  vector<vector<ll>> dp(600, vector<ll>(600));

  REP(i, 1, 600) {
    dp[i][1] = 1;
    REP(j, 2, 600) {
      (dp[i][j] = dp[i-1][j-1] + j*dp[i-1][j]) %= MOD;
    }
  }

  vector<vector<ll>> comb(600, vector<ll>(600));
  rep(i, 600) {
    comb[i][0] = 1;
    REP(j, 1, i+1) {
      (comb[i][j] = comb[i-1][j] + comb[i-1][j-1]) %= MOD;
    }
  }

  // (i*(i-1)) ^ j を計算
  vector<vector<ll>> prepow(600, vector<ll>(600));
  rep(i, 600) {
    prepow[i][0] = 1;
    rep(j, 600) {
      (prepow[i][j+1] = prepow[i][j] * i*(i-1)) %= MOD;
    }
  }

  ll ans = 0;
  REP(i, 1, N+1) REP(j, 1, i+1) {
    (ans += comb[N][i] * dp[i][j] % MOD * prepow[j][N-i] % MOD) %= MOD;
  }

  cout << ans << endl;

  return 0;
}
0