题目大意:
现在让我们求一个子序列(逆序数/子序列长度)的最大值。
思路:
经典的最大密度子图问题
①首先用到01分数规划。
因为要求X/Y的最大值,那么我们可以定义:F(Mid)=X-Mid*Y
如果有F(MID)>=0,那么就有:X-mid*Y>=0,那么就有X/Y>=Mid
那么观察到问题具有单调性,那么Mid的值,我们二分处理即可。
②接下来我们就要Check 这个Mid值。
对于最大密度子图问题,我们将边设定为左集合,将点设定为右集合,做一个最大权闭包模型即可。
对于这个问题,我们将能够构成逆序对的边构成左集合,将每个点设定为右集合,然后建图方式如下:
①建立源点,连入各个左集合的边,流为1.
②建立汇点,将各个右集合的边连入汇点,流为Mid.
③将构成逆序对的边连入两个相关联的点上,流为INF.
然后跑出最大流==最小割,那么问题的利益点就收获了:Tot-Maxflow,这里Tot表示构成的逆序对的边数、
过程维护一下,注意精度就行了。
另外,数组别开太大,开太大喜欢用memset的同学,会TLE= =
Ac代码:
#include<stdio.h>
#include<string.h>
#include<iostream>
#include<queue>
using namespace std;
double eps=1e-7;
struct node
{
int from,to,next;
double w;
} e[1500000];
int head[15050];
int cur[15050];
int divv[15050];
int a[15050];
int n,ss,tt,cont;
void add(int from,int to,double w)
{
e[cont].to=to;
e[cont].w=w;
e[cont].next=head[from];
head[from]=cont++;
}
int makedivv()
{
memset(divv,0,sizeof(divv));
divv[ss]=1;
queue<int >s;
s.push(ss);
while(!s.empty())
{
int u=s.front();
if(u==tt)return 1;
s.pop();
for(int i=head[u];i!=-1;i=e[i].next)
{
double w=e[i].w;
int v=e[i].to;
if(divv[v]==0&&w>eps)
{
divv[v]=divv[u]+1;
s.push(v);
}
}
}
return 0;
}
double Dfs(int u,double maxflow,int tt)
{
if(u==tt)return maxflow;
double ret=0;
for(int &i=cur[u];i!=-1;i=e[i].next)
{
int v=e[i].to;
double w=e[i].w;
if(divv[v]==divv[u]+1&&w>eps)
{
double f=Dfs(v,min(maxflow-ret,w),tt);
e[i].w-=f;
e[i^1].w+=f;
ret+=f;
if(ret==maxflow)return ret;
}
}
return ret;
}
double Dinic()
{
double ans=0;
while(makedivv()==1)
{
memcpy(cur,head,sizeof(head));
ans+=Dfs(ss,(double)100000000,tt);
}
return ans;
}
int Check(double mid)
{
ss=0;
cont=0;
int cnt=0;
memset(head,-1,sizeof(head));
for(int i=1; i<=n; i++)
{
for(int j=1; j<i; j++)
{
if(a[j]>a[i])
{
++cnt;
add(ss,n+cnt,1);
add(n+cnt,ss,0);
add(n+cnt,i,(double)10000000);
add(i,n+cnt,0);
add(n+cnt,j,(double)10000000);
add(j,n+cnt,0);
}
}
}
tt=cnt+n+1;
for(int i=1;i<=n;i++)
{
add(i,tt,mid);
add(tt,i,0);
}
double ans=Dinic();
double temp=cnt;
if(temp-ans<=0)return 0;
else return 1;
}
int main()
{
int kase=0;
int t;
scanf("%d",&t);
while(t--)
{
scanf("%d",&n);
for(int i=1; i<=n; i++)scanf("%d",&a[i]);
double l=0;
double r=n*n;
double ans=0;
while(r-l>eps)
{
double mid=(l+r)/2;
if(Check(mid)==1)
{
ans=max(ans,mid);
l=mid;
}
else r=mid;
}
printf("Case #%d: ",++kase);
printf("%.7f\n",ans);
}
}