您的位置:首页 > 其它

luogu P2766 最长不下降子序列问题

2018-03-08 14:04 274 查看
关于构图:
首先,第一问dp(O(n^2) or O(nlogn)这里就不用单调队列优化了)稳稳的。
很明显用网络流做,明明是看到标签。
很明显第三问是第二问的升级版,只需将超级源点和汇点到相应点的边权改为INF即可。
构图(p[i]表示原序列第i个的数,dp[i]表示以i为终点的最长不下降子序列的长度):
1.超级源点和汇点:st=0,ed=n*2+1。
2.从i'向i建一条流量为1的边,这样它就可以流回去,进行下一次选择(另一个最长不下降子序列)。
3.判断,当dp[i]为1时,表示它是某些(个)最长不下降子序列的起点,由源点向它建一条流量为1的边。
4.判断,当dp[i]为maxx时,表示它是某些(个)最长不下降子序列的终点,由它向终点建一条流量为1的边。
5.判断,当p[i]<=p[j]&&(dp[j]==dp[i]+1||ans==1)时,表示它是某些(个)最长不下降子序列的中间点(不为起点和终点),由i向j'建一条流量为1的边。

注意:第二问和第三问分别求最大流。

关于代码:#include<cstdio>
#include<cstring>
#include<queue>
#include<algorithm>
#define INF 2147483647
using namespace std;
queue<int> f;
struct node{int x,y,z,next;} a[100000];
int last[10000],p[1000],dp[1000];
int n,m,len,st,ed,ans=0;
void ins(int x,int y,int z)
{
a[++len].x=x;a[len].y=y;a[len].z=z;a[len].next=last[x];last[x]=len;
}
void init(int x)
{
len=-1;
memset(last,-1,sizeof(last));
for(int i=1;i<=n;i++)
{
ins(i+n,i,1),ins(i+n,i,0);
if(dp[i]==1) ins(st,i,i==1?x:1),ins(i,st,0);
if(dp[i]==ans) ins(i,ed,i==n?x:1),ins(ed,i,0);
for(int j=i+1;j<=n;j++)
if(p[i]<=p[j]&&(dp[j]==dp[i]+1||ans==1)) ins(i,n+j,1),ins(n+j,i,0);
}
}
int h[100000];
bool bfs()
{
memset(h,0,sizeof(h));
h[st]=1;
f.push(st);
while(!f.empty())
{
int x=f.front();
for(int i=last[x];i>=0;i=a[i].next)
{
int y=a[i].y;
if(a[i].z>0&&h[y]==0)
{
h[y]=h[x]+1;
f.push(y);
}
}
f.pop();
}
if(h[ed]) return true; else return false;
}
int dfs(int x,int f)
{
int s=0,t;
if(x==ed) return f;
for(int i=last[x];i>=0;i=a[i].next)
{
int y=a[i].y;
if(a[i].z>0&&h[y]==h[x]+1&&f>s)
{
s+=(t=(dfs(y,min(f-s,a[i].z))));
a[i].z-=t;
a[i^1].z+=t;
}
}
if(!s) h[x]=0;
return s;
}
int dinic()
{
int sum=0;
while(bfs())
sum+=dfs(st,INF);
return sum;
}
int main()
{
scanf("%d",&n);
st=0,ed=2*n+1;
for(int i=1;i<=n;i++)
{
scanf("%d",&p[i]);
dp[i]=1;
for(int j=1;j<i;j++)
if(p[j]<=p[i]) dp[i]=max(dp[i],dp[j]+1);
ans=max(ans,dp[i]);
}
printf("%d\n",ans);
init(1);
printf("%d\n",dinic());
init(INF);
printf("%d\n",dinic());
}
内容来自用户分享和网络整理,不保证内容的准确性,如有侵权内容,可联系管理员处理 点击这里给我发消息
标签: