QOJ.ac

QOJ

ID题目提交者结果用时内存语言文件大小提交时间测评时间
#554949#9111. Zayin and StringCryingWA 418ms67532kbC++142.7kb2024-09-09 18:17:352024-09-09 18:17:35

Judging History

你现在查看的是最新测评结果

  • [2024-09-09 18:17:35]
  • 评测
  • 测评结果:WA
  • 用时:418ms
  • 内存:67532kb
  • [2024-09-09 18:17:35]
  • 提交

answer

#include<bits/stdc++.h>
using namespace std;
typedef long long ll; typedef double db;
const int N = 1010,M = 4010,K = 26,mod = 998244353;
const db eps = 1e-6;
int qpow(int a,int n){
    int res = 1; a%=mod;
    while(n){
        if(n&1)res = 1ll*res*a%mod;
        n>>=1; a = 1ll*a*a%mod;
    }
    return res;
}
int qinv(int a){return qpow(a,mod-2);}
//
int T,n,m; string s;

int ch[M][K],sum[M],fail[M],tot;

struct Info{
    db w; int cnt;
    Info(db vw = -1e18,int vc = 1e9){w = vw,cnt = vc;}
}; Info Z;
bool operator>(Info x,Info y){
    if(fabs(x.w-y.w) < eps)return x.cnt < y.cnt;
    else return x.w > y.w;
}
bool operator==(Info x,Info y){
    return fabs(x.w-y.w) < eps && x.cnt == y.cnt;
}

void build(){
    queue<int> qu;
    for(int i=0;i<K;i++)if(ch[0][i])qu.push(ch[0][i]);
    while(qu.size()){
        int u = qu.front(); qu.pop(); sum[u] += sum[fail[u]];
        for(int i=0;i<K;i++){
            if(ch[u][i])fail[ch[u][i]] = ch[fail[u]][i],qu.push(ch[u][i]);
            else ch[u][i] = ch[fail[u]][i];
        }
    }
}

Info dp[N][M];
Info calc(db w){
    for(int i=0;i<=n;i++)for(int j=0;j<=tot;j++)dp[i][j] = Z;
    for(int i=1;i<=n;i++){
        for(int j=0;j<=tot;j++){
            Info pre = dp[i-1][j]; if(pre == Z)continue;
            if(pre > dp[i][j])dp[i][j] = pre;
            int k = ch[j][s[i]-'a']; 
            pre.w += sum[k] + w; pre.cnt++;
            if(pre > dp[i][k])dp[i][k] = pre;
        }
        int k = ch[0][s[i]-'a']; Info now(sum[k]+w,1);
        if(now > dp[i][k])dp[i][k] = now;
    }
    Info res = Z;
    for(int j=0;j<=tot;j++)if(dp[n][j] > res)res = dp[n][j];
    return res;
}

void solve(){
    memset(ch,0,sizeof ch); memset(fail,0,sizeof fail); memset(sum,0,sizeof sum); tot = 0;

    cin>>n>>m>>s; s=" "+s;
    for(int i=1;i<=m;i++){
        string tmp; cin>>tmp; int w; cin>>w;
        int u = 0,len = tmp.length();
        for(int j=0;j<len;j++){
            int& v = ch[u][tmp[j]-'a'];
            if(!v)v = ++tot; u = v;
        }
        sum[u] += w;
    }
    build();
    //
    db L = 0,R = 1e11; int len = 0;
    while(R-L > eps){
        db mid = (L+R)/2;
        Info rev = calc(-mid);
        if(rev.w >= 0){
            len = rev.cnt; L = mid;
        }else{
            R = mid;
        }
    }
    //
    ll vL = 1,vR = 1e11,vres = 0;
    while(vL<=vR){
        ll mid = (vL+vR)>>1;
        db w = 1.00*mid/len;
        Info rev = calc(-w);
        if(rev.w >= 0){
            vres = mid,vL = mid+1;
        }else{
            vR = mid-1;
        }
    }
    //
    vres %= mod;
    int ans = 1ll*vres*qinv(len)%mod; cout<<ans<<"\n";
}

int main(){
    cin>>T;
    while(T--)solve();

    return 0;
}

詳細信息

Test #1:

score: 0
Wrong Answer
time: 418ms
memory: 67532kb

input:

80
4 7
ggeg
g 92946
d 65678
gg 50828
wrj 93954
gge 53780
a 94179
geg 40837
5 6
hiiii
ii 67984
foh 69185
hi 88153
pj 39431
i 32219
wfmve 96834
8 12
wvxxvwww
xw 1522
rzl 16262
wx 77596
v 69622
vw 82225
nykkmkv 19236
xv 83470
o 16781
w 2405
m 98319
ww 13889
qggbvty 16331
8 14
jjjiiijh
i 96670
z 74397
i...

output:

332874949
599030808
249640519
150055
665548146
81272
61572
67112
499196918
748779217
88888
831949361
74929
665552405
499139737
665543594
332830174
60785
748771076
63646
103574
101389
128437
39266
249703944
89874
110460
499215461
665540622
41750
782592345
96189
665612195
94537
83443
111657
665560098
...

result:

wrong answer 4th lines differ - expected: '332898173', found: '150055'