結果

問題 No.140 みんなで旅行
ユーザー sune232002sune232002
提出日時 2015-01-30 00:14:40
言語 C++11
(gcc 11.4.0)
結果
AC  
実行時間 63 ms / 5,000 ms
コード長 4,642 bytes
コンパイル時間 1,315 ms
コンパイル使用メモリ 161,932 KB
実行使用メモリ 50,400 KB
最終ジャッジ日時 2024-06-23 04:19:18
合計ジャッジ時間 3,396 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge2
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 58 ms
50,272 KB
testcase_01 AC 57 ms
50,272 KB
testcase_02 AC 57 ms
50,176 KB
testcase_03 AC 59 ms
50,304 KB
testcase_04 AC 55 ms
50,304 KB
testcase_05 AC 57 ms
50,392 KB
testcase_06 AC 57 ms
50,304 KB
testcase_07 AC 58 ms
50,176 KB
testcase_08 AC 59 ms
50,176 KB
testcase_09 AC 57 ms
50,304 KB
testcase_10 AC 59 ms
50,304 KB
testcase_11 AC 63 ms
50,304 KB
testcase_12 AC 58 ms
50,304 KB
testcase_13 AC 57 ms
50,264 KB
testcase_14 AC 62 ms
50,176 KB
testcase_15 AC 63 ms
50,392 KB
testcase_16 AC 59 ms
50,304 KB
testcase_17 AC 57 ms
50,304 KB
testcase_18 AC 60 ms
50,304 KB
testcase_19 AC 60 ms
50,268 KB
testcase_20 AC 58 ms
50,400 KB
testcase_21 AC 57 ms
50,268 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
using namespace std;
#define REP(i,n) for(int i=0;i<(int)(n);++i)
#define FOR(i,c) for(__typeof((c).begin())i=(c).begin();i!=(c).end();++i)
#define ALL(c) (c).begin(), (c).end()
#define valid(y,x,h,w) (0<=y&&y<h&&0<=x&&x<w)
#define tpl(...) make_tuple(__VA_ARGS__)
const int INF = 0x3f3f3f3f;
const double EPS = 1e-8;
const double PI = acos(-1);
typedef long long ll;
typedef pair<int,int> pii;
template<class T>bool chmax(T &a, const T &b) { if (a<b) { a=b; return 1; } return 0; }
template<class T>bool chmin(T &a, const T &b) { if (b<a) { a=b; return 1; } return 0; }
template<typename Ch,typename Tr,typename C,typename=decltype(begin(C()))>basic_ostream<Ch,Tr>& operator<<(basic_ostream<Ch,Tr>&os,
const C& c){os<<'[';for(auto i=begin(c);i!=end(c);++i)os<<(i==begin(c)?"":" ")<<*i;return os<<']';}
template<class S,class T>ostream&operator<<(ostream &o,const pair<S,T>&t){return o<<'('<<t.first<<','<<t.second<<')';}
template<int N,class Tp>void output(ostream&,const Tp&){}
template<int N,class Tp,class,class ...Ts>void output(ostream &o,const Tp&t){if(N)o<<',';o<<get<N>(t);output<N+1,Tp,Ts...>(o,t);}
template<class ...Ts>ostream&operator<<(ostream&o,const tuple<Ts...>&t){o<<'(';output<0,tuple<Ts...>,Ts...>(o,t);return o<<')';}
template<class T>void output(T t,char z=10){if(t<0)t=-t,putchar(45);int c[20];
int k=0;while(t)c[k++]=t%10,t/=10;for(k||(c[k++]=0);k;)putchar(c[--k]^48);putchar(z);}
template<class T>void outputs(T t){output(t);}
template<class S,class ...T>void outputs(S a,T...t){output(a,32);outputs(t...);}
template<class T>void output(T *a,int n){REP(i,n)cout<<a[i]<<(i!=n-1?',':'\n');}
template<class T>void output(T *a,int n,int m){REP(i,n)output(a[i],m);}
template<class T>bool input(T &t){int n=1,c;for(t=0;!isdigit(c=getchar())&&~c&&c-45;);
if(!~c)return 0;for(c-45&&(n=0,t=c^48);isdigit(c=getchar());)t=10*t+c-48;t=n?-t:t;return 1;}
template<class S,class ...T>bool input(S&a,T&...t){input(a);return input(t...);}
template<class T>bool inputs(T *a, int n) { REP(i,n) if(!input(a[i])) return 0; return 1;}
template<int MOD>
struct ModInt {
  int x;
  ModInt() : x(0) {}
  ModInt(ll s) { if ((x = s % MOD) < 0) x += MOD; }
  ModInt operator+=(ModInt rhs) { if ((x+=rhs.x) >= MOD) x -= MOD; return *this; }
  ModInt operator-=(ModInt rhs) { if ((x-=rhs.x) < 0) x += MOD; return *this; }
  ModInt operator*=(ModInt rhs) { x = (ll)x*rhs.x % MOD; return *this; }
  ModInt operator+(ModInt rhs) const { return ModInt(*this) += rhs; }
  ModInt operator-(ModInt rhs) const { return ModInt(*this) -= rhs; }
  ModInt operator*(ModInt rhs) const { return ModInt(*this) *= rhs; }
  ModInt operator/=(ModInt rhs) {
    static const ll inv2 = ModInt(2).inv().x; // 2で割るのは特別に保存してみる
    ll i = (rhs.x == 2 ? inv2 : rhs.inv().x);
    x = x*i%MOD; return *this;
  }
  ModInt operator/(ModInt rhs) const { return ModInt(*this) /= rhs; }
  ModInt inv() { return pow(MOD-2); }
  ModInt pow(ll n) const {
    ModInt r = 1, t = x;
    for (;n;n>>=1,t*=t) if (n&1) r *= t;
    return r;
  }
  friend ostream& operator<<(ostream &os, const ModInt rhs){ return os << rhs.x; }
};
typedef ModInt<int(1e9+7)> mint;

// 順番を考慮した分割数
const int MAX_N = 2000;
mint part[MAX_N][MAX_N];
// void calc_part(int n, int k) {
//   memset(part, 0, sizeof(part));
//   part[0][0] = 1;
//   for (int i=0; i<=n; ++i)
//     for (int j=1; j<=k; ++j)
//       part[i][j] = part[i][j-1] + (i>0?part[i-1][j]:0);
// }
void calc_part(int n, int k) {
  memset(part, 0, sizeof(part));
  part[0][0] = 1;
  for (int i=0; i<=n; ++i)
    for (int j=1; j<=k; ++j)
      part[i][j] = part[i][j-1] + (i-j>=0?part[i-j][j]:0);
}

const int NUM = 2000;
mint combi[NUM][NUM];
void calcCombi() {
  combi[0][0] = 1;
  REP(i, NUM-1) {
    combi[i][0] = combi[i][i] = 1;
    REP(j, i)
      combi[i+1][j+1] = combi[i][j] + combi[i][j+1];
  }
}
mint dp[2000][2000];
int main() {
  memset(dp,0,sizeof(dp));
  dp[0][0] = 1;
  REP(i,2000) {
    REP(j,2000) {
      if (i>0&&j>0) dp[i][j] += dp[i-1][j-1];
      if (i>0) dp[i][j] += dp[i-1][j]*j;
    }
  }
  calc_part(1500,1500);
  calcCombi();
  int n;
  while(cin>>n) {
    mint ans = 0;
    for (int g=1; g<=n; ++g) {
      for (int a=g; a<=n; ++a) {
        // mint tmp = part[a][g] - part[a][g-1];
        mint tmp = dp[a][g];
        // cout << tmp << endl;
        tmp *= combi[n][a];
        // cout << a << " "<< g << " " << part[a][g] << " " << part[a][g-1] << " " << tmp << " " << endl;
        tmp *= mint(g*(g-1)).pow(n-a);
        ans += tmp;
        // cout << g << " " << a << " " << tmp << endl;
      }

    }
    cout << ans << endl;
  }

}
0