結果

問題 No.599 回文かい
ユーザー momoyuumomoyuu
提出日時 2022-09-30 19:10:28
言語 C++23
(gcc 12.3.0 + boost 1.83.0)
結果
WA  
実行時間 -
コード長 4,228 bytes
コンパイル時間 2,658 ms
コンパイル使用メモリ 256,268 KB
実行使用メモリ 6,944 KB
最終ジャッジ日時 2024-06-02 02:43:20
合計ジャッジ時間 5,194 ms
ジャッジサーバーID
(参考情報)
judge4 / judge3
このコードへのチャレンジ
(要ログイン)

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
6,816 KB
testcase_01 AC 1 ms
6,944 KB
testcase_02 AC 2 ms
6,944 KB
testcase_03 AC 1 ms
6,944 KB
testcase_04 AC 1 ms
6,944 KB
testcase_05 WA -
testcase_06 AC 1 ms
6,940 KB
testcase_07 AC 2 ms
6,940 KB
testcase_08 AC 1 ms
6,940 KB
testcase_09 WA -
testcase_10 AC 3 ms
6,940 KB
testcase_11 AC 3 ms
6,940 KB
testcase_12 AC 4 ms
6,940 KB
testcase_13 WA -
testcase_14 WA -
testcase_15 WA -
testcase_16 WA -
testcase_17 WA -
testcase_18 AC 2 ms
6,944 KB
testcase_19 WA -
testcase_20 WA -
evil_0.txt WA -
権限があれば一括ダウンロードができます

ソースコード

diff #

#include <bits/stdc++.h>
using namespace std;

using ll = long long;
using ull = unsigned long long;
using ld = long double;
template<class t> using vc = vector<t>;
template<class t> using vvc = vc<vc<t>>;
using pi = pair<int,int>;
using pl = pair<ll,ll>;
using vi = vc<int>;
using vvi = vvc<int>;
using vl = vc<ll>;
using vvl = vvc<ll>;

#define rep(i,a,b) for (int i = (int)(a); i < (int)(b); i++)
#define irep(i,a,b) for (int i = (int)(a); i > (int)(b); i--)
#define all(a) a.begin(),a.end()
#define print(n) cout << n << '\n'
#define pritn(n) print(n)
#define printv(n,a) {copy(all(n),ostream_iterator<a>(cout," ")); cout<<"\n";}
#define printvv(n,a) {for(auto itr:n) printv(itr,a);}
#define rup(a,b) (a+b-1)/b
#define input(A,N) rep(i,0,N) cin>>A[i]
#define chmax(a,b) a = max(a,b)
#define chmin(a,b) a = min(a,b)

//const int mod = 998'244'353;
//const ll mod = 1'000'000'009;
//const ll mod = 67'280'421'310'721;
template< ll mod >
struct mint{
    long long x;
    mint(long long x=0):x((x%mod+mod)%mod){}
    mint operator-() const{
        return mint(-x);
    }
    mint& operator+=(const mint& a){
        if((x+=a.x)>=mod)x-=mod;
        return *this;
    }
    mint& operator-=(const mint& a){
        if((x+=mod-a.x)>=mod)x-=mod;
        return *this;
    }
    mint& operator*=(const  mint& a){
        (x *= a.x) %= mod;
        return *this;
    }
    mint operator+(const mint& a) const{
        mint res(*this);
        return res+=a;
    }
    mint operator-(const mint& a) const{
        mint res(*this);
        return res-=a;
    }
    mint operator*(const mint& a) const{
        mint res(*this);
        return res*=a;
    }
    mint pow(long long n) const {
        assert(0 <= n);
        mint a = *this, r = 1;
        while (n) {
            if (n & 1) r *= a;
            a *= a;
            n >>= 1;
        }
        return r;
    }
    mint inv() const{
        return pow(mod-2);
    }
    mint& operator/=(const mint& a){
        return (*this)*=a.inv();
    }
    mint operator/(const mint& a) const {
        mint res(*this);
        return res/=a;
    }
    friend ostream& operator<<(ostream& os, const mint& m){
        os << m.x;
        return os;
    }
    bool operator==(const mint& a) const {
        return x == a.x;
    }
    bool operator<(const mint& a) const{
        return x < a.x;
    }
};

//const ll B = 111;
//const ll m2 = 67'280'421'310'721;
//const ll m2 = 998244353;
//const ll m2 = (1LL<<61)-1LL; 

//https://gist.github.com/privet-kitty/295ac9202b7abb3039b493f8238bf40f#file-modulus-random-base-pair64-txt
struct rhash{
    string s;
    int n;
    ll mod,B;
    vc<ll> sum;
    vc<ll> powb;
    rhash(string s,ll mod,ll B):s(s),mod(mod),B(B){
        n = s.size();
        sum = vc<ll>(n+1,0);
        powb = vc<ll>(n+1,1);
        for(int i = 0;i<n;i++){
            sum[i+1] = sum[i];
            sum[i+1] *=  B;
            sum[i+1] %= mod;
            sum[i+1] += s[i];
            sum[i+1] %= mod;
            powb[i+1] = powb[i]*B;
            powb[i+1] %= mod;
        }
    }
    ll get(int l,int r){
        ll tmp = sum[r] - sum[l]*powb[r-l];
        tmp %= mod;
        if(tmp<0) tmp += mod;
        return tmp;
    }
};
const ll m1 = 1'000'000'009;

mint<m1> dfs(int l,int r,vi&vis,vc<rhash>&hash,vc<mint<m1>>&dp){
    if(vis[l])return dp[l];
    if(r<l) return 1;
    mint<m1> now(1);
    int L = l,R = r;
    while(L<R){
        bool p = true;
        rep(i,0,hash.size()){
            if(hash[i].get(l,L+1)==hash[i].get(R,r+1)) continue;
            p = false;
            break;
        }
        if(p){
            now += dfs(L+1,R-1,vis,hash,dp);
        }
        L++;R--;
    }
    vis[l] = 1;
    dp[l] = now;
    return now;
}

const ll m2 = 2147483647;
const ll B2 = 2147483634;
const ll m3 = 2147483629;
const ll B3 = 2147483627;
const ll m4 = 2147483587;
const ll B4 = 2147483583;


int main(){
    cout << fixed << setprecision(15);
    string s;
    cin>>s;
    int n = s.size();
    
    vc<rhash> hash;
    hash.push_back(rhash(s,m2,B2));
    hash.push_back(rhash(s,m3,B3));
    hash.push_back(rhash(s,m4,B4));
    
    vc<mint<m1>> dp(n,0);
    vi vis(n,0);
    print(dfs(0,n-1,vis,hash,dp));


    //system("pause");
    return 0;
}
0