zoukankan      html  css  js  c++  java
  • 省选模拟赛day4

    怎么说?发现自己越来越菜了 到了不写题解写不出来题目的地步了。。

    这次题目我都有认真思考 尽管思考的时候状态不太好 但是 我想 再多给我时间也思考不出来什么吧

    所以写一份题解。

    T1 n个点的有根树 1号点为根 第i个点的点权为(w_i) (1leq w_ileq L) 定义一个长度为L的序列是完美的 当且仅当(forall 1leq ileq L)
    都有(w_{a_i}=i)(forall 1leq k<i) (a_k)(a_i)的父亲。定义(f_i)为第i个点的父亲.

    有m次修改 对于第i次修改会把(u=(i-1)%mod+n+1)的点权修改为(v_i)

    对于每次修改过后的答案乘上i 对1e9+7取模 询问m次修改的答案和。(Lleq nleq 1e6) (mleq 2e6) (1leq f_i<i)

    显然有一个nm的dp 只不过放到了树上。

    不难发现这个L很像我们的最长上升子序列 如果该序列中存在最长上升子序列的话 那这样求的就是最长上升子序列的方案数了。

    不过每次有单点修改。往这个方向思考 是无果的 基本上我们最长上升子序列方案数不支持修改 查询。

    但是可以发现 每次修改的点是有序的 1 2 3 4 5...n...1.... 同时看题目中关键的条件 (f_i<i)自己的父亲小于自己 也就是说对于每个点在修改的时候 自己的父亲已经被改过了 这是这道题的关键之处 修改并非无序修改。

    做题的时候抓住题目中的条件来解决问题 把题目中的条件都理解透了 再思考 问题的解决。

    考虑对于链的怎么做 其实就是序列了 容易设出设(f_i)表示以i为结尾的方案数 (g_i)表示以i为开头的方案数。

    一个点将会被修改掉 我们需要-其在原本的答案中所占方案数+新产生的方案。

    所占方案 那显然是(f_icdot g_i) 考虑一个新的数字(v_i) 其产生的方案为 ((sumlimits_{1leq k<i,w_k=v_i-1}f_k)cdot (sumlimits_{i+1leq kleq L,w_k==v_i+1}g_k))

    可以发现这正是我们新的f值 和g值的相乘。

    但是不可做的一点是 由于存在修改 我们一个点的f 和g容易得出 但是其他点的f和g的值将会被扰乱导致下一个点再这样做的时候结果不正确。

    当然这种修改是很难维护的 尽管我们维护值域线段树 想做出这样的修改也很不容易。

    接下来利用题目中的条件 自己的父亲比自己先修改这一条件。

    可以发现我们当前点在修改过后 对于新的f值我们利用父亲修改过后的再累计一个方案数维护。

    对于g的维护可以发现还是旧的就行 因为儿子不可能发生修改。

    但是有多次轮回修改怎么解决。其实可以直接莽了 每次树形dp处理一轮的问题。时间复杂度O(n) 同时也处理O(n)个问题。

    那么有O(m)个问题 显然 我们时间复杂度就是O(m)的 简单来讲就是 m/n轮处理 每次时间复杂度O(n)故总时间复杂度O(m).

    考虑在树上怎么做我们还是套用刚才的模型来做 对于f显然可以这样做。

    考虑g的维护 发现处理完一个子树的g 这个子树的g由于对我们的父亲还有贡献我们不能将其清理掉 所以可以进行桶差分一下即可 超级好写。

    综上我们解决了这个问题 时间复杂度O(m).(真的不是什么LIS的方案数啊啊啊。

    可以发现当自己认为问题不可被解决时 80%是自己没有认真读题。

    ···
    int n,m,L,len,top;
    int a[MAXN],v[MAXN];ll cf[MAXN],w[MAXN],b[MAXN],g[MAXN],cg[MAXN],f[MAXN],cnt,ans;
    int lin[MAXN],ver[MAXN],nex[MAXN];
    inline void add(int x,int y)
    {
    ver[++len]=y;
    nex[len]=lin[x];
    lin[x]=len;
    }
    inline void dp(int x)
    {
    ll ww=cg[v[x]+1];
    ll w1=cg[a[x]+1];
    //求出f数组 f数组表示以某个点为结尾的方案数 更新过的v数组来求出
    f[x]=cf[v[x]-1];
    if(v[x]1)f[x]=1;
    w[x]=cf[a[x]-1];
    if(a[x]
    1)w[x]=1;
    cf[v[x]]=(cf[v[x]]+f[x])%mod;
    go(x)dp(tn);
    cf[v[x]]=(cf[v[x]]-f[x])%mod;
    //求出g数组 g数组表示以某个点为开头的方案数 通过a数组求出
    g[x]=(cg[v[x]+1]-ww)%mod;
    if(v[x]L)g[x]=1;
    b[x]=(cg[a[x]+1]-w1)%mod;
    if(a[x]
    L)b[x]=1;
    cg[a[x]]=(cg[a[x]]+b[x])%mod;
    if(a[x]1)ans=(ans+b[x])%mod;
    }
    int main()
    {
    freopen("1.in","r",stdin);
    //freopen("tree.out","w",stdout);
    get(n);get(m);get(L);
    rep(2,n,i)add(read(),i);
    rep(1,n,i)v[i]=get(a[i]);
    int st;
    rep(1,m,i)
    {
    int x=(i-1)%n+1;
    if(x
    1)st=i;
    a[x]=v[x];
    get(v[x]);
    if(xn||im)
    {
    rep(1,n,j)cg[j]=cf[j]=0;
    rep(x+1,n,j)a[j]=v[j];
    ans=0;dp(1);
    rep(1,x,j)
    {
    ans=(ans-w[j]b[j]+f[j]g[j])%mod;
    cnt=(cnt+ans*st)%mod;++st;
    }
    }
    }
    printf("%lld ",(cnt+mod)%mod);
    return 0;
    }

    可以说这道题并非很难 失误之处在于没有抓住题目种给的条件。
  • 相关阅读:
    人脸识别算法初次了解
    白话经典算法系列之二 直接插入排序的三种实现
    常见浏览器兼容性问题与解决方式
    MP3的频率、比特率、码率与音质的关系
    程序猿接私活经验总结,来自csdn论坛语录
    Java虚拟机工作原理具体解释
    簡單SQL存儲過程實例
    全栈JavaScript之路(七)学习 Comment 类型节点.
    strdup函数的使用方法
    POJ 2823 Sliding Window 【单调队列】
  • 原文地址:https://www.cnblogs.com/chdy/p/12499080.html
Copyright © 2011-2022 走看看