結果
| 問題 | 
                            No.93 ペガサス
                             | 
                    
| コンテスト | |
| ユーザー | 
                             vjudge1
                         | 
                    
| 提出日時 | 2025-09-06 12:58:55 | 
| 言語 | C++14  (gcc 13.3.0 + boost 1.87.0)  | 
                    
| 結果 | 
                             
                                AC
                                 
                             
                            
                         | 
                    
| 実行時間 | 19 ms / 5,000 ms | 
| コード長 | 896 bytes | 
| コンパイル時間 | 1,500 ms | 
| コンパイル使用メモリ | 162,120 KB | 
| 実行使用メモリ | 34,436 KB | 
| 最終ジャッジ日時 | 2025-09-06 12:58:58 | 
| 合計ジャッジ時間 | 2,786 ms | 
| 
                            ジャッジサーバーID (参考情報)  | 
                        judge3 / judge4 | 
(要ログイン)
| ファイルパターン | 結果 | 
|---|---|
| sample | AC * 4 | 
| other | AC * 16 | 
ソースコード
#include<bits/stdc++.h>
#define int long long
using namespace std;
const int N=1010,mod=1e9+7;
int n,f[N][N][2][2];
signed main(){
   ios::sync_with_stdio(0),cin.tie(0),cout.tie(0);
	cin>>n,f[1][0][0][0]=1,f[2][0][0][0]=2;
	for(int i=3;i<=n;i++){
		for(int j=0;j<i;j++){
			f[i][j][0][0]=((i-j-2)*f[i-1][j][0][0]%mod+(i-j-1)*f[i-1][j][0][1]%mod+(j+1)*f[i-1][j+1][0][0]%mod+j*f[i-1][j+1][0][1]+f[i-1][j+1][1][0]+f[i-1][j+1][1][1])%mod;
			f[i][j][0][1]=((i-j-2)*f[i-1][j][1][0]%mod+(i-j-1)*f[i-1][j][1][1]%mod+j*f[i-1][j+1][1][0]%mod+(j-1)*f[i-1][j+1][1][1]%mod)%mod;
			f[i][j][1][0]=(f[i-1][j][0][1]+(j?f[i-1][j-1][0][0]*2+f[i-1][j-1][0][1]:0)%mod)%mod;
			f[i][j][1][1]=(f[i-1][j][1][1]+(j?f[i-1][j-1][1][0]*2+f[i-1][j-1][1][1]:0)%mod)%mod;
		}
	}
	cout<<f[n][0][0][0]<<"\n";
	return 0;
}
/*
????? ??? ?
?????| ?_? _|
? ????`? _x ?
?? ? /??? ? |
??? /?  ??? ?
???|?? |?|?|
?| (???_?_)_)
????
*/
            
            
            
        
            
vjudge1