結果
問題 | No.93 ペガサス |
ユーザー | omochana1 |
提出日時 | 2014-12-15 01:11:41 |
言語 | C++11 (gcc 13.3.0) |
結果 |
AC
|
実行時間 | 37 ms / 5,000 ms |
コード長 | 2,104 bytes |
コンパイル時間 | 924 ms |
コンパイル使用メモリ | 87,824 KB |
実行使用メモリ | 22,656 KB |
最終ジャッジ日時 | 2024-06-11 21:10:42 |
合計ジャッジ時間 | 1,422 ms |
ジャッジサーバーID (参考情報) |
judge2 / judge4 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 4 |
other | AC * 16 |
ソースコード
#include <iostream> #include <cmath> #include <cstdio> #include <algorithm> #include <string> #include <vector> #include <queue> #include <string.h> #include <map> #include <fstream> #include <functional> #include <bitset> #include <iomanip> #include <stack> #include <set> #include <climits> #define MAX_N 200010 #define PI 3.141592653589 #define ESP 1e-20 #define BS 10 #define MOD 1000000007 #define ZERO 10001 #define YJSNPI 810 #define INF (1LL << 50) #define ADD(a, b) a = (a + (ll)b) % MOD #define MUL(a, b) a = (a * (ll)b) % MOD #define MAX(a, b) a = max(a, b) #define MIN(a, b) a = min(a, b) using namespace std; typedef pair<int, int> pi; typedef long long ll; int N; ll dp[1010][1010][2][2]; int main() { cin >> N; dp[2][0][0][0] = 2; dp[1][0][0][0] = 1; for(int i = 2; i < N; i++) { for(int j = 0; j < N; j++) { for(int k = 0; k < 2; k++) { for(int l = 0; l < 2; l++) { if(dp[i][j][k][l] == 0) continue; if(k == 0 && l == 0) { ADD(dp[i + 1][j][0][1], dp[i][j][k][l] * 2); if(j > 0) ADD(dp[i + 1][j - 1][0][0], dp[i][j][k][l] * j); ADD(dp[i + 1][j][0][0], dp[i][j][k][l] * (i - j - 1)); } if(k == 1 && l == 0) { ADD(dp[i + 1][j][0][1], dp[i][j][k][l]); ADD(dp[i + 1][j + 1][0][1], dp[i][j][k][l]); if(j > 0) ADD(dp[i + 1][j][0][0], dp[i][j][k][l] * j); ADD(dp[i + 1][j + 1][0][0], dp[i][j][k][l] * (i - j - 1)); } if(k == 0 && l == 1) { ADD(dp[i + 1][j][1][1], dp[i][j][k][l] * 2); ADD(dp[i + 1][j][0][0], dp[i][j][k][l]); if(j > 0) ADD(dp[i + 1][j - 1][1][0], dp[i][j][k][l] * j); ADD(dp[i + 1][j][1][0], dp[i][j][k][l] * (i - j - 2)); } if(k == 1 && l == 1) { ADD(dp[i + 1][j][1][1], dp[i][j][k][l]); ADD(dp[i + 1][j + 1][0][0], dp[i][j][k][l]); ADD(dp[i + 1][j + 1][1][1], dp[i][j][k][l]); if(j > 0) ADD(dp[i + 1][j][1][0], dp[i][j][k][l] * j); ADD(dp[i + 1][j + 1][1][0], dp[i][j][k][l] * (i - j - 2)); } } } } } cout << dp[N][0][0][0] << endl; }