結果

問題 No.93 ペガサス
コンテスト
ユーザー vjudge1
提出日時 2025-09-08 11:17:05
言語 C++23
(gcc 13.3.0 + boost 1.89.0)
結果
AC  
実行時間 33 ms / 5,000 ms
コード長 1,021 bytes
記録
記録タグの例:
初AC ショートコード 純ショートコード 純主流ショートコード 最速実行時間
コンパイル時間 2,718 ms
コンパイル使用メモリ 275,244 KB
実行使用メモリ 23,040 KB
最終ジャッジ日時 2025-09-08 11:17:10
合計ジャッジ時間 3,871 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge4
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 4
other AC * 16
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

#include<bits/stdc++.h>
#define mm 1010
#define mod 1000000007
#define int long long
using namespace std;
int read()
{
	int x=0,f=1;
	char ch=getchar();
	while(ch<'0' || ch>'9')
	{
		if(ch=='-') f=-1;
		ch=getchar();
	}
	while(ch>='0' && ch<='9')
	{
		x=(x<<3)+(x<<1)+(ch^48);
		ch=getchar();
	}
	return x*f;
}
int n;
int f[mm][mm][2][2];
signed main()
{
	int n=read();
	f[1][0][0][0]=1,f[2][0][0][0]=2;
	for(int i=2;i<=n-1;i++)for(int j=0;j<i;j++)
	{
		if(j>0)
		{
			(f[i+1][j-1][0][0]+=f[i][j][0][0]*j+f[i][j][1][0]*(j-1)+f[i][j][0][1]+f[i][j][1][1])%=mod;
			(f[i+1][j][0][1]+=f[i][j][1][0])%=mod;
			(f[i+1][j-1][1][0]+=f[i][j][0][1]*(j-1)+f[i][j][1][1]*(j-2))%=mod;
		}
		(f[i+1][j][0][0]+=f[i][j][0][0]*(i-j-1)+f[i][j][1][0]*(i-j))%=mod;
		(f[i+1][j+1][0][1]+=f[i][j][0][0]*2+f[i][j][1][0])%=mod;
		(f[i+1][j][1][0]+=f[i][j][0][1]*(i-j-1)+f[i][j][1][1]*(i-j))%=mod;
		(f[i+1][j+1][1][1]+=f[i][j][0][1]*2+f[i][j][1][1])%=mod;
		(f[i+1][j][1][1]+=f[i][j][1][1])%=mod;
	}
	printf("%lld\n",f[n][0][0][0]);
	return 0;
}
0