传送门


考虑对于每一个位置\(i\),计算所有可能的结果出现的概率。

定义一个区间\([l,r]\)为对于\(x\)的极大区间,当且仅当\(\max \limits _{i=l}^r \{a_i\} \leq x < min\{a_{l-1},a_{r+1}\}\)(其中\(a_0=a_{N+1} = MOD\))

那么我们需要考虑极大区间是如何缩小的。

考虑到设某一个点恰好为\(x\)的方案数转移还需容斥,不妨这么设状态:

设\(f_{x,p,l,r}\)表示已经进行了\(p\)次操作,且区间\([l,r]\)为对于\(x\)的极大区间的方案数。如果初始区间\([l,r]\)是对于\(x\)的极大区间,\(f_{x,0,l,r}=1\),否则为\(0\)。

考虑转移:

①第\(p\)次修改的区间对于\([l,r]\)是否为极大区间没有影响,也就是区间同时位于\([1,l-1]\)或者\([l,r]\)或者\([r+1,N]\),即\(f_{x,p,l,r} += f_{x,p-1,l,r} \times g_{l,r}\),其中\(g_{l,r} = C_l^2 + C_{r-l+2}^2 + C_{N-r+1}^2\)

②第\(p\)次修改的区间的左端点小于\(l\),右端点大于\(r\),这种直接会让\([l,r]\)变为非极大区间,贡献为\(0\)

③第\(p\)次修改的区间有一部分在之前的\([l,r]\),有一部分伸出去,即\(f_{x,p,l,r} += \sum\limits_{j=1}^{l-1} (j-1) \times f_{x,p-1,j,r} + \sum\limits_{j=r+1}^N (N-j) \times f_{x,p-1,l,j}\)

对于③显然是可以前缀和优化的,①中的\(g_{i,j}\)直接预处理,可以做到\(O(N^4)\)的\(DP\)

考虑计算答案,设第\(i\)个点的答案为\(ans_i\),则有\(ans_i = \sum\limits_{j =1}^{max} j \times (\sum\limits_{m=1} ^ i \sum\limits_{p=i}^N f_{j,Q,m,p} - \sum\limits_{m=1} ^ i \sum\limits_{p=i}^N f_{j-1,Q,m,p})\)

其实到这里似乎卡卡常就能过了,因为原题数据随机

但是\(O(N^4)\)过\(400\)未免太假,所以我们需要优化算法

我们把\(ans_i\)的式子拆开

考虑对\(\sum\limits_{m=1} ^ i \sum\limits_{p=i}^N f_{j,Q,m,p}\)进行合并,除了\(j=max\)时前面的系数为\(max\)以外,其他项前面的系数是\(-1\)。也就是说每一个\(\sum\limits_{m=1} ^ i \sum\limits_{p=i}^N f_{j,Q,m,p}\)对答案的贡献为\(-1\)。又考虑到上面\(DP\)的转移式子,\(x\)的值不会影响转移方程,所以我们可以将所有\(x\)的值一起转移

故设\(g_{i,l,r}\)表示进行了\(i\)次操作,区间\([l,r]\)对答案产生的贡献的总和,转移与上面\(f\)的转移是一样的

考虑初始值,对于一个区间\([l,r]\),如果它对于\(x\)是极大区间,则有\(\max \limits _{i=l}^r \{a_i\} \leq x < min\{a_{l-1},a_{r+1}\}\),也就是说合法的\(x\)共有\(min\{a_{l-1},a_{r+1}\}-\max \limits _{i=l}^r \{a_i\}\)个,而每一个\(f\)对答案的贡献为\(-1\),所以\(g_{0,l,r}=-(min\{a_{l-1},a_{r+1}\}-\max \limits _{i=l}^r \{a_i\})\)。

最后考虑一个细节:\(a_0=a_{N+1}=MOD\)。首先很显然地,\([1,N]\)需要是一个极大区间,否则会少算\(x=max\)的情况;然后对于\(f_{max,Q,i,j}\)来说,它的系数在\(ans\)中为\(max\),所以要求\(g_{0,1,N}=max\),也就是\(min\{a_0,a_{N+1}\}=MOD\),所以直接赋值为\(a_0=a_{N+1}=MOD\)即可。

总复杂度\(O(N^3)\),注意需要滚动数组压缩空间

#include<bits/stdc++.h>
//This code is written by Itst
using namespace std;

inline int read(){
    int a = 0;
    char c = getchar();
    bool f = 0;
    while(!isdigit(c) && c != EOF){
        if(c == '-')
            f = 1;
        c = getchar();
    }
    if(c == EOF)
        exit(0);
    while(isdigit(c)){
        a = (a << 3) + (a << 1) + (c ^ '0');
        c = getchar();
    }
    return f ? -a : a;
}

const int MOD = 1e9 + 7;
int num[410] , g[410][410] , dp[2][410][410] , C[410][410] , s1[2][410][410] , s2[2][410][410] , N , Q;

void input(){
    N = read();
    Q = read();
    for(int i = 1 ; i <= N ; ++i)
        num[i] = read();
    num[0] = num[N + 1] = 1e9 + 7;
}

inline int add(int a , int b){
    return a + b >= MOD ? a + b - MOD : a + b;
}

inline void calc(int x){
    for(int i = 1 ; i <= N ; ++i)
        for(int j = 1 ; j <= i ; ++j)
            s1[x][j][i] = add(s1[x][j - 1][i] , 1ll * (j - 1) * dp[x][j][i] % MOD);
    for(int i = N ; i ; --i)
        for(int j = N ; j >= i ; --j)
            s2[x][i][j] = add(s2[x][i][j + 1] , 1ll * (N - j) * dp[x][i][j] % MOD);
}

void init(){
    for(int i = 1 ; i <= N ; ++i)
        for(int j = i ; j <= N ; ++j){
            int maxN = 0;
            for(int k = i ; k <= j ; ++k)
                maxN = max(maxN , num[k]);
            if(maxN < min(num[i - 1] , num[j + 1]))
                dp[0][i][j] = maxN - min(num[i - 1] , num[j + 1]) + MOD;
        }
    for(int i = 0 ; i <= N + 1 ; ++i){
        C[i][0] = 1;
        for(int j = 1 ; j <= i ; ++j)
            C[i][j] = (C[i - 1][j - 1] + C[i - 1][j]) % MOD;
    }
    for(int i = 1 ; i <= N ; ++i)
        for(int j = i ; j <= N ; ++j)
            g[i][j] = add(C[i][2] , add(C[N - j + 1][2] , C[j - i + 2][2]));
    calc(0);
}

void work(){
    for(int i = 1 ; i <= Q ; ++i){
        for(int j = 1 ; j <= N ; ++j)
            for(int k = j ; k <= N ; ++k)
                dp[i & 1][j][k] = (1ll * dp[~i & 1][j][k] * g[j][k] + s1[~i & 1][j - 1][k] + s2[~i & 1][j][k + 1]) % MOD;
        calc(i & 1);
    }
}

void output(){
    for(int i = 1 ; i <= N ; ++i){
        int ans = 0;
        for(int j = 1 ; j <= i ; ++j)
            for(int k = i ; k <= N ; ++k)
                ans = (ans + dp[Q & 1][j][k]) % MOD;
        cout << ans << ' ';
    }
}

int main(){
#ifndef ONLINE_JUDGE
    freopen("in" , "r" , stdin);
    //freopen("out" , "w" , stdout);
#endif
    input();
    init();
    work();
    output();
    return 0;
}

Luogu3352 ZJOI2016 线段树 概率、区间DP的更多相关文章

  1. 二叉搜索树TREE(线段树,区间DP)

    前言 线段树+区间DP题,线段树却不是优化DP的,是不是很意外? 题面 二叉搜索树是一种二叉树,每个节点都有一个权值,并且一个点的权值比其左子树里的点权值都大,比起右子树里的点权值都小. 一种朴素的向 ...

  2. Codeforces Round #278 (Div. 1) Strip (线段树 二分 RMQ DP)

    Strip time limit per test 1 second memory limit per test 256 megabytes input standard input output s ...

  3. xdoj-1324 (区间离散化-线段树求区间最值)

    思想 : 1 优化:题意是覆盖点,将区间看成 (l,r)转化为( l-1,r) 覆盖区间 2 核心:dp[i]  覆盖从1到i区间的最小花费 dp[a[i].r]=min (dp[k])+a[i]s; ...

  4. Codeforces Round #222 (Div. 1) D. Developing Game 线段树有效区间合并

    D. Developing Game   Pavel is going to make a game of his dream. However, he knows that he can't mak ...

  5. hiho一下20周 线段树的区间修改

    线段树的区间修改 时间限制:10000ms 单点时限:1000ms 内存限制:256MB 描述 对于小Ho表现出的对线段树的理解,小Hi表示挺满意的,但是满意就够了么?于是小Hi将问题改了改,又出给了 ...

  6. 2016年湖南省第十二届大学生计算机程序设计竞赛---Parenthesis(线段树求区间最值)

    原题链接 http://acm.csu.edu.cn/OnlineJudge/problem.php?id=1809 Description Bobo has a balanced parenthes ...

  7. poj 3468:A Simple Problem with Integers(线段树,区间修改求和)

    A Simple Problem with Integers Time Limit: 5000MS   Memory Limit: 131072K Total Submissions: 58269   ...

  8. hdu 1556:Color the ball(线段树,区间更新,经典题)

    Color the ball Time Limit: 9000/3000 MS (Java/Others)    Memory Limit: 32768/32768 K (Java/Others)To ...

  9. hdu 1698:Just a Hook(线段树,区间更新)

    Just a Hook Time Limit: 4000/2000 MS (Java/Others)    Memory Limit: 32768/32768 K (Java/Others)Total ...

随机推荐

  1. TurboLinux系统管理习题一

    TurboLinux系统管理习题一 1. 使用vi编辑文本只读时,强制存盘并退出的命令是?(单选题)A :w!        B :q!       C   :wq!        D   :e!答案 ...

  2. ArcGIS JavaScript API动态图层

      矢量动态图层     <!DOCTYPE HTML> <html>   <head> <meta http-equiv="Content-Typ ...

  3. eclipse没有server选项解决方法

    eclipse是是一个开放源代码的.基于Java的可扩展开发平台.就其本身而言,它只是一个框架和一组服务,用于通过插件组件构建开发环境. 它使用频率十分高,然而当使用它配置weblogic的时候,经常 ...

  4. Jenkins 开启用户注册机制及用户权限设置

    Jenkins 开启用户注册机制及用户权限设置   by:授客 QQ:1033553122 步骤 1.  系统管理-Configure Global Security 2.  设置

  5. loadrunner 运行脚本-Run-time Settings之Preferences设置

    运行脚本-Run-time Settings之Preferences设置 by:授客 QQ:1033553122 打开Preferences设置对话框,这里提供了对运行时的参数选择设置 Enable ...

  6. [20180317]12c TABLE ACCESS BY INDEX ROWID BATCHED3.txt

    [20180317]12c TABLE ACCESS BY INDEX ROWID BATCHED3.txt --//简单探究12c TABLE ACCESS BY INDEX ROWID BATCH ...

  7. jmeter利用自身代理录制脚本

    在利用代理录制脚本时一定要安装java jdk,不然不能录制的. 没有安装过java jdk安装jmeter后打开时会提示安装jdk,但是mac系统中直接打开提示安装jdk页面后下载的java并不是j ...

  8. 自动化测试基础篇--Selenium Python环境搭建

    学习selenium python需要的工具: 1.浏览器 2.Python 3.Selenium 4.FireBug(Firefox) 5.chromedriver.IEDriverServer.g ...

  9. Chrome 百度搜索热点过滤插件 - 开源软件

    学习时,为了搜集最全的中文资料,有时候不得不使用Baidu搜索引擎.在你还是个小菜鸡的时候你可能会花费大量时间在百度上! 但是,时间久了你会发现,你总会被网络上一些奇奇怪怪或者有趣的事情吸引过去而逐渐 ...

  10. [MapReduce_add_5] MapReduce 实现标签的生成与聚合

    0. 说明 MapReduce 实现标签的生成与聚合 介绍 && 流程图 && 程序编写 1. 介绍 [1.1 原始有效数据] 86913510 {"revi ...