結果

問題 No.599 回文かい
ユーザー confconf
提出日時 2017-03-18 01:12:05
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 709 ms / 4,000 ms
コード長 1,201 bytes
コンパイル時間 541 ms
コンパイル使用メモリ 65,872 KB
実行使用メモリ 6,944 KB
最終ジャッジ日時 2024-07-04 13:50:55
合計ジャッジ時間 5,279 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge5
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 22
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <iostream>
using namespace std;

int64_t DP[5002];

uint64_t mod_pow(uint64_t a, uint64_t e, uint64_t p){
    uint64_t res = 1;
    for(;e>0;e>>=1){
        if(e&1)res=(res*a)%p;
        a=(a*a)%p;
    }
    return res;
}

int main(){
    string S;
    cin>>S;
    int64_t ans = 1;
    DP[0]=1;
    uint64_t base=1007,mod=1000000009,hl,hr;
    for(int i=0;i<S.size()/2;i++){
        hl=0;hr=0;
        for(int j=i;j>=0;j--){
            hl=(hl*base+S[j])%mod;
            hr=(hr+mod_pow(base,i-j,mod)*S[S.size()-j-1])%mod;
            // cout<<S.substr(j,i-j+1)<<' '<<S.substr(S.size()-i-1,i-j+1)<<' '<<hl<<' '<<hr<<' '<<S[j]<<' '<<S[S.size()-j-1]<<endl;
            if(hl==hr){
                DP[i+1]+=DP[j];
                DP[i+1]%=1000000007;
            }
        }
        
        /*
        for(int j=0;j<=i;j++){
            cout<<S.substr(j,i-j+1)<<' '<<S.substr(S.size()-i-1,i-j+1)<<endl;
            if(S.substr(j,i-j+1)==S.substr(S.size()-i-1,i-j+1)){
                cout<<'A'<<endl;
                DP[i+1]+=DP[j];
                DP[i+1]%=1000000007;
            }
        }
        */
        ans+=DP[i+1];
        ans%=1000000007;
    }
    cout<<ans<<endl;
    return 0;
}
0