「SDOI2014」Lis 解题报告
「SDOI2014」Lis
题目描述
给定序列 \(A\),序列中的每一项 \(A_i\) 有删除代价 \(B_i\) 和附加属性 \(C_i\)。
请删除若干项,使得 \(A\) 的最长上升子序列长度减少至少 \(1\),且付出的代价之和最小,并输出方案。
如果有多种方案,请输出将删去项的附加属性排序之后,字典序最小的一种。
\(T\le 5,1\le n\le 700,1\le A_i,B_i,C_i\le 10^9\)
上午想这个题的时候,有一种迷之直觉,让我去做网络流24题的最长不下降子序列那题,然后我就去了,然后建模果然是一样的,太神奇啦
按\(dp\)值把点分层,然后把入度为\(0\)的和\(s\)连,把出度为\(0\)的和\(t\)连,发现问题转化成割掉权值和最小的点,使\(s,t\)不连通,可以用拆点最小割做。
然后考虑如何做字典序最小,我们可以从小到大枚举\(C\),看这条边可不可以成为最小割上的边。
一条边可以是割集上的边\((u,v)\)得在残余网络上满足
- 已经流满了
- \(u\)不可以通过残余网络到\(v\)
可以从\(s\)遍历一遍看看\(u,v\)是否都被遍历或者都没被遍历。
然而这个题在枚举这个边合法后得钦定它,所以得删了它重新建图重新跑,但这肯定不能搞。
考虑一个退流操作,如果在残余网络上删去边\((u,v)\),可以先跑两边最大流\((t,v)\)和\((u,s)\)的,然后把这条边流量清空,这样复杂度就是对的了(
Code:
#include <cstdio>
#include <cctype>
#include <cstring>
#include <algorithm>
using std::min;
using std::max;
int read()
{
int x=0;char c=getchar();
while(!isdigit(c)) c=getchar();
while(isdigit(c)) x=x*10+c-'0',c=getchar();
return x;
}
const int N=1500;
const int M=3e5;
const int inf=0x3f3f3f3f;
int head[N],to[M],Next[M],edge[M],cnt;
void add(int u,int v,int w)
{
to[++cnt]=v,edge[cnt]=w,Next[cnt]=head[u],head[u]=cnt;
to[++cnt]=u,edge[cnt]=0,Next[cnt]=head[v],head[v]=cnt;
}
int q[N],dep[N],l,r;
bool bfs(int s,int t)
{
memset(dep,0,sizeof dep);
dep[q[l=r=1]=s]=1;
while(l<=r)
{
int now=q[l++];
for(int v,i=head[now];i;i=Next[i])
if(edge[i]&&!dep[v=to[i]])
{
dep[v]=dep[now]+1;
if((q[++r]=v)==t) return true;
}
}
return false;
}
int dfs(int now,int flow,int t)
{
if(now==t) return flow;
int res=flow,yuu;
for(int v,i=head[now];i&&res;i=Next[i])
if(edge[i]&&dep[v=to[i]]==dep[now]+1)
{
yuu=dfs(v,min(res,edge[i]),t);
if(!yuu){dep[v]=0;continue;}
edge[i]-=yuu,edge[i^1]+=yuu;
res-=yuu;
}
return flow-res;
}
int Dinic(int s,int t)
{
int ret=0,flow;
while(bfs(s,t)) while(flow=dfs(s,inf,t)) ret+=flow;
return ret;
}
int dp[N],sta[N],tot,a[N];
struct node
{
int id,x;
bool friend operator <(node a,node b){return a.x<b.x;}
}yuu[N];
void work()
{
cnt=1;memset(head,0,sizeof head);
int n=read(),ans=0;
for(int i=1;i<=n;i++)
{
a[i]=read();
yuu[i].id=i;
dp[i]=1;
for(int j=1;j<i;j++)
if(a[j]<a[i])
dp[i]=max(dp[i],dp[j]+1);
ans=max(ans,dp[i]);
}
for(int i=1;i<=n;i++)
add(i,i+n,read());
int s=n<<1|1,t=s+1;
for(int i=1;i<=n;i++)
{
yuu[i].x=read();
if(dp[i]==1) add(s,i,inf);
if(dp[i]==ans) add(i+n,t,inf);
for(int j=1;j<i;j++)
if(a[j]<a[i]&&dp[j]+1==dp[i])
add(j+n,i,inf);
}
printf("%d ",Dinic(s,t));
std::sort(yuu+1,yuu+1+n);
for(int i=1;i<=n;i++)
{
int u=yuu[i].id;
if(edge[u<<1]||bfs(u,u+n)) continue;
sta[++tot]=u;
Dinic(u,s),Dinic(t,u+n);
edge[u<<1]=edge[u<<1|1]=0;
}
std::sort(sta+1,sta+1+tot);
printf("%d\n",tot);
for(int i=1;i<=tot;i++) printf("%d ",sta[i]);
}
int main()
{
int T=read();
while(T--) work();
return 0;
}
2019.2.21
「SDOI2014」Lis 解题报告的更多相关文章
- 「SDOI2014」重建 解题报告
「SDOI2014」重建 题意 给一个图\(G\),两点\((u,v)\)有边的概率是\(p_{u,v}\),求有\(n-1\)条边通行且组成了一颗树的概率是多少. 抄了几个矩阵树定理有趣的感性说法 ...
- 「ZJOI2016」旅行者 解题报告
「ZJOI2016」旅行者 对网格图进行分治. 每次从中间选一列,然后枚举每个这一列的格子作为起点跑最短路,进入子矩形时把询问划分一下,有点类似整体二分 至于复杂度么,我不会阿 Code: #incl ...
- 「HNOI2016」树 解题报告
「HNOI2016」树 事毒瘤题... 我一开始以为每次把大树的子树再接给大树,然后死活不知道咋做,心想怕不是个神仙题哦 然后看题解后才发现是把模板树的子树给大树,虽然思维上难度没啥了,但是还是很难写 ...
- 「HNOI2016」序列 解题报告
「HNOI2016」序列 有一些高妙的做法,懒得看 考虑莫队,考虑莫队咋移动区间 然后你在区间内部找一个最小值的位置,假设现在从右边加 最小值左边区间显然可以\(O(1)\),最小值右边的区间是断掉的 ...
- 「HNOI2016」网络 解题报告
「HNOI2016」网络 我有一个绝妙的可持久化树套树思路,可惜的是,它的空间是\(n\log^2 n\)的... 注意到对一个询问,我们可以二分答案 然后统计经过这个点大于当前答案的路径条数,如果这 ...
- 「HAOI2018」染色 解题报告
「HAOI2018」染色 是个套路题.. 考虑容斥 则恰好为\(k\)个颜色恰好为\(c\)次的贡献为 \[ \binom{m}{k}\sum_{i\ge k}(-1)^{i-k}\binom{m-k ...
- 「HNOI2016」最小公倍数 解题报告
「HNOI2016」最小公倍数 考虑暴力,对每个询问,处理出\(\le a,\le b\)的与询问点在一起的联通块,然后判断是否是一个联通块,且联通块\(a,b\)最大值是否满足要求. 然后很显然需要 ...
- 「SCOI2016」围棋 解题报告
「SCOI2016」围棋 打CF后困不拉基的,搞了一上午... 考虑直接状压棋子,然后发现会t 考虑我们需要上一行的状态本质上是某个位置为末尾是否可以匹配第一行的串 于是状态可以\(2^m\)压住了, ...
- 「SCOI2016」妖怪 解题报告
「SCOI2016」妖怪 玄妙...盲猜一个结论,然后过了,事后一证,然后假了,数据真水 首先要最小化 \[ \max_{i=1}^n (1+k)x_i+(1+\frac{1}{k})y_i \] \ ...
随机推荐
- 百度地图开发者API学习笔记二
一,地图上多个覆盖物(Marker). 当有多个覆盖物时,我们需要获取每个点的信息.如下图,每个Marker的经度都不相同 二,代码: <!DOCTYPE html> <html&g ...
- Linux 典型应用之缓存服务
memcached 安装和简单使用 yum install memcached 启动 -d 表示以守护进程的方式启动 memcached -d 安装telnet 它可以检测某个端口是否是通的,可以发送 ...
- CentOS云厂商清单
Download CentOShttps://www.centos.org/download/ Download - CentOS Wikihttps://wiki.centos.org/Downlo ...
- js 翻牌活动效果
直接上代码 html: <div class="index_main"> <ul class="index_card"> <li ...
- Vue+iview实现添加删除类
<style> .tab-warp{ border-bottom: solid 1px #e0e0e0; overflow: hidden; margin-top: 30px; posit ...
- C# Note31: 如何使用Visual Studio做单元测试
待更! 使用Visual Studio 2013进行单元测试--初级篇 带你玩转Visual Studio——单元测试(C++例)
- Hbase表结构模型
- 1、通过eureka创建注册中心
第一个demo(用户需要调用电影服务) 1.创建项目 new starter project 勾选上Eureka Server 2.编写application.yml #配置端口 server: po ...
- docker 列出每个容器的IP
抄来的...找不到出处了. 常用方法有两种 docker inspect 容器ID | grep IPAddress 方法二 查看docker name: sudo docker inspect ...
- 使用javaWeb的二大(Listener、Filter)组件实现分IP统计访问次数
分析: 统计工作需要在所有资源之前都执行,那么就可以放到Filter中. 我们这个过滤器不打算做拦截操作!因为我们只是用来做统计 用什么东西来装载统计的数据.Map<String,Integer ...