bzoj 1014 [JSOI2008]火星人prefix

填一波坑。。当初懒得挑这题。。(还是我太弱。。写的时候脑子不好使)。

这题很裸。。但我的错误太傻了。。

首先不要相信什么长度一定小于100000,不开二倍会死的。。。

然后每次增加点的时候要把全长+1.。。。(我就错这了。。调了1h)。

之后就是很裸的啦。

维护hash值,正常的splay,对于增加点,就是把前驱转上去,后继转到前驱的右子节点,新的点插入左子节点。

修改就是正常splay上去修改就好。

查询的时候利用二分,对于每一个答案进行验证就好了。

#include<cstdio>
#include<cmath>
#include<algorithm>
#include<cstring>
#define seed 13131
#define N 100005
using namespace std;
typedef unsigned long long ull;
struct node
{
    int son[2],siz,fa;
    ull has;
    char val;
}tr[2*N];
ull cj[2*N];
char ss[2*N],ins[10];
int root,tot,n,m;
void update(int x)
{
    tr[x].siz=tr[tr[x].son[0]].siz+tr[tr[x].son[1]].siz+1;
    tr[x].has=tr[tr[x].son[0]].has*cj[tr[tr[x].son[1]].siz+1]+cj[tr[tr[x].son[1]].siz]*tr[x].val+tr[tr[x].son[1]].has;
}
void rotate(int x)
{
    int y=tr[x].fa,z=tr[y].fa;
    int typ=(x==tr[y].son[1]);
    tr[y].son[typ]=tr[x].son[typ^1],tr[tr[x].son[typ^1]].fa=y;
    tr[x].son[typ^1]=y,tr[y].fa=x;
    tr[x].fa=z,tr[z].son[tr[z].son[1]==y]=x;
    update(y);update(x);
}
int findkth(int x,int k)
{
    int t=tr[tr[x].son[0]].siz;
    if(k<=t)return findkth(tr[x].son[0],k);
    else if(k==t+1)return x;
    else return findkth(tr[x].son[1],k-t-1);
}
void splay(int x,int goal)
{
    for(int y;(y=tr[x].fa)!=goal;rotate(x))
    {
        if(tr[y].fa!=goal)rotate((x==tr[y].son[0])==(y==tr[tr[y].fa].son[0])?y:x);
    }
    if(!goal)root=x;
}
int biuld(int l,int r,int f)
{
    if(l>r)return 0;
    int mid=(l+r)/2;
    int rt=++tot;
    tr[rt].val=ss[mid];tr[rt].siz=1;tr[rt].fa=f;//has会自己更新进去
 	tr[rt].son[0]=biuld(l,mid-1,rt);
 	tr[rt].son[1]=biuld(mid+1,r,rt);
 	update(rt);
 	return rt;
}
void insert(int po,char s)//x
{
    int x=findkth(root,po+1);
    int y=findkth(root,po+2);
    splay(x,0);splay(y,x);
    int rt=++tot;n++;
    tr[y].son[0]=rt;
    tr[rt].val=s;tr[rt].siz=1;tr[rt].fa=y;/*tr[rt].has=s*/
    update(rt);update(y);update(x);
}
void change(int k,char v)
{
    int x=findkth(root,k+1);
    splay(x,0);
    tr[x].val=v;
    update(x);
}
int LCQ(int xx,int yy)//x->x,x+mid+2
{
    int l=0,r=tot-2,mid;
    while(l<=r)
    {
        mid=(l+r)/2;
        if(yy+mid-1>n)
        {
            r=mid-1;continue;
        }
        int x=findkth(root,xx),y=findkth(root,xx+mid+1);
        splay(x,0);splay(y,x);
        ull tt=tr[tr[y].son[0]].has;
        x=findkth(root,yy),y=findkth(root,yy+mid+1);
        splay(x,0);splay(y,x);
        if(tt==tr[tr[y].son[0]].has)l=mid+1;
        else r=mid-1;
    }
    return r;
}
int main()
{
    scanf("%s",ss+2);n=strlen(ss+2);scanf("%d",&m);
    cj[0]=1;
    for(int i=1;i<=2*N-2;i++)
    {
        cj[i]=cj[i-1]*seed;
    }
    root=biuld(1,n+2,0);
    int t1,t2;char t3[2];
    for(int i=1;i<=m;i++)
    {
        scanf("%s",ins);
        if(ins[0]=='Q')
        {
            scanf("%d%d",&t1,&t2);
            if(t1>t2)swap(t1,t2);
        //	if(t1==t2)printf("%d\n",n-t1+1);
            printf("%d\n",LCQ(t1,t2));
        }else if(ins[0]=='R')
        {
            scanf("%d%s",&t1,t3);change(t1,t3[0]);
        }else 
        {
            scanf("%d%s",&t1,t3);insert(t1,t3[0]);
        }
    }
}

 

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
题目描述 有一个 $n$ 个点的棋盘,每个点上有一个数字 $a_i$,你需要从 $(1,1)$ 走到 $(n,n)$,每次只能往右或往下走,每个格子只能经过一次,路径上的数字和为 $S$。定义一个点 $(x,y)$ 的权值为 $a_x+a_y$,求所有满足条件的路径中,所有点的权值和的最小值。 输入格式 第一行一个整数 $n$。 接下来 $n$ 行,每行 $n$ 个整数,表示棋盘上每个点的数字。 输出格式 输出一个整数,表示所有满足条件的路径中,所有点的权值和的最小值。 数据范围 $1\leq n\leq 300$ 输入样例 3 1 2 3 4 5 6 7 8 9 输出样例 25 算法1 (树形dp) $O(n^3)$ 我们可以先将所有点的权值求出来,然后将其看作是一个有权值的图,问题就转化为了在这个图中求从 $(1,1)$ 到 $(n,n)$ 的所有路径中,所有点的权值和的最小值。 我们可以使用树形dp来解决这个问题,具体来说,我们可以将这个图看作是一棵树,每个点的父节点是它的前驱或者后继,然后我们从根节点开始,依次向下遍历,对于每个节点,我们可以考虑它的两个儿子,如果它的两个儿子都被遍历过了,那么我们就可以计算出从它的左儿子到它的右儿子的路径中,所有点的权值和的最小值,然后再将这个值加上当前节点的权值,就可以得到从根节点到当前节点的路径中,所有点的权值和的最小值。 时间复杂度 树形dp的时间复杂度是 $O(n^3)$。 C++ 代码 算法2 (动态规划) $O(n^3)$ 我们可以使用动态规划来解决这个问题,具体来说,我们可以定义 $f(i,j,s)$ 表示从 $(1,1)$ 到 $(i,j)$ 的所有路径中,所有点的权值和为 $s$ 的最小值,那么我们就可以得到如下的状态转移方程: $$ f(i,j,s)=\min\{f(i-1,j,s-a_{i,j}),f(i,j-1,s-a_{i,j})\} $$ 其中 $a_{i,j}$ 表示点 $(i,j)$ 的权值。 时间复杂度 动态规划的时间复杂度是 $O(n^3)$。 C++ 代码

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值