【Luogu】P2766最长不下降子序列问题(暴力网络流)
2024-09-30 00:54:09
水题qwq,数据都那么水。
我要是出数据的人我就卡$n^3$建图。
qwq。
然而这么水的题我!居!然!没!有!1!A!!还!提!交!了!五!遍!!!
md从现在开始要锻炼1A率了
看我从今往后做完一道题之后至少检查TM十分钟
可恶qwq。
第一问$n^2$sbDP可解。然而你们知道我提交五遍TM是错在哪里了吗?????
我TM就错在这个pj-,sb到不能再sb的sb暴力DP上!!!
气死我了!!!
关于第二问和第三问,先拆点再拆点qwq。
先把每个点拆成入点和出点用来限制流量,然后把出点拆成s个点为的是跑最大流的时候统计比较方便。
能跑到汇点就是一种方式,跑不到就说明没有那么长的最长不下降子序列。
第三问暴力修改边容量再跑一遍即可。
qwqqqqqqqqq我网络流都没写错栽在一个入门DP上
qwqqqqqq
#include<cstdio>
#include<cctype>
#include<algorithm>
#include<cstring>
#include<cstdlib>
#include<queue>
#define maxm 200010
using namespace std;
inline long long read(){
long long num=,f=;
char ch=getchar();
while(!isdigit(ch)){
if(ch=='-') f=-;
ch=getchar();
}
while(isdigit(ch)){
num=num*+ch-'';
ch=getchar();
}
return num*f;
} int n,m;
int q[];
int s[];
int ne[][],tot[];
int Start,End;
inline int count(int i){ return i&?i+:i-; } struct Edge{
int next,to,val;
}edge[maxm*];
int head[maxm],num;
inline void addedge(int from,int to,int val){
edge[++num]=(Edge){head[from],to,val};
head[from]=num;
}
inline void add(int from,int to,int val){
addedge(from,to,val);
addedge(to,from,);
} bool vis[maxm];
int dfn[maxm];
int list[maxm]; bool bfs(){
memset(vis,,sizeof(vis)); vis[Start]=; dfn[Start]=;
queue<int>q; q.push(Start);
while(!q.empty()){
int from=q.front();q.pop();
for(int i=head[from];i;i=edge[i].next){
int to=edge[i].to;
if(edge[i].val==||vis[to]) continue;
vis[to]=; dfn[to]=dfn[from]+;
q.push(to);
}
}
return vis[End];
} int dfs(int x,int val){
if(x==End||val==) return val;
int flow=; vis[x]=;
for(int &i=list[x];i;i=edge[i].next){
int to=edge[i].to;
if(vis[to]||edge[i].val==||dfn[to]!=dfn[x]+) continue;
int now=dfs(to,min(val,edge[i].val));
val-=now; flow+=now; edge[i].val-=now; edge[count(i)].val+=now;
if(val<=) break;
}
if(val!=flow) dfn[x]=-;
return flow;
} int maxflow(){
int ans=;
while(bfs()){
memset(vis,,sizeof(vis));
for(int i=Start;i<=End;++i) list[i]=head[i];
int now=dfs(Start,0x7fffffff);
if(!now) break;
ans+=now;
}
return ans;
} int main(){
int n=read();
for(int i=;i<=n;++i){
q[i]=read();s[i]=;
}
m=;
for(int i=n-;i>=;--i){
for(int j=i+;j<=n;++j){
if(q[i]>q[j]) continue;
if(s[i]<s[j]+) s[i]=s[j]+;
if(m<s[i]) m=s[i];
}
}
printf("%d\n",m);
End=(m+)*n+n+;
for(int i=;i<=n;++i){
add(Start,i,0x7fffffff);
add(i,i+n,);
add(m*n+i,End,);
for(int j=i+;j<=n;++j){
if(q[i]>q[j]) continue;
for(register int k=;k<m;++k) add(k*n+i,(k+)*n+j,);
}
}
int ans=maxflow();
printf("%d\n",ans);
int ret=;
for(int j=head[ret];j;j=edge[j].next){
int to=edge[j].to;
if(to<ret) continue;
edge[j].val+=;
}
/*for(int j=head[1];j;j=edge[j].next){
int to=edge[j].to;
edge[j].val+=1234567;
}*/
ret=m*n+n;
for(int j=head[ret];j;j=edge[j].next){
int to=edge[j].to;
if(to!=End) continue;
edge[j].val+=;
}
printf("%d",ans+maxflow());
return ;
}
最新文章
- 7.让网站支持http和https的访问方式
- Laxcus大数据管理系统2.0(13)- 总结
- webservice之XFire的使用(java调用java)
- Log.i()的用法
- const和#define的区别
- 强大的CImage类
- Unity3D 之UGUI 滑动条(Slider)
- [C#]『Barrier』任务并行库使用小计
- [Java] LinkedList / Queue - 源代码学习笔记
- HDU 1757 A Simple Math Problem(矩阵高速幂)
- C/C++ 结构体成员在内存中的对齐规则
- Java经典编程题50道之五十
- There is no action xxxFun defined for api controller api/subitem
- OO第三次博客作业---透过代码看设计
- BZOJ2190 [SDOI2008]仪仗队 [欧拉函数]
- CentOS7使用打开关闭防火墙与端口
- bzoj 1288: Neighbours
- 《Two Dozen Short Lessons in Haskell》(二十二)递归
- 安装apache+php记录
- 如果天空不死博客java阅读列表整理
热门文章
- Codeforces Round #321 (Div. 2) D Kefa and Dishes(dp)
- CDOJ 485 UESTC 485 	Game (八数码变形,映射,逆cantor展开)
- 2018.4.22 深入理解Java的接口和抽象类
- 【转】Java8学习笔记(1) -- 从函数式接口说起
- mysql中的空值问题
- linux - mysql 安装教程
- 为什么方差的分母有时是n,有时是n-1 源于总体方差和样本方差的不同
- 【dsu || 线段树合并】bzoj4756: [Usaco2017 Jan]Promotion Counting
- intellij idea 下载安装破解教程
- (70)zabbix telnet监控类型