【AtCoder】【DP】【思维】Prefix Median(AGC012)
模的是这位神犇的代码:Atcoder AGC012F : Prefix Median
题意:
在动态中位数那道题上做了一些改动。给你一个序列a,可以将a重新任意排序,然后对于a序列构造出b序列。
假设a序列有2*n-1个元素,b序列有n个元素。
其中b[i]=Median(a[1],a[2],a[3]...a[2i-1])。求能够构造出多少个不同的b序列。
数据范围:
1<=N<=50,1<=ai<=2N-1
思路:
这道题真的是究极神题...虽然说代码实现比较简单,但是分析的过程是恶心到吐的。
首先,我们需要分析一堆性质:
首先将a序列排序,便于讨论
- \(a[i]<=b[i]<=a[2*N-i]\)
- \(当i<j时,不存在i使得:b[j]<b[i]<b[j+1]或者是b[j+1]<=b[i]<=b[j]\)
先抛开性质不谈,首先有很明显的一个性质:i从左往右走的时候,假设已经选取的元素的集合为c(c是排好了序的)。每加入两个元素,b[i]的取值的下标在c中只会向左或者是向右移动一格,或者是保持不变。(可以通过枚举法简单证得,即新加入的两个元素x,y都小于b[i-1],或者是都大于b[i+1],或者是一个大于,一个小于的情况)。
然后来看上面的性质。因为求b[i]时的序列中小于等于b[i]至少有i个元素,大于等于b[i]的至少有i个元素,所以说在a数组中的b[i]至少在a[i]~a[2N-i]之间。性质1是没有问题的;
然后是性质2。因为最开始所说的“明显的性质”,所以说b[j]与b[j+1]在已经选取的数字的序列之的位置相差不会超过1,也就是说中间不会间隔任何的数字。但同时还需要注意一点,就是这里是'<'和'>',也就是说是可以取等的,也就是位置保持不变的情况。那么性质2也有了。
接下来就是怎么实现了。
由于开头赋的题解对于代码的部分的具体细节解释的不是很详细,我这里就贸然做一些补充了(〃'▽'〃)
首先定义状态dp[i][L][R]表示现在确定的是b[i],可选元素中小于等于b[i+1]的有L个,大于b[i+1]的有R个。然后考虑转移。转移的话就是枚举b[i]选取的位置,然后删去相应的元素就可以了。然而并没有这么简单!!!具体细节将在代码中进行详细的说明。
代码:
#include<cstdio>
#include<cstring>
#include<algorithm>
#define MAXN 100
#define MO 1000000007
using namespace std;
typedef long long LL;
int a[MAXN+5],n;
LL f[MAXN+5][MAXN+5][MAXN+5];
int main()
{
scanf("%d",&n);
for(int i=1;i<=2*n-1;i++)
scanf("%d",&a[i]);
sort(a+1,a+2*n);
f[n][1][0]=1;
//初始状态中,小于等于b[n]的有1个,大于b[n]的有0个(虽然稍微有些违背定义,但是暂且先这样定义吧)
LL t;
for(int i=n-1;i>=1;i--)
{
int al=(a[i]!=a[i+1]),ar=(a[2*n-i]!=+a[2*n-i-1]);
//这里是向两边进行扩展,也就是运用性质1进行扩展
//因为相同的数进行的扩展是没有意义的,所以说需要这样进行能否进行有效的扩展的判断
for(int l=0;l<=2*n-1;l++)
for(int r=0;l+r<=2*n-1;r++)
if(f[i+1][l][r])
{
t=f[i+1][l][r];
for(int dl=1;dl<=l+al;dl++)//选择小于b[i+1]的元素
{
f[i][l+al-dl+1][r+ar+(dl>1)]+=t;
//+al是进行对于L的拓展
//-dl是删去选了的元素与b[i+1]之间的元素
//+1是因为删多了一个,那就是已经选取的b[i],所以说要加回来
//右边+(dl>1)是因为:
//当dl==1的时候,选取的就是b[i+1]这个元素,那就是什么都不会改变的
//当dl>1的时候,选取的就是<b[i+1]中的一个元素,那么b[i+1]就会成为>b[i]中的一个元素,又因为可以取整,所以说要加上去
f[i][l+al-dl+1][r+ar+(dl>1)]%=MO;
}
for(int dr=1;dr<=r+ar;dr++)//选择大于等于b[i+2]的元素
{
f[i][l+al+1][r+ar-dr]+=t;
//左边+1是因为加入b[i]这个可选元素,与上面的比较类似
//右边就是正常的转移,这还比较简单
f[i][l+al+1][r+ar-dr]%=MO;
}
}
}
LL ans=0;
for(int l=0;l<=2*n-1;l++)
for(int r=0;l+r<=2*n-1;r++)
ans=(1LL*ans+1LL*f[1][l][r])%MO;//对于每一个可能的情况都需要计算答案
printf("%lld\n",ans);
return 0;
}
以上只是个人的理解,如果出了什么偏差,还请读者自行脑补ヽ(・ω・´メ)
【AtCoder】【DP】【思维】Prefix Median(AGC012)的更多相关文章
- Atcoder Grand Contest 024 E - Sequence Growing Hard(dp+思维)
题目传送门 典型的 Atcoder 风格的计数 dp. 题目可以转化为每次在序列中插入一个 \([1,k]\) 的数,共操作 \(n\) 次,满足后一个序列的字典序严格大于前一个序列,问有多少种操作序 ...
- cf1153D 树形dp+思维
一千八的题也不会做了呜呜呜 size[u]表示结点u下的叶子结点, 思维:可以想到一个子树对其父亲会有一个消耗值 考虑一个点如果是max,那么其最大值可以是size[u]-p,p是消耗值最小的子树 一 ...
- E. The Contest ( 简单DP || 思维 + 贪心)
传送门 题意: 有 n 个数 (1 ~ n) 分给了三个人 a, b, c: 其中 a 有 k1 个, b 有 k2 个, c 有 k3 个. 现在问最少需要多少操作,使得 a 中所有数 是 1 ~ ...
- 7月15日考试 题解(链表+状压DP+思维题)
前言:蒟蒻太弱了,全打的暴力QAQ. --------------------- T1 小Z的求和 题目大意:求$\sum\limits_{i=1}^n \sum\limits_{j=i}^n kth ...
- codeforces 1140D(区间dp/思维题)
D. Minimum Triangulation time limit per test 2 seconds memory limit per test 256 megabytes input sta ...
- POJ 1390 Blocks(DP + 思维)题解
题意:有一排颜色的球,每次选择一个球消去,那么这个球所在的同颜色的整段都消去(和消消乐同理),若消去k个,那么得分k*k,问你消完所有球最大得分 思路:显然这里我们直接用二位数组设区间DP行不通,我们 ...
- “玲珑杯”ACM比赛 Round #18---图论你先敲完模板(DP+思维)
题目链接 DESCRIPTION INPUT OUTPUT SAMPLE INPUT 2 3 2 3 5 7 3 10 3 5 7 SAMPLE OUTPUT 12 26 HINT 官方题解: 代码如 ...
- HDU - 5117 Fluorescent(状压dp+思维)
原题链接 题意 有N个灯和M个开关,每个开关控制着一些灯,如果按下某个开关,就会让对应的灯切换状态:问在每个开关按下与否的一共2^m情况下,每种状态下亮灯的个数的立方的和. 思路1.首先注意到N< ...
- Codeforces 407B Long Path(好题 DP+思维)
题目链接:http://codeforces.com/problemset/problem/407/B 题目大意:一共n+1个房间,一个人从1走到n+1,每次经过房间都会留下一个标记,每个房间有两扇门 ...
随机推荐
- Redux Todos Example
此项目模板是使用Create React App构建的,它提供了一种简单的方法来启动React项目而无需构建配置. 使用Create-React-App构建的项目包括对ES6语法的支持,以及几种非官方 ...
- Kubernetes之ServiceAccount
ServiceAccount 是什么 Service Account为Pod中的进程和外部用户提供身份信息.所有的kubernetes集群中账户分为两类,Kubernetes管理的serviceacc ...
- kubernetes云平台管理实战: 高级资源deployment-滚动升级(八)
一.通过文件创建deployment 1.创建deployment文件 [root@k8s-master ~]# cat nginx_deploy.yml apiVersion: extensions ...
- DirectX11 With Windows SDK--11 混合状态与光栅化状态
前言 虽然这一部分的内容主要偏向于混合(Blending),但这里还需提及一下,关于渲染管线可以绑定的状态主要有如下四种: 光栅化状态(光栅化阶段) 采样器状态(像素着色阶段) 混合状态(输出合并阶段 ...
- 第六节:框架搭建之EF的Fluent Api模式的使用流程
一. 前言 沉寂了约一个月的时间,今天用一篇简单的文章重新回归博客,主要来探讨一下Fluent Api模式在实际项目中的使用流程. 1. Fluent API属于EF CodeFirst模式的一种,E ...
- JavaScript 正则表达式基础语法
前言 正则表达式在人们的印象中可能是一堆无法理解的字符,但就是这些符号却实现了字符串的高效操作.通常的情况是,问题本身并不复杂,但没有正则表达式就成了大问题.javascript中的正则表达式作为相当 ...
- mask rcnn
Mask RCNN Mask RCNN 中主要改进是在faster rcnn中box regression 的branch 上加入mask prediction branch,能够得到点到点的预测. ...
- 目标检测网络之 YOLOv3
本文逐步介绍YOLO v1~v3的设计历程. YOLOv1基本思想 YOLO将输入图像分成SxS个格子,若某个物体 Ground truth 的中心位置的坐标落入到某个格子,那么这个格子就负责检测出这 ...
- Linux Shell远程执行命令
1.问题描述 经常需要远程到其他节点上执行一些shell命令,如果分别ssh到每台主机上再去执行很麻烦,下边介绍shell命令远程执行的方法. 前提: 远程电脑之间已经配置ssh免密码登陆 2.脚本方 ...
- LeetCode第十二题-将数字转化为罗马数字
Integer to Roman 问题简介:将输入的int类型数字转化为罗马数字 问题详解:罗马数字由七个不同的符号表示:I,V,X,L,C,D和M 符号-数值 I - 1 V - 5 X -10 L ...