結果
問題 | No.140 みんなで旅行 |
ユーザー | eQe |
提出日時 | 2021-11-11 22:49:43 |
言語 | C++17 (gcc 12.3.0 + boost 1.83.0) |
結果 |
AC
|
実行時間 | 1,506 ms / 5,000 ms |
コード長 | 3,175 bytes |
コンパイル時間 | 5,165 ms |
コンパイル使用メモリ | 273,412 KB |
実行使用メモリ | 6,944 KB |
最終ジャッジ日時 | 2024-05-03 01:13:55 |
合計ジャッジ時間 | 13,007 ms |
ジャッジサーバーID (参考情報) |
judge4 / judge3 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 2 ms
5,248 KB |
testcase_01 | AC | 2 ms
5,248 KB |
testcase_02 | AC | 17 ms
5,376 KB |
testcase_03 | AC | 2 ms
5,376 KB |
testcase_04 | AC | 2 ms
5,376 KB |
testcase_05 | AC | 2 ms
5,376 KB |
testcase_06 | AC | 2 ms
5,376 KB |
testcase_07 | AC | 2 ms
5,376 KB |
testcase_08 | AC | 2 ms
5,376 KB |
testcase_09 | AC | 2 ms
5,376 KB |
testcase_10 | AC | 2 ms
5,376 KB |
testcase_11 | AC | 1,487 ms
5,632 KB |
testcase_12 | AC | 10 ms
5,376 KB |
testcase_13 | AC | 3 ms
5,376 KB |
testcase_14 | AC | 1,506 ms
5,632 KB |
testcase_15 | AC | 1,480 ms
6,940 KB |
testcase_16 | AC | 354 ms
6,940 KB |
testcase_17 | AC | 112 ms
6,940 KB |
testcase_18 | AC | 1,010 ms
6,944 KB |
testcase_19 | AC | 1,183 ms
6,944 KB |
testcase_20 | AC | 79 ms
6,940 KB |
testcase_21 | AC | 2 ms
6,944 KB |
ソースコード
#include <atcoder/all> using namespace atcoder; #include <bits/stdc++.h> #define all(a) (a.begin()),(a.end()) #define rall(a) (a.rbegin()),(a.rend()) using namespace std; using ll=long long; using ld=long double; using pll=pair<ll,ll>; 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;} istream& operator >> (istream& is, modint1000000007& x) { unsigned int t; is >> t; x=t; return is; } istream& operator >> (istream& is, modint998244353& x) { unsigned int t; is >> t; x=t; return is; } ostream& operator << (ostream& os, const modint1000000007& x) { os << x.val(); return os; } ostream& operator << (ostream& os, const modint998244353& x) { os << x.val(); return os; } template<class T1, class T2>istream& operator >> (istream& is, pair<T1,T2>& p) { is >> p.first >> p.second; return is; } template<class T1, class T2>ostream& operator << (ostream& os, const pair<T1,T2>& p) { os << p.first << " " << p.second; return os;} template<class T>istream& operator >> (istream& is, vector<T>& v) { for(T& x:v) is >> x; return is; } template<class T>ostream& operator << (ostream& os, const vector<T>& v) {for(int i=0;i<v.size();i++) {os << v[i] << (i+1 == v.size() ? "":" ");} return os;} template<class... A> void pt() { std::cout << "\n"; } template<class... A> void pt_rest() { std::cout << "\n"; } template<class T, class... A> void pt_rest(const T& first, const A&... rest) { std::cout << " " << first; pt_rest(rest...); } template<class T, class... A> void pt(const T& first, const A&... rest) { std::cout << first; pt_rest(rest...); } static const ll INF=1e18+7; static const double PI=acos(-1); static const double EPS=1e-10; using v1=vector<ll>; using v2=vector<v1>; using v3=vector<v2>; template<class T>struct Combination_mod{ vector<T>fa,rfa,bsum; vector<vector<pair<bool,T>>>mp; ll LIM=5000; Combination_mod(ll n){ fa.resize(n+1); rfa.resize(n+1); fa[0]=rfa[0]=1; for(ll i=1;i<=n;i++)fa[i]=fa[i-1]*i; for(ll i=1;i<=n;i++)rfa[i]=rfa[i-1]/i; bsum.resize(n+1); bsum[0]=1; for(ll i=1;i<=n;i++)bsum[i]=bsum[i-1]+T(-1).pow(i&1)*rfa[i]; if(n<=LIM)mp.resize(n+1,vector<pair<bool,T>>(n+1)); } T C(ll n,ll k){ if(n<k||n<0||k<0) return 0; return fa[n]*rfa[k]*rfa[n-k]; } T H(ll n,ll k){ return C(n-1+k,k); } T S(ll n,ll k){ T res=0; for(ll i=0;i<=k;i++)res+=T(-1).pow((k-i)&1)*C(k,i)*T(i).pow(n); res*=rfa[k]; return res; } T B(ll n,ll k){ T res=0; for(ll i=0;i<=k;i++)res+=T(i).pow(n)*rfa[i]*bsum[k-i]; return res; } T P(ll n,ll k){ if(n==0 and k==0)return 1; if(k==0)return 0; if(mp[n][k].first)return mp[n][k].second; if(n-k>=0)return (mp[n][k]={1,P(n,k-1)+P(n-k,k)}).second; return (mp[n][k]={1,P(n,k-1)}).second; } }; int main(void) { cin.tie(nullptr); ios::sync_with_stdio(false); ll N;cin>>N; using mint=modint1000000007; Combination_mod<mint>comb(N); mint ans=0; for(ll k=1;k<=N;k++)for(ll i=k;i<=N;i++) ans+=comb.C(N,i)*comb.S(i,k)*mint(k*(k-1)).pow(N-i); pt(ans); }