[复习]动态dp

其他 2019-03-21 13:51:19 阅读次数: 0

[复习]动态dp

你还是可以认为我原来写的动态dp就是在扯蛋。

[Luogu4719]【模板】动态dp

首先作为一个\(dp\)题，我们显然可以每次修改之后都进行暴力\(dp\)，设\(f[i][0/1]\)表示当前考虑\(i\)及其子树内的点，当前这个点是选还是不选时能够得到的最大权值，那么我们可以得到转移：\(f[i][0]+=\max\{f[v][0],f[v][1]\},f[i][1]+=f[v][0]\)，其中\(v\)是\(i\)的一个儿子。
那么这样子的复杂度就是\(O(qn)\)。
仔细想想我们这样子为什么慢？因为我们重复计算了大量相同的转移，那么我们发现我们可以每次只需要修改当前的修改点到达根节点的这条链的所有点的\(dp\)值就好了，然而这样子复杂度最差还是\(O(nq)\)的。
那么我们还有哪里慢呢？
仔细想想，我们的转移都是一模一样的，区别只是在于每次的取值不同而已。
先考虑一模一样的转移可以怎么处理，似乎可以矩阵乘法？那么转移可以写成：
\[ \begin{bmatrix}0&0\\V_u&-\infty\end{bmatrix}\times \begin{bmatrix}f_{v,0}\\f_{v,1}\end{bmatrix}=\begin{bmatrix}f_{u,0}\\f_{u,1}\end{bmatrix} \]
然而这样子的转移只有在这个点转移第一个儿子的时候是对的。
如果有多个儿子的话，我们把\(V_u\)看成除了这个儿子之外的\(f[u][1]\),\(0\)看成除了这个儿子之外的\(f[u][0]\)。
那么等价于是我们要把儿子给分成两个部分，第一部分表示的是这个特殊转移的儿子，第二部分是剩下的儿子。似乎可以重链剖分？
那么对于重儿子我们可以直接维护这个矩阵，这样子用线段树维护矩阵乘法就可以算出重儿子部分的贡献，其他轻儿子暴力往上贡献。
看起来这个东西复杂度就很对了啊，单次修改只会造成\(log\)次暴力修改。

#include<iostream>
#include<cstdio>
using namespace std;
#define ll long long
#define MAX 100100
#define lson (now<<1)
#define rson (now<<1|1)
const ll inf=1e17;
inline int read()
{
    int x=0;bool t=false;char ch=getchar();
    while((ch<'0'||ch>'9')&&ch!='-')ch=getchar();
    if(ch=='-')t=true,ch=getchar();
    while(ch<='9'&&ch>='0')x=x*10+ch-48,ch=getchar();
    return t?-x:x;
}
struct Line{int v,next;}e[MAX<<1];
int h[MAX],cnt=1;
inline void Add(int u,int v){e[cnt]=(Line){v,h[u]};h[u]=cnt++;}
int n,Q,V[MAX];
int fa[MAX],top[MAX],bot[MAX],size[MAX],hson[MAX],dfn[MAX],tim,ln[MAX];
void dfs1(int u,int ff)
{
    fa[u]=ff;size[u]=1;
    for(int i=h[u];i;i=e[i].next)
    {
        int v=e[i].v;if(v==ff)continue;
        dfs1(v,u);size[u]+=size[v];
        if(size[v]>size[hson[u]])hson[u]=v;
    }
}
void dfs2(int u,int tp)
{
    top[u]=tp;dfn[u]=++tim,ln[tim]=u;
    if(hson[u])dfs2(hson[u],tp),bot[u]=bot[hson[u]];
    else bot[u]=u;
    for(int i=h[u];i;i=e[i].next)
        if(e[i].v!=fa[u]&&e[i].v!=hson[u])
            dfs2(e[i].v,e[i].v);
}
ll f[MAX][2];
void dfs(int u,int ff)
{
    f[u][0]=0;f[u][1]=V[u];
    for(int i=h[u];i;i=e[i].next)
    {
        int v=e[i].v;if(v==ff)continue;
        dfs(v,u);
        f[u][0]+=max(f[v][0],f[v][1]);
        f[u][1]+=f[v][0];
    }
}
struct Matrix
{
    ll a[2][2];
    ll*operator[](int x){return a[x];}
}t[MAX<<2],tmp[MAX];
Matrix operator*(Matrix a,Matrix b)
{
    Matrix c;
    c[0][0]=max(a[0][0]+b[0][0],a[0][1]+b[1][0]);
    c[0][1]=max(a[0][0]+b[0][1],a[0][1]+b[1][1]);
    c[1][0]=max(a[1][0]+b[0][0],a[1][1]+b[1][0]);
    c[1][1]=max(a[1][0]+b[0][1],a[1][1]+b[1][1]);
    return c;
}
void Build(int now,int l,int r)
{
    if(l==r)
    {
        int u=ln[l];ll g0=0,g1=V[u];
        for(int i=h[u];i;i=e[i].next)
            if(e[i].v!=fa[u]&&e[i].v!=hson[u])
                g0+=max(f[e[i].v][0],f[e[i].v][1]),g1+=f[e[i].v][0];
        tmp[l]=t[now]=(Matrix){g0,g0,g1,-inf};
        return;
    }
    int mid=(l+r)>>1;
    Build(lson,l,mid);Build(rson,mid+1,r);
    t[now]=t[lson]*t[rson];
}
void Modify(int now,int l,int r,int p)
{
    if(l==r){t[now]=tmp[l];return;}
    int mid=(l+r)>>1;
    if(p<=mid)Modify(lson,l,mid,p);
    else Modify(rson,mid+1,r,p);
    t[now]=t[lson]*t[rson];
}
Matrix Query(int now,int l,int r,int L,int R)
{
    if(L==l&&r==R)return t[now];
    int mid=(l+r)>>1;
    if(R<=mid)return Query(lson,l,mid,L,R);
    if(L>mid)return Query(rson,mid+1,r,L,R);
    return Query(lson,l,mid,L,mid)*Query(rson,mid+1,r,mid+1,R);
}
void Modify(int u,int w)
{
    tmp[dfn[u]][1][0]+=w-V[u];V[u]=w;
    while(u)
    {
        Matrix a=Query(1,1,n,dfn[top[u]],dfn[bot[u]]);
        Modify(1,1,n,dfn[u]);
        Matrix b=Query(1,1,n,dfn[top[u]],dfn[bot[u]]);
        u=fa[top[u]];if(!u)break;int x=dfn[u];
        ll g0=a[0][0],g1=a[1][0],f0=b[0][0],f1=b[1][0];
        tmp[x][0][0]=tmp[x][0][1]=tmp[x][0][0]+max(f0,f1)-max(g0,g1);
        tmp[x][1][0]=tmp[x][1][0]+f0-g0;
    }
}
int main()
{
    n=read();Q=read();
    for(int i=1;i<=n;++i)V[i]=read();
    for(int i=1,u,v;i<n;++i)u=read(),v=read(),Add(u,v),Add(v,u);
    dfs1(1,0);dfs2(1,1);dfs(1,0);Build(1,1,n);
    while(Q--)
    {
        int x=read(),w=read();Modify(x,w);
        Matrix ans=Query(1,1,n,dfn[1],dfn[bot[1]]);
        printf("%lld\n",max(ans[0][0],ans[1][0]));
    }
    return 0;
}

[NOIP2018]保卫王国

要求最小点覆盖。这个东西等于权值和减去最大独立集。
而最大独立集就是上面的那个东西。
考虑如何强制某个点必须选或者不选，如果必须选则将其权值减去一个\(inf\)，如果必须不选则将其权值加上一个\(inf\)。
然后就做完了。
cpp #include<iostream> #include<cstdio> using namespace std; #define ll long long #define MAX 100100 #define lson (now<<1) #define rson (now<<1|1) const ll inf=1e17; inline int read() { int x=0;bool t=false;char ch=getchar(); while((ch<'0'||ch>'9')&&ch!='-')ch=getchar(); if(ch=='-')t=true,ch=getchar(); while(ch<='9'&&ch>='0')x=x*10+ch-48,ch=getchar(); return t?-x:x; } struct Line{int v,next;}e[MAX<<1]; int h[MAX],cnt=1; inline void Add(int u,int v){e[cnt]=(Line){v,h[u]};h[u]=cnt++;} int n,Q;ll V[MAX]; int fa[MAX],top[MAX],bot[MAX],size[MAX],hson[MAX],dfn[MAX],tim,ln[MAX]; void dfs1(int u,int ff) { fa[u]=ff;size[u]=1; for(int i=h[u];i;i=e[i].next) { int v=e[i].v;if(v==ff)continue; dfs1(v,u);size[u]+=size[v]; if(size[v]>size[hson[u]])hson[u]=v; } } void dfs2(int u,int tp) { top[u]=tp;dfn[u]=++tim,ln[tim]=u; if(hson[u])dfs2(hson[u],tp),bot[u]=bot[hson[u]]; else bot[u]=u;#[复习]动态dp 你还是可以认为我原来写的动态dp就是在扯蛋。 ##[[Luogu4719]【模板】动态dp](https://www.luogu.org/problemnew/show/P4719) 首先作为一个$dp$题，我们显然可以每次修改之后都进行暴力$dp$，设$f[i][0/1]$表示当前考虑$i$及其子树内的点，当前这个点是选还是不选时能够得到的最大权值，那么我们可以得到转移：$f[i][0]+=\max\{f[v][0],f[v][1]\},f[i][1]+=f[v][0]$，其中$v$是$i$的一个儿子。 那么这样子的复杂度就是$O(qn)$。 仔细想想我们这样子为什么慢？因为我们重复计算了大量相同的转移，那么我们发现我们可以每次只需要修改当前的修改点到达根节点的这条链的所有点的$dp$值就好了，然而这样子复杂度最差还是$O(nq)$的。 那么我们还有哪里慢呢？ 仔细想想，我们的转移都是**一模一样**的，区别只是在于每次的取值不同而已。 先考虑一模一样的转移可以怎么处理，似乎可以矩阵乘法？那么转移可以写成： $$ \begin{bmatrix}0&0\\V_u&-\infty\end{bmatrix}\times \begin{bmatrix}f_{v,0}\\f_{v,1}\end{bmatrix}=\begin{bmatrix}f_{u,0}\\f_{u,1}\end{bmatrix} $$ 然而这样子的转移只有在这个点转移第一个儿子的时候是对的。 如果有多个儿子的话，我们把$V_u$看成除了这个儿子之外的$f[u][1]$,$0$看成除了这个儿子之外的$f[u][0]$。 那么等价于是我们要把儿子给分成两个部分，第一部分表示的是这个特殊转移的儿子，第二部分是剩下的儿子。似乎可以重链剖分？ 那么对于重儿子我们可以直接维护这个矩阵，这样子用线段树维护矩阵乘法就可以算出重儿子部分的贡献，其他轻儿子暴力往上贡献。 看起来这个东西复杂度就很对了啊，单次修改只会造成$log$次暴力修改。cpp

include

include

using namespace std;

define ll long long

define MAX 100100

define lson (now<<1)

define rson (now<<1|1)

const ll inf=1e17;
inline int read()
{
int x=0;bool t=false;char ch=getchar();
while((ch<'0'||ch>'9')&&ch!='-')ch=getchar();
if(ch=='-')t=true,ch=getchar();
while(ch<='9'&&ch>='0')x=x10+ch-48,ch=getchar();
return t?-x:x;
}
struct Line{int v,next;}e[MAX<<1];
int h[MAX],cnt=1;
inline void Add(int u,int v){e[cnt]=(Line){v,h[u]};h[u]=cnt++;}
int n,Q,V[MAX];
int fa[MAX],top[MAX],bot[MAX],size[MAX],hson[MAX],dfn[MAX],tim,ln[MAX];
void dfs1(int u,int ff)
{
fa[u]=ff;size[u]=1;
for(int i=h[u];i;i=e[i].next)
{
int v=e[i].v;if(v==ff)continue;
dfs1(v,u);size[u]+=size[v];
if(size[v]>size[hson[u]])hson[u]=v;
}
}
void dfs2(int u,int tp)
{
top[u]=tp;dfn[u]=++tim,ln[tim]=u;
if(hson[u])dfs2(hson[u],tp),bot[u]=bot[hson[u]];
else bot[u]=u;
for(int i=h[u];i;i=e[i].next)
if(e[i].v!=fa[u]&&e[i].v!=hson[u])
dfs2(e[i].v,e[i].v);
}
ll f[MAX][2];
void dfs(int u,int ff)
{
f[u][0]=0;f[u][1]=V[u];
for(int i=h[u];i;i=e[i].next)
{
int v=e[i].v;if(v==ff)continue;
dfs(v,u);
f[u][0]+=max(f[v][0],f[v][1]);
f[u][1]+=f[v][0];
}
}
struct Matrix
{
ll a[2][2];
lloperator{return a[x];}
}t[MAX<<2],tmp[MAX];
Matrix operator(Matrix a,Matrix b)
{
Matrix c;
c[0][0]=max(a[0][0]+b[0][0],a[0][1]+b[1][0]);
c[0][1]=max(a[0][0]+b[0][1],a[0][1]+b[1][1]);
c[1][0]=max(a[1][0]+b[0][0],a[1][1]+b[1][0]);
c[1][1]=max(a[1][0]+b[0][1],a[1][1]+b[1][1]);
return c;
}
void Build(int now,int l,int r)
{
if(l==r)
{
int u=ln[l];ll g0=0,g1=V[u];
for(int i=h[u];i;i=e[i].next)
if(e[i].v!=fa[u]&&e[i].v!=hson[u])
g0+=max(f[e[i].v][0],f[e[i].v][1]),g1+=f[e[i].v][0];
tmp[l]=t[now]=(Matrix){g0,g0,g1,-inf};
return;
}
int mid=(l+r)>>1;
Build(lson,l,mid);Build(rson,mid+1,r);
t[now]=t[lson]t[rson];
}
void Modify(int now,int l,int r,int p)
{
if(l==r){t[now]=tmp[l];return;}
int mid=(l+r)>>1;
if(p<=mid)Modify(lson,l,mid,p);
else Modify(rson,mid+1,r,p);
t[now]=t[lson]t[rson];
}
Matrix Query(int now,int l,int r,int L,int R)
{
if(L==l&&r==R)return t[now];
int mid=(l+r)>>1;
if(R<=mid)return Query(lson,l,mid,L,R);
if(L>mid)return Query(rson,mid+1,r,L,R);
return Query(lson,l,mid,L,mid)Query(rson,mid+1,r,mid+1,R);
}
void Modify(int u,int w)
{
tmp[dfn[u]][1][0]+=w-V[u];V[u]=w;
while(u)
{
Matrix a=Query(1,1,n,dfn[top[u]],dfn[bot[u]]);
Modify(1,1,n,dfn[u]);
Matrix b=Query(1,1,n,dfn[top[u]],dfn[bot[u]]);
u=fa[top[u]];if(!u)break;int x=dfn[u];
ll g0=a[0][0],g1=a[1][0],f0=b[0][0],f1=b[1][0];
tmp[x][0][0]=tmp[x][0][1]=tmp[x][0][0]+max(f0,f1)-max(g0,g1);
tmp[x][1][0]=tmp[x][1][0]+f0-g0;
}
}
int main()
{
n=read();Q=read();
for(int i=1;i<=n;++i)V[i]=read();
for(int i=1,u,v;i<n;++i)u=read(),v=read(),Add(u,v),Add(v,u);
dfs1(1,0);dfs2(1,1);dfs(1,0);Build(1,1,n);
while(Q--)
{
int x=read(),w=read();Modify(x,w);
Matrix ans=Query(1,1,n,dfn[1],dfn[bot[1]]);
printf("%lld\n",max(ans[0][0],ans[1][0]));
}
return 0;
}
##[[NOIP2018]保卫王国](https://www.luogu.org/problemnew/show/P5024) 要求最小点覆盖。这个东西等于权值和减去最大独立集。 而最大独立集就是上面的那个东西。 考虑如何强制某个点必须选或者不选，如果必须选则将其权值减去一个$inf$，如果必须不选则将其权值加上一个$inf$。 然后就做完了。cpp

include

include

using namespace std;

define ll long long

define MAX 100100

define lson (now<<1)

define rson (now<<1|1)

const ll inf=1e17;
inline int read()
{
int x=0;bool t=false;char ch=getchar();
while((ch<'0'||ch>'9')&&ch!='-')ch=getchar();
if(ch=='-')t=true,ch=getchar();
while(ch<='9'&&ch>='0')x=x10+ch-48,ch=getchar();
return t?-x:x;
}
struct Line{int v,next;}e[MAX<<1];
int h[MAX],cnt=1;
inline void Add(int u,int v){e[cnt]=(Line){v,h[u]};h[u]=cnt++;}
int n,Q;ll V[MAX];
int fa[MAX],top[MAX],bot[MAX],size[MAX],hson[MAX],dfn[MAX],tim,ln[MAX];
void dfs1(int u,int ff)
{
fa[u]=ff;size[u]=1;
for(int i=h[u];i;i=e[i].next)
{
int v=e[i].v;if(v==ff)continue;
dfs1(v,u);size[u]+=size[v];
if(size[v]>size[hson[u]])hson[u]=v;
}
}
void dfs2(int u,int tp)
{
top[u]=tp;dfn[u]=++tim,ln[tim]=u;
if(hson[u])dfs2(hson[u],tp),bot[u]=bot[hson[u]];
else bot[u]=u;
for(int i=h[u];i;i=e[i].next)
if(e[i].v!=fa[u]&&e[i].v!=hson[u])
dfs2(e[i].v,e[i].v);
}
ll f[MAX][2];
void dfs(int u,int ff)
{
f[u][0]=0;f[u][1]=V[u];
for(int i=h[u];i;i=e[i].next)
{
int v=e[i].v;if(v==ff)continue;
dfs(v,u);
f[u][0]+=max(f[v][0],f[v][1]);
f[u][1]+=f[v][0];
}
}
struct Matrix
{
ll a[2][2];
lloperator{return a[x];}
}t[MAX<<2],tmp[MAX];
Matrix operator(Matrix a,Matrix b)
{
Matrix c;
c[0][0]=max(a[0][0]+b[0][0],a[0][1]+b[1][0]);
c[0][1]=max(a[0][0]+b[0][1],a[0][1]+b[1][1]);
c[1][0]=max(a[1][0]+b[0][0],a[1][1]+b[1][0]);
c[1][1]=max(a[1][0]+b[0][1],a[1][1]+b[1][1]);
return c;
}
void Build(int now,int l,int r)
{
if(l==r)
{
int u=ln[l];ll g0=0,g1=V[u];
for(int i=h[u];i;i=e[i].next)
if(e[i].v!=fa[u]&&e[i].v!=hson[u])
g0+=max(f[e[i].v][0],f[e[i].v][1]),g1+=f[e[i].v][0];
tmp[l]=t[now]=(Matrix){g0,g0,g1,-inf};
return;
}
int mid=(l+r)>>1;
Build(lson,l,mid);Build(rson,mid+1,r);
t[now]=t[lson]t[rson];
}
void Modify(int now,int l,int r,int p)
{
if(l==r){t[now]=tmp[l];return;}
int mid=(l+r)>>1;
if(p<=mid)Modify(lson,l,mid,p);
else Modify(rson,mid+1,r,p);
t[now]=t[lson]t[rson];
}
Matrix Query(int now,int l,int r,int L,int R)
{
if(L==l&&r==R)return t[now];
int mid=(l+r)>>1;
if(R<=mid)return Query(lson,l,mid,L,R);
if(L>mid)return Query(rson,mid+1,r,L,R);
return Query(lson,l,mid,L,mid)Query(rson,mid+1,r,mid+1,R);
}
void Modify(int u,ll w)
{
tmp[dfn[u]][1][0]+=w-V[u];V[u]=w;
while(u)
{
Matrix a=Query(1,1,n,dfn[top[u]],dfn[bot[u]]);
Modify(1,1,n,dfn[u]);
Matrix b=Query(1,1,n,dfn[top[u]],dfn[bot[u]]);
u=fa[top[u]];if(!u)break;int x=dfn[u];
ll g0=a[0][0],g1=a[1][0],f0=b[0][0],f1=b[1][0];
tmp[x][0][0]=tmp[x][0][1]=tmp[x][0][0]+max(f0,f1)-max(g0,g1);
tmp[x][1][0]=tmp[x][1][0]+f0-g0;
}
}
int main()
{
n=read();Q=read();scanf("%*s");ll sum=0;
for(int i=1;i<=n;++i)sum+=(V[i]=read());
for(int i=1,u,v;i<n;++i)u=read(),v=read(),Add(u,v),Add(v,u);
dfs1(1,0);dfs2(1,1);dfs(1,0);Build(1,1,n);
while(Q--)
{
int a=read(),x=read(),b=read(),y=read(),pa=V[a],pb=V[b];
ll pls=0;
if(!x)Modify(a,pa+inf),pls+=inf;else Modify(a,pa-inf);
if(!y)Modify(b,pb+inf),pls+=inf;else Modify(b,pb-inf);
Matrix ret=Query(1,1,n,dfn[1],dfn[bot[1]]);
ll ans=max(ret[0][0],ret[1][0])-pls;
if(ans<0)puts("-1");
else printf("%lld\n",sum-ans);
Modify(a,pa);Modify(b,pb);
}
return 0;
}


##[[SDOI2017]切树游戏](https://www.luogu.org/problemnew/show/P3781)
[戳这里](https://www.cnblogs.com/cjyyb/p/10570636.html)
    for(int i=h[u];i;i=e[i].next)
        if(e[i].v!=fa[u]&&e[i].v!=hson[u])
            dfs2(e[i].v,e[i].v);
}
ll f[MAX][2];
void dfs(int u,int ff)
{
    f[u][0]=0;f[u][1]=V[u];
    for(int i=h[u];i;i=e[i].next)
    {
        int v=e[i].v;if(v==ff)continue;
        dfs(v,u);
        f[u][0]+=max(f[v][0],f[v][1]);
        f[u][1]+=f[v][0];
    }
}
struct Matrix
{
    ll a[2][2];
    ll*operator[](int x){return a[x];}
}t[MAX<<2],tmp[MAX];
Matrix operator*(Matrix a,Matrix b)
{
    Matrix c;
    c[0][0]=max(a[0][0]+b[0][0],a[0][1]+b[1][0]);
    c[0][1]=max(a[0][0]+b[0][1],a[0][1]+b[1][1]);
    c[1][0]=max(a[1][0]+b[0][0],a[1][1]+b[1][0]);
    c[1][1]=max(a[1][0]+b[0][1],a[1][1]+b[1][1]);
    return c;
}
void Build(int now,int l,int r)
{
    if(l==r)
    {
        int u=ln[l];ll g0=0,g1=V[u];
        for(int i=h[u];i;i=e[i].next)
            if(e[i].v!=fa[u]&&e[i].v!=hson[u])
                g0+=max(f[e[i].v][0],f[e[i].v][1]),g1+=f[e[i].v][0];
        tmp[l]=t[now]=(Matrix){g0,g0,g1,-inf};
        return;
    }
    int mid=(l+r)>>1;
    Build(lson,l,mid);Build(rson,mid+1,r);
    t[now]=t[lson]*t[rson];
}
void Modify(int now,int l,int r,int p)
{
    if(l==r){t[now]=tmp[l];return;}
    int mid=(l+r)>>1;
    if(p<=mid)Modify(lson,l,mid,p);
    else Modify(rson,mid+1,r,p);
    t[now]=t[lson]*t[rson];
}
Matrix Query(int now,int l,int r,int L,int R)
{
    if(L==l&&r==R)return t[now];
    int mid=(l+r)>>1;
    if(R<=mid)return Query(lson,l,mid,L,R);
    if(L>mid)return Query(rson,mid+1,r,L,R);
    return Query(lson,l,mid,L,mid)*Query(rson,mid+1,r,mid+1,R);
}
void Modify(int u,ll w)
{
    tmp[dfn[u]][1][0]+=w-V[u];V[u]=w;
    while(u)
    {
        Matrix a=Query(1,1,n,dfn[top[u]],dfn[bot[u]]);
        Modify(1,1,n,dfn[u]);
        Matrix b=Query(1,1,n,dfn[top[u]],dfn[bot[u]]);
        u=fa[top[u]];if(!u)break;int x=dfn[u];
        ll g0=a[0][0],g1=a[1][0],f0=b[0][0],f1=b[1][0];
        tmp[x][0][0]=tmp[x][0][1]=tmp[x][0][0]+max(f0,f1)-max(g0,g1);
        tmp[x][1][0]=tmp[x][1][0]+f0-g0;
    }
}
int main()
{
    n=read();Q=read();scanf("%*s");ll sum=0;
    for(int i=1;i<=n;++i)sum+=(V[i]=read());
    for(int i=1,u,v;i<n;++i)u=read(),v=read(),Add(u,v),Add(v,u);
    dfs1(1,0);dfs2(1,1);dfs(1,0);Build(1,1,n);
    while(Q--)
    {
        int a=read(),x=read(),b=read(),y=read(),pa=V[a],pb=V[b];
        ll pls=0;
        if(!x)Modify(a,pa+inf),pls+=inf;else Modify(a,pa-inf);
        if(!y)Modify(b,pb+inf),pls+=inf;else Modify(b,pb-inf);
        Matrix ret=Query(1,1,n,dfn[1],dfn[bot[1]]);
        ll ans=max(ret[0][0],ret[1][0])-pls;
        if(ans<0)puts("-1");
        else printf("%lld\n",sum-ans);
        Modify(a,pa);Modify(b,pb);
    }
    return 0;
}

[SDOI2017]切树游戏

戳这里

猜你喜欢

转载自www.cnblogs.com/cjyyb/p/10571095.html
今日推荐
数学建模Matlab之数据预处理方法
充电桩---ISO15118协议详细介绍
对话Kaldi之父、小米首席语音科学家Daniel Povey：开源环境比金钱和荣誉更吸引我 | AGI技术50人...
Hugging Face全攻略：轻松下载Llama 3模型，探索NLP的无限可能！【实操】
阅读送书抽奖？玩转抽奖游戏，js-tool-big-box工具库新上抽奖功能
百度发布Comate代码知识增强2.0，国内首个支持实时检索智能代码助手
黑客利用扫雷游戏 Python 克隆隐藏恶意脚本，攻击欧洲和美国金融机构
微软对开源字体 Cascadia Code 进行重大更新
好书推荐《ChatGPT原理与架构：大模型的预训练、迁移和中间件编程 》
Baidu Comate 智能编码助手：编程新伙伴，效率新飞跃
AI时代：人工智能大模型引领科技创造新时代
百篇博客 · 千里之行
周排行
WebSocket、HTTP 与 TCP
private,public,protected的区别
Python用了这么多年，总结出超实用的功能和特点
dgwp笔记
ModuleNotFoundError: No module named 'gdbm'
数组的去重方法
Ternsorflow 学习：005-MNIST 实现模型
SpringBoot 2 源码学习笔记（二）
jaxws-spring 搭建Web Services笔记
读取properties文件并获取属性值
每日归档
更多
2024-05-27(56)
2024-05-26(6)
2024-05-25(68)
2024-05-24(65)
2024-05-23(9)
2024-05-22(41)
2024-05-21(8)
2024-05-20(36)
2024-05-19(0)
2024-05-18(4)

代码天地 © 2018 codetd.com

境外服务器   最新电视剧2022