結果

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

ソースコード

diff #
raw source code 

#include<bits/stdc++.h>
#define int long long
using namespace std;
const int mo=1e9+7;
int n,f[2001][2001][2][2];
signed main(){
    cin>>n;
    f[1][0][0][0]=1,f[2][0][0][0]=2;
    for(int i=2;i<=n;i++){
        for(int j=0;j<=i;j++){
            for(int k=0;k<=1;k++){
                for(int l=0;l<=1;l++){
                    if(j==0&&(l==1||k==1)) continue;
                    if(k==1){
                        f[i+1][j][l][0]+=f[i][j][k][l]*(i+1-j-1);//???
                        f[i+1][j][l][0]%=mo;
                        f[i+1][j][l][1]+=f[i][j][k][l];//??-1?i-3??
                        f[i+1][j][l][1]%=mo;
                        f[i+1][j+1][l][1]+=f[i][j][k][l];//?i-1??
                        f[i+1][j+1][l][1]%=mo;
                        if(l==1){
                            f[i+1][j-1][l][0]+=f[i][j][k][l]*(j-2);
                            f[i+1][j-1][l][0]%=mo;
                            f[i+1][j-1][0][0]+=f[i][j][k][l];
                            f[i+1][j-1][0][0]%=mo;
                        }
                        else {
                            f[i+1][j-1][l][0]+=f[i][j][k][l]*(j-1);//?|p_x-p_x-1|=2??(???i-1?i-3)
                            f[i+1][j-1][l][0]%=mo;
                        }
                    }
                    else{
                        f[i+1][j][l][0]+=f[i][j][k][l]*(i+1-j-2);//???
                        f[i+1][j][l][0]%=mo;
                        f[i+1][j+1][l][1]+=f[i][j][k][l]*2;//?i-1??
                        f[i+1][j+1][l][1]%=mo;
                        if(j!=0)
                        if(l==1){
                            f[i+1][j-1][0][0]+=f[i][j][k][l];
                            f[i+1][j-1][0][0]%=mo;
                            f[i+1][j-1][l][0]+=f[i][j][k][l]*(j-1);
                            f[i+1][j-1][l][0]%=mo;
                        }
                        else {
                            f[i+1][j-1][l][0]+=f[i][j][k][l]*j;//?|p_x-p_x-1|=2??
                            f[i+1][j-1][l][0]%=mo;
                        }
                    }
                    //cout<<f[i][j][k][l]<<" ";
                }
            }
            //cout<<endl;
        }
        //cout<<endl;
    }
    cout<<f[n][0][0][0]<<endl;
    return 0;
}
0