結果
問題 | No.122 傾向と対策:門松列(その3) |
ユーザー | 沙耶花 |
提出日時 | 2021-10-26 22:35:48 |
言語 | C++17 (gcc 12.3.0 + boost 1.83.0) |
結果 |
AC
|
実行時間 | 100 ms / 5,000 ms |
コード長 | 1,629 bytes |
コンパイル時間 | 5,766 ms |
コンパイル使用メモリ | 270,048 KB |
実行使用メモリ | 22,068 KB |
最終ジャッジ日時 | 2024-10-06 04:38:07 |
合計ジャッジ時間 | 6,611 ms |
ジャッジサーバーID (参考情報) |
judge1 / judge2 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 96 ms
21,936 KB |
testcase_01 | AC | 96 ms
22,068 KB |
testcase_02 | AC | 99 ms
21,936 KB |
testcase_03 | AC | 97 ms
21,940 KB |
testcase_04 | AC | 96 ms
21,936 KB |
testcase_05 | AC | 98 ms
21,936 KB |
testcase_06 | AC | 97 ms
21,932 KB |
testcase_07 | AC | 100 ms
21,932 KB |
ソースコード
#include <stdio.h> #include <bits/stdc++.h> #include <atcoder/all> using namespace atcoder; using mint = modint1000000007; using namespace std; #define rep(i,n) for (int i = 0; i < (n); ++i) #define Inf 1000000000 vector<mint> ika(vector<int> mini,vector<int> maxi){ int n = mini.size(); vector dp(20001,vector(1<<n,vector<mint>(2,0))); dp[0][0][0] = 1; for(int i=0;i<20000;i++){ rep(j,1<<n){ rep(k,2){ dp[i+1][j][0] += dp[i][j][k]; rep(l,n){ if((j>>l)&1)continue; if(mini[l]<=i+1&&i+1<=maxi[l]){ dp[i+1][j|(1<<l)][1] += dp[i][j][k]; } } } } } vector<mint> ret(20001,0); rep(i,dp.size()){ ret[i] += dp[i][(1<<n)-1][1]; } return ret; } vector<mint> more(vector<int> mini,vector<int> maxi){ int n = mini.size(); vector dp(20001,vector(1<<n,vector<mint>(2,0))); dp.back()[0][0] = 1; for(int i=20000;i>=1;i--){ rep(j,1<<n){ rep(k,2){ dp[i-1][j][0] += dp[i][j][k]; rep(l,n){ if((j>>l)&1)continue; if(mini[l]<=i&&i<=maxi[l]){ dp[i-1][j|(1<<l)][1] += dp[i][j][k]; } } } } } vector<mint> ret(20001,0); rep(i,dp.size()){ rep(j,2){ ret[i] += dp[i][(1<<n)-1][j]; } } return ret; } int main(){ vector<int> m(7),M(7); rep(i,7)cin>>m[i]>>M[i]; vector<int> x = {m[0],m[2],m[4],m[6]}; vector<int> y = {M[0],M[2],M[4],M[6]}; vector<int> z = {m[1],m[3],m[5]}; vector<int> w = {M[1],M[3],M[5]}; auto a = ika(x,y); auto b = more(z,w); auto c = more(x,y); auto d = ika(z,w); mint ans = 0; rep(i,a.size()){ ans += a[i] * b[i]; ans += c[i] * d[i]; } cout<<ans.val()<<endl; return 0; }