結果

問題 No.93 ペガサス
コンテスト
ユーザー koyumeishi
提出日時 2015-09-09 16:21:26
言語 C++11（廃止可能性あり）
(gcc 13.3.0)
結果
AC  
実行時間 176 ms / 5,000 ms
コード長 2,114 bytes
コンパイル時間 732 ms
コンパイル使用メモリ 86,384 KB
実行使用メモリ 5,376 KB
最終ジャッジ日時 2024-07-19 05:04:28
合計ジャッジ時間 2,247 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge2
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 4
other AC * 16
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <iostream>
#include <vector>
#include <cstdio>
#include <sstream>
#include <map>
#include <string>
#include <algorithm>
#include <queue>
#include <cmath>
#include <set>
using namespace std;

#define MOD 1000000007

void resize_(long long& vec, vector<int>::iterator begin_, vector<int>::iterator end_, const long long& val){
	vec = val;
	return;
}

template<class T, class U>
void resize_(vector<T>& vec, vector<int>::iterator begin_, vector<int>::iterator end_, const U& val){
	vec.resize(*begin_);
	for(int i=0; i<vec.size(); i++){
		resize_(vec[i], begin_+1, end_, val);
	}
}

template<class T, class U>
void resize_(vector<T>& vec, vector<int>&& sz, const U& val){
	resize_(vec, sz.begin(), sz.end(), val);
}


int main(){
	int n;
	cin >> n;

	if(n==1){
		cout << 1 << endl;
		return 0;
	}else if(n==2){
		cout << 2 << endl;
		return 0;
	}

	vector<vector<vector<long long>>> dp;
	resize_(dp, {n+2, 2, 2}, 0);
	dp[0][0][0] = 2;

	for(int i=2; i<n; i++){
		vector<vector<vector<long long>>> dp_;
		resize_(dp_, {n+2,2,2}, 0);

		for(int k=0; k<i; k++)
		for(int a=0; a<2; a++)
		for(int b=0; b<2; b++){
			if(a && b){
				if(k>0) dp_[k-1][b][0] += dp[k][a][b] * (k-2);
						dp_[k  ][b][1] += dp[k][a][b] * 1;
						dp_[k+1][b][1] += dp[k][a][b] * 1;
				if(k>0) dp_[k-1][0][0] += dp[k][a][b] * 1;
						dp_[k  ][b][0] += dp[k][a][b] * (i+1-1-k);

			}else if(a && !b){
				if(k>0) dp_[k-1][b][0] += dp[k][a][b] * (k-1);
						dp_[k  ][b][1] += dp[k][a][b] * 1;
						dp_[k+1][b][1] += dp[k][a][b] * 1;
						dp_[k  ][b][0] += dp[k][a][b] * (i+1-1-k);

			}else if(!a && b){
				if(k>0) dp_[k-1][b][0] += dp[k][a][b] * (k-1);
						dp_[k+1][b][1] += dp[k][a][b] * 2;
				if(k>0) dp_[k-1][0][0] += dp[k][a][b] * 1;
						dp_[k  ][b][0] += dp[k][a][b] * (i+1-2-k);

			}else if(!a && !b){
				if(k>0) dp_[k-1][b][0] += dp[k][a][b] * k;
						dp_[k+1][b][1] += dp[k][a][b] * 2;
						dp_[k  ][b][0] += dp[k][a][b] * (i+1-2-k);
			}
		}

		for(int k=0; k<=i; k++)
		for(int a=0; a<2; a++)
		for(int b=0; b<2; b++){
			dp_[k][a][b] %= MOD;
		}
		swap(dp, dp_);
	}

	cout << dp[0][0][0] << endl;
	return 0;
}
0