結果

問題 No.93 ペガサス
ユーザー Yang33Yang33
提出日時 2018-04-22 03:37:24
言語 C++14
(gcc 13.2.0 + boost 1.83.0)
結果
AC  
実行時間 22 ms / 5,000 ms
コード長 2,589 bytes
コンパイル時間 1,514 ms
コンパイル使用メモリ 166,332 KB
実行使用メモリ 34,660 KB
最終ジャッジ日時 2023-09-09 12:28:04
合計ジャッジ時間 2,574 ms
ジャッジサーバーID
（参考情報） 		judge13 / judge11
このコードへのチャレンジ（β）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
4,380 KB
testcase_01 AC 2 ms
4,376 KB
testcase_02 AC 1 ms
4,380 KB
testcase_03 AC 2 ms
4,376 KB
testcase_04 AC 7 ms
14,632 KB
testcase_05 AC 4 ms
10,292 KB
testcase_06 AC 3 ms
8,168 KB
testcase_07 AC 15 ms
27,592 KB
testcase_08 AC 9 ms
19,028 KB
testcase_09 AC 21 ms
34,296 KB
testcase_10 AC 3 ms
5,876 KB
testcase_11 AC 3 ms
8,060 KB
testcase_12 AC 19 ms
32,152 KB
testcase_13 AC 8 ms
18,888 KB
testcase_14 AC 5 ms
12,560 KB
testcase_15 AC 17 ms
29,868 KB
testcase_16 AC 5 ms
12,508 KB
testcase_17 AC 2 ms
4,376 KB
testcase_18 AC 2 ms
4,376 KB
testcase_19 AC 22 ms
34,660 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
using namespace std;

using VS = vector<string>;    using LL = long long;
using VI = vector<int>;       using VVI = vector<VI>;
using PII = pair<int, int>;   using PLL = pair<LL, LL>;
using VL = vector<LL>;        using VVL = vector<VL>;

#define ALL(a)  begin((a)),end((a))
#define RALL(a) (a).rbegin(), (a).rend()
#define PB push_back
#define EB emplace_back
#define MP make_pair
#define SZ(a) int((a).size())
#define SORT(c) sort(ALL((c)))
#define RSORT(c) sort(RALL((c)))
#define UNIQ(c) (c).erase(unique(ALL((c))), end((c)))
#define FOR(i, s, e) for (int(i) = (s); (i) < (e); (i)++)
#define FORR(i, s, e) for (int(i) = (s); (i) > (e); (i)--)
#define debug(x) cerr << #x << ": " << x << endl
const int INF = 1e9;                          const LL LINF = 1e16;
const LL MOD = 1000000007;                    const double PI = acos(-1.0);
int DX[8] = { 0, 0, 1, -1, 1, 1, -1, -1 };    int DY[8] = { 1, -1, 0, 0, 1, -1, 1, -1 };

/* -----  2018/04/19  Problem: yukicoder 093  / Link: http://yukicoder.me/problems/no/093  ----- */
/* ------問題------



-----問題ここまで----- */
/* -----解説等-----

挿入DPのお気持ちはわかったけど

----解説ここまで---- */

LL N;

LL ans = 0LL;

LL dp[1003][1003][2][2];

int main() {
	cin.tie(0);
	ios_base::sync_with_stdio(false);

	cin >> N;
	dp[1][0][0][0] = 1;
	dp[2][0][0][0] = 2;
	FOR(i, 2, N) {
		FOR(j, 0, N) {


			{
				LL dp00 = dp[i][j][0][0];
				dp[i + 1][j + 1][0][1] += (2)*dp00; // o i-2 o
		   if(j)dp[i + 1][j - 1][0][0] += (j)*dp00; // [j]
				dp[i + 1][j    ][0][0] += (i - j - 1)*dp00; // other
			}
		    if(j){
				LL dp01 = dp[i][j][0][1];
				dp[i + 1][j - 1][0][0] += (1)*dp01; // i-1 o i-3
				dp[i + 1][j + 1][1][1] += (2)*dp01; // o i-2 o
				dp[i + 1][j - 1][1][0] += (j - 1)*dp01; // [j]
				dp[i + 1][j    ][1][0] += (i - j - 1)*dp01; // other 
			}
			if(j){

				LL dp10 = dp[i][j][1][0];
				dp[i + 1][j-1+1][0][1] += (1)*dp10; // i-2 o i-4
				dp[i + 1][j + 1][0][1] += (1)*dp10; // o i-2 i-4
				dp[i + 1][j - 1][0][0] += (j - 1)*dp10; // [j]
				dp[i + 1][j   ][0][0] += (i - j)*dp10; // other 
			}
			if(j){
				LL dp11 = dp[i][j][1][1];
				dp[i + 1][j-1+1][1][1] += (1)*dp11; // i-2 o i-4
				dp[i + 1][j + 1][1][1] += (1)*dp11; // o i-2 i-4
				dp[i + 1][j - 1][0][0] += (1)*dp11; // i-1 o i-3
				dp[i + 1][j - 1][1][0] += (j - 2)*dp11; // [j]
				dp[i + 1][j    ][1][0] += (i - j)*dp11; // other 
			}
		}
		// mod
		FOR(j, 0, N + 1)FOR(x, 0, 2)FOR(y, 0, 2)dp[i + 1][j][x][y] %= MOD;
	}




	ans = dp[N][0][0][0]%MOD;
	cout << ans << "\n";

	return 0;
}
0