[刷题#1][FJOI2015]所有公共子序列问题

题意简述:
给定两个长度分别为n,m的字符串X,Y以及一个整数k。
若k=0,输出X,Y的公共子序列个数
若k=1,输出X,Y的所有公共子序列以及个数,空集算在内。(字符只有英文大小写字母)
这一道题我刚开始没什么思路,然后看了一下题解,发现做法不是按照一位一位来,而是按照字符来建造Trie树。
于是自己思考了一下做法,以0位根节点,对于每个(i,j)(i<j)to[i][s[j]]=j,这样在dfs的时候省去了很多步骤,建造Trie的时间是O(n2)的,那么我们就可dfs了,由于k=1要全部输出,我们不用优化了,过不了还要打出题人呢。。。
k=0时,f[i][j]表示做到(i,j)的答案,记忆化搜索就可以过了。
在做这题的时候,我遇到了两个问题,随时都可能在比赛上置我于死地。
1.我将’a’~’z’,’A’~’Z’从小写到大写离散化,导致输出不会按照大小写大写优先的字典序,还好数据弱有80分,考场上写挂了真的太可怕了。。。
2.差点调吐血,我的下标从1开始,可是我输入的时候是scanf(“%s”,s)而不是scanf(“%s”,s+1)导致我调了1个小时,断点输出了好几段(代码里有)

代码:

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cstdlib>
int table[300],val[300];
struct trie{
    char s[3010];
    int size;
    int t[3010][60];
    void insert(int pos,int x){
        if(t[pos][table[s[x]]])return;
        t[pos][table[s[x]]]=x;
        //printf("pos = %d, x = %d, s[x] = %c, table[s[x]] = %d\n", pos, x, s[x], table[s[x]]);
    } 
    void build(){
        for(int i=0;i<=size;i++)
            for(int j=i+1;j<=size;j++)
                insert(i,j); 
    }
    void out(){
        for(int i=1;i<=size;i++)
            for(int j=1;j<=52;j++)
                if(t[i][j])printf("t[%d , %c] = %d\n",i,val[j],t[i][j]);
    }
    void scan(){
        scanf("%s",s+1);
    }
}X,Y;
int n,m,k;
char stk[3010];
long long f[3010][3010];
long long dfs(int i,int j,int top=0){
    //printf("now = (%d , %d)\n",i,j);
    if(!k&&~f[i][j]) return f[i][j];
    long long &res=f[i][j];res=1;
    if(k){
        for(int o=1;o<=top;o++)putchar(stk[o]);
        puts("");
    }
    top++;
    for(int c=1;c<=52;c++)
        if(X.t[i][c]&&Y.t[j][c])
    //      printf("next = (%d , %d)\n",X.t[i][c],Y.t[j][c]),
            stk[top]=val[c],res+=dfs(X.t[i][c],Y.t[j][c],top);
    return res;
}
int main(){
    memset(f,-1,sizeof f);
    for(int i='a';i<='z';i++)val[i-'a'+27]=i,table[i]=i-'a'+27;//,printf("%d : %c\n",table[i],i);
    for(int i='A';i<='Z';i++)val[i-'A'+1]=i,table[i]=i-'A'+1;//,printf("%d : %c\n",table[i],i);
    scanf("%d%d",&n,&m);
    X.scan(),Y.scan();
    X.size=n,Y.size=m;
    X.build(),Y.build();
    //X.out(),Y.out();
    scanf("%d",&k);
    printf("%lld\n",dfs(0,0));
    return 0;
}
为什么会从第三个点开始TLE? # P8868 [NOIP2022] 比赛 ## 目描述 小 N 和小 O 会在 2022 年 11 月参加一场盛大的程序设计大赛 NOIP!小 P 会作为裁判主持竞赛。小 N 和小 O 各自率领了一支 $n$ 个人的队伍,选手在每支队伍内都是从 $1$ 到 $n$ 编号。每一个选手都有相应的程序设计水平。具体的,小 N 率领的队伍中,编号为 $i$($1 \leq i \leq n$)的选手的程序设计水平为 $a _ i$;小 O 率领的队伍中,编号为 $i$($1 \leq i \leq n$)的选手的程序设计水平为 $b _ i$。特别地,$\{a _ i\}$ 和 $\{b _ i\}$ 还分别构成了从 $1$ 到 $n$ 的排列。 每场比赛前,考虑到路途距离,选手连续参加比赛等因素,小 P 会选择两个参数 $l, r$($1 \leq l \leq r \leq n$),表示这一场比赛会邀请两队中编号属于 $[l, r]$ 的所有选手来到现场准备比赛。在比赛现场,小 N 和小 O 会以掷骰子的方式挑选出参数 $p, q$($l \leq p \leq q \leq r$),只有编号属于 $[p, q]$ 的选手才能参赛。为了给观众以最精彩的比赛,两队都会派出编号在 $[p, q]$ 内的、程序设计水平值最大的选手参加比赛。假定小 N 派出的选手水平为 $m _ a$,小 O 派出的选手水平为 $m _ b$,则比赛的精彩程度为 $m _ a \times m _ b$。 NOIP 总共有 $Q$ 场比赛,每场比赛的参数 $l, r$ 都已经确定,但是 $p, q$ 还没有抽取。小 P 想知道,对于每一场比赛,在其所有可能的 $p, q$($l \leq p \leq q \leq r$)参数下的比赛的精彩程度之和。由于答案可能非常之大,你只需要对每一场答案输出结果对 $2 ^ {64}$ 取模的结果即可。 ## 输入格式 第一行包含两个正整数 $T, n$,分别表示测试点编号和参赛人数。如果数据为样例则保证 $T = 0$。 第二行包含 $n$ 个正整数,第 $i$ 个正整数为 $a _ i$,表示小 N 队伍中编号为 $i$ 的选手的程序设计水平。 第三行包含 $n$ 个正整数,第 $i$ 个正整数为 $b _ i$,表示小 O 队伍中编号为 $i$ 的选手的程序设计水平。 第四行包含一个正整数 $Q$,表示比赛场数。 接下来的 $Q$ 行,第 $i$ 行包含两个正整数 $l _ i, r _ i$,表示第 $i$ 场比赛的参数 $l, r$。 ## 输出格式 输出 $Q$ 行,第 $i$ 行包含一个非负整数,表示第 $i$ 场比赛中所有可能的比赛的精彩程度之和对 $2 ^ {64}$ 取模的结果。 ## 输入输出样例 #1 ### 输入 #1 ``` 0 2 2 1 1 2 1 1 2 ``` ### 输出 #1 ``` 8 ``` ## 输入输出样例 #2 ### 输入 #2 ``` 见附件下的 match/match2.in。 ``` ### 输出 #2 ``` 见附件下的 match/match2.ans。 ``` ## 输入输出样例 #3 ### 输入 #3 ``` 见附件下的 match/match3.in。 ``` ### 输出 #3 ``` 见附件下的 match/match3.ans。 ``` ## 说明/提示 **【样例 1 解释】** 当 $p = 1, q = 2$ 的时候,小 N 会派出 $1$ 号选手,小 O 会派出 $2$ 号选手,比赛精彩程度为 $2 \times 2 = 4$。 当 $p = 1, q = 1$ 的时候,小 N 会派出 $1$ 号选手,小 O 会派出 $1$ 号选手,比赛精彩程度为 $2 \times 1 = 2$。 当 $p = 2, q = 2$ 的时候,小 N 会派出 $2$ 号选手,小 O 会派出 $2$ 号选手,比赛精彩程度为 $1 \times 2 = 2$。 **【样例 2】** 该样例满足测试点 $1 \sim 2$ 的限制。 **【样例 3】** 该样例满足测试点 $3 \sim 5$ 的限制。 **【数据范围】** 对于所有数据,保证:$1 \leq n, Q \leq 2.5 \times 10 ^ 5$,$1 \leq l _ i \leq r _ i \leq n$,$1 \leq a _ i, b _ i \leq n$ 且 $\{a _ i\}$ 和 $\{b _ i\}$ 分别构成了从 $1$ 到 $n$ 的排列。 ::cute-table{tuack} | 测试点 | $n$ | $Q$ | 特殊性质 A | 特殊性质 B | | :----------: | :----------: | :----------: | :----------: | :----------: | | $1, 2$ | $\leq 30$ | $\leq 30$ | 是 | 是 | | $3, 4, 5$ | $\leq 3,000$ | $\leq 3,000$ | ^ | ^ | | $6, 7$ | $\leq 10 ^ 5$ | $\leq 5$ | ^ | ^ | | $8, 9$ | $\leq 2.5 \times 10 ^ 5$ | ^ | ^ | ^ | | $10, 11$ | $\leq 10 ^ 5$ | ^ | 否 | 否 | | $12, 13$ | $\leq 2.5 \times 10 ^ 5$ | ^ | ^ | ^ | | $14, 15$ | $\leq 10 ^ 5$ | $\leq 10 ^ 5$ | 是 | 是 | | $16, 17$ | $\leq 2.5 \times 10 ^ 5$ | $\leq 2.5 \times 10 ^ 5$ | ^ | ^ | | $18, 19$ | $\leq 10 ^ 5$ | $\leq 10 ^ 5$ | ^ | 否 | | $20, 21$ | $\leq 2.5 \times 10 ^ 5$ | $\leq 2.5 \times 10 ^ 5$ | ^ | ^ | | $22, 23$ | $\leq 10 ^ 5$ | $\leq 10 ^ 5$ | 否 | ^ | | $24, 25$ | $\leq 2.5 \times 10 ^ 5$ | $\leq 2.5 \times 10 ^ 5$ | ^ | ^ | 特殊性质 A:保证 $a$ 是均匀随机生成的 $1 \sim n$ 的排列。 特殊性质 B:保证 $b$ 是均匀随机生成的 $1 \sim n$ 的排列。 #include<bits/stdc++.h> using namespace std; int a[250005],b[250005]; struct node{ int l,r; int maxn; }; node ta[250005*4]; node tb[250005*4]; void builda(int i,int le,int ri){ ta[i].l=le; ta[i].r=ri; if(le==ri){ ta[i].maxn=a[le]; return; } int mid=(ri-le)/2+le; builda(i*2,le,mid); builda(i*2+1,mid+1,ri); ta[i].maxn=max(ta[i*2].maxn,ta[i*2+1].maxn); } void buildb(int i,int le,int ri){ tb[i].l=le; tb[i].r=ri; if(le==ri){ tb[i].maxn=b[le]; return; } int mid=(ri-le)/2+le; buildb(i*2,le,mid); buildb(i*2+1,mid+1,ri); tb[i].maxn=max(tb[i*2].maxn,tb[i*2+1].maxn); } long long querya(int i, int le, int ri){ if(ta[i].l>=le&&ta[i].r<=ri){ return ta[i].maxn; } int mid=(ta[i].l+ta[i].r)/2; if(ri<=mid){ return querya(i*2,le,ri); } else if(le>mid){ return querya(i*2+1,le,ri); } else{ return max(querya(i*2,le,ri),querya(i*2+1,le,ri)); } } long long queryb(int i, int le, int ri){ if(tb[i].l>=le&&tb[i].r<=ri){ return tb[i].maxn; } int mid=(tb[i].l+tb[i].r)/2; if(ri<=mid){ return queryb(i*2,le,ri); } else if(le>mid){ return queryb(i*2+1,le,ri); } else{ return max(queryb(i*2,le,ri),queryb(i*2+1,le,ri)); } } int main(){ int t,n; cin>>t>>n; for(int i=1;i<=n;i++){ cin>>a[i]; } for(int i=1;i<=n;i++){ cin>>b[i]; } builda(1,1,n); buildb(1,1,n); int q; cin>>q; while(q--){ int l,r; cin>>l>>r; long long ans=0; for(int i=l;i<=r;i++){ for(int j=i;j<=r;j++){ long long ma=querya(1,i,j); long long mb=queryb(1,i,j); ans=(ans+ma*mb); } } cout<<ans<<endl; } return 0; }
最新发布
11-27
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值