2014西安现场赛C(最大密度子图->分数规划+最大权封闭子图)

最新推荐文章于 2021-11-04 15:32:01 发布

qkoqhh

最新推荐文章于 2021-11-04 15:32:01 发布

阅读量294

点赞数

分类专栏：网络流分数规划

本文链接：https://blog.csdn.net/qkoqhh/article/details/82117500

版权

网络流同时被 2 个专栏收录

34 篇文章 0 订阅

订阅专栏

分数规划

4 篇文章 0 订阅

订阅专栏

题目链接：http://codeforces.com/gym/100548

这个题还是比较好的。。首先把能形成逆序对的点连边，发现题目要求的，其实是一个最大的子图密度。。

然后很容易联想到封闭子图，然而密度是用除法表示的。。

因此需要分数规划一下，令二分答案为g，要求的是ΣV-gΣE=0的解。。(V代表边集，E代表点集)

这样可以构造出封闭子图了，把每个逆序对代表的边看成点u，将源点与每个u连边，容量为1，把原图的点v与汇点连边，容量为g*1，然后u向其两个顶点连边，容量为inf

虽然有n^2个点还带个log然而跑起来挺快的。。而且时限也有7s。。

另一种做法来自胡伯涛的论文《最小割模型在信息学竞赛中的应用》，这里将这个方法复述一遍。。

其思路是将ΣV-gΣE做如下变换：

每枚举一个g，我们需要找到ΣV-gΣE的最大值来作为check的依据，即找gΣE-ΣV的最小值

考虑与E中的一点u相连的所有的边，共有 $d_u$ 条(设 $d_u$ 为u的度数)，然后再去掉V中不包含的边，而V其实满足其顶点都在E内，设V'为V的补集，将 $d_u$ 减去V与V'的割集(设为c[V,V'])后除以2就是V的边数了(因为每条边都被其顶点算了一次)，即

$g\sum V-\sum E=g\sum V-(\sum_{u\in V}d_u-c[V,V'])/2$

同乘以2得

$\sum_{u \in V} (2g-d_u)+c[V,V']$

接下来利用最小割求解时为了避免出现容量为负的情况，我们手动加个大正数U下去，让S给每个点连边，容量为U，每个点与T连边，容量为U+2g- $d_u$ ，然后原图的无向边直接拆成2个容量为1的边，跑最小割c[S,T]，而事实上T集合正是我们要求的V，于是有

$c[S,T]=\sum_{u \in S}U+\sum_{v \in T}(U+2g-d_u)+\sum_{u \in s ,v\in T}1=Un+\sum_{v\in T}(2g-d_u)+c[V,V']$

所以将c[S,T]减去Un就能求得最小值，而相应地，ΣV-gΣE的最大值为(Un-c[S,T])/2

然后问题就解决了。。。

这样相比起解法一，这个网络流的点只有n个，可以放心地跑。。

然而事实上跑起来只比解法一快个700 ms？？(cf评测机是1700ms)

解法一：

/**
 *　　　　　　  　　┏┓　　 　┏┓
 * 　　　　　  　　┏┛┗━━━━━━━┛┗━━━┓
 * 　　　　　  　　┃　　　　　　　┃ 　
 * 　　　　　  　　┃　　　━　　 　┃
 * 　　　　　  　　┃　＞　　　＜　┃
 * 　　　　　  　　┃　　　　　　　┃
 * 　　　　　  　　┃...　⌒　... 　┃
 * 　　　　　　  　┃              ┃
 * 　　　　　　  　┗━┓          ┏━┛
 * 　　　　　　　　　┃          ┃　Code is far away from bug with the animal protecting　　　　　　　　　　
 * 　　　　　　　　　┃          ┃   神兽保佑,代码无bug
 * 　　　　　　　　　┃          ┃　　　　　　　　　　　
 * 　　　　　　　　　┃          ┃  　　　　　　
 * 　　　　　　　　　┃          ┃
 * 　　　　　　　　　┃          ┃　　　　　　　　　　　
 * 　　　　　　　　　┃          ┗━━━┓
 * 　　　　　　　　　┃              ┣┓
 * 　　　　　　　　　┃              ┏┛
 * 　　　　　　　　　┗┓┓┏━━━━━━━━┳┓┏┛
 * 　　　　　　　　　　┃┫┫       ┃┫┫
 * 　　　　　　　　　　┗┻┛       ┗┻┛
 */
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<iostream>
#include<queue>
#include<cmath>
#include<map>
#include<stack>
#include<set>
#include<bitset>
#define inc(i,l,r) for(int i=l;i<=r;i++)
#define dec(i,l,r) for(int i=l;i>=r;i--)
#define link(x) for(edge *j=h[x];j;j=j->next)
#define mem(a) memset(a,0,sizeof(a))
#define ll long long
#define eps 1e-10
#define succ(x) (1LL<<(x))
#define lowbit(x) (x&(-x))
#define sqr(x) ((x)*(x))
#define mid ((x+y)/2.0)
#define NM 6005
#define nm 30005
#define N 1000005
#define M(x,y) x=max(x,y)
const double pi=acos(-1);
const int inf=1e9+7;
using namespace std;
ll read(){
    ll x=0,f=1;char ch=getchar();
    while(!isdigit(ch)){if(ch=='-')f=-1;ch=getchar();}
    while(isdigit(ch))x=x*10+ch-'0',ch=getchar();
    return f*x;
}
  




struct edge{int t;double v;edge*next,*rev;}e[nm],*h[NM],*o=e,*p[NM],*tmp[NM];
void _add(int x,int y,double v){o->t=y;o->v=v;o->next=h[x];h[x]=o++;}
void add(int x,int y,double v){_add(x,y,v);_add(y,x,0);h[x]->rev=h[y];h[y]->rev=h[x];}
struct tmp{int x,y;}b[NM];
int n,a[NM],m,tot,ca;
int d[NM],cnt[NM];
double ans;

double maxflow(){
    double flow=0,s=inf;edge*j;
    mem(cnt);mem(d);mem(p);mem(tmp);
    inc(i,0,n)tmp[i]=h[i];
    cnt[0]=n+1;
    for(int x=0;d[x]<n+1;){
	for(j=tmp[x];j;j=j->next)if(j->v>eps&&d[j->t]+1==d[x])break;
	if(j){
	    s=min(s,j->v);tmp[x]=p[j->t]=j;
	    if((x=j->t)==n){
		for(;x;x=p[x]->rev->t)p[x]->v-=s,p[x]->rev->v+=s;
		flow+=s;s=inf;
	    }
	}else{
	    if(!--cnt[d[x]])break;d[x]=n+1;
	    link(x)if(j->v>eps&&d[x]>d[j->t]+1)tmp[x]=j,d[x]=d[j->t]+1;
	    cnt[d[x]]++;
	    if(p[x])x=p[x]->rev->t;
	}
    }
    return flow;
}


bool check(double t){
    double s=tot;
    mem(e);mem(h);o=e;
    inc(i,1,tot)add(0,i,1);
    inc(i,1,m)add(i+tot,n,t);
    inc(i,1,tot)add(i,b[i].x+tot,inf),add(i,b[i].y+tot,inf);
    s-=maxflow();
    return s>eps;
}

int main(){
    int _=read();while(_--){
	m=read();tot=0;ans=0;
	inc(i,1,m)a[i]=read();
	inc(i,1,m)inc(j,i+1,m)if(a[i]>a[j])b[++tot].x=i,b[tot].y=j;
	n=tot+1+m;
	for(double x=0,y=1000;x+eps<y;)
	    if(check(mid))ans=mid,x=mid;else y=mid;
	printf("Case #%d: %.10lf\n",++ca,ans);
    }
    return 0;
}

解法二：

/**
 *　　　　　　  　　┏┓　　 　┏┓
 * 　　　　　  　　┏┛┗━━━━━━━┛┗━━━┓
 * 　　　　　  　　┃　　　　　　　┃ 　
 * 　　　　　  　　┃　　　━　　 　┃
 * 　　　　　  　　┃　＞　　　＜　┃
 * 　　　　　  　　┃　　　　　　　┃
 * 　　　　　  　　┃...　⌒　... 　┃
 * 　　　　　　  　┃              ┃
 * 　　　　　　  　┗━┓          ┏━┛
 * 　　　　　　　　　┃          ┃　Code is far away from bug with the animal protecting　　　　　　　　　　
 * 　　　　　　　　　┃          ┃   神兽保佑,代码无bug
 * 　　　　　　　　　┃          ┃　　　　　　　　　　　
 * 　　　　　　　　　┃          ┃  　　　　　　
 * 　　　　　　　　　┃          ┃
 * 　　　　　　　　　┃          ┃　　　　　　　　　　　
 * 　　　　　　　　　┃          ┗━━━┓
 * 　　　　　　　　　┃              ┣┓
 * 　　　　　　　　　┃              ┏┛
 * 　　　　　　　　　┗┓┓┏━━━━━━━━┳┓┏┛
 * 　　　　　　　　　　┃┫┫       ┃┫┫
 * 　　　　　　　　　　┗┻┛       ┗┻┛
 */
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<iostream>
#include<queue>
#include<cmath>
#include<map>
#include<stack>
#include<set>
#include<bitset>
#define inc(i,l,r) for(int i=l;i<=r;i++)
#define dec(i,l,r) for(int i=l;i>=r;i--)
#define link(x) for(edge *j=h[x];j;j=j->next)
#define mem(a) memset(a,0,sizeof(a))
#define ll long long
#define eps 1e-8
#define succ(x) (1LL<<(x))
#define lowbit(x) (x&(-x))
#define sqr(x) ((x)*(x))
#define mid ((x+y)/2.0)
#define NM 6005
#define nm 30005
#define N 1000005
#define M(x,y) x=max(x,y)
const double pi=acos(-1);
const int inf=1e9+7;
using namespace std;
ll read(){
    ll x=0,f=1;char ch=getchar();
    while(!isdigit(ch)){if(ch=='-')f=-1;ch=getchar();}
    while(isdigit(ch))x=x*10+ch-'0',ch=getchar();
    return f*x;
}
  




struct edge{int t;double v;edge*next,*rev;}e[nm],*h[NM],*o=e,*p[NM],*tmp[NM];
void _add(int x,int y,double v){o->t=y;o->v=v;o->next=h[x];h[x]=o++;}
void add(int x,int y,double v){_add(x,y,v);_add(y,x,0);h[x]->rev=h[y];h[y]->rev=h[x];}
struct tmp{int x,y;}b[NM];
int n,a[NM],m,tot,ca,c[NM];
int d[NM],cnt[NM];
double ans;

double maxflow(){
    double flow=0,s=inf;edge*j;
    mem(cnt);mem(d);mem(p);mem(tmp);
    inc(i,0,n)tmp[i]=h[i];
    cnt[0]=n+1;
    for(int x=0;d[x]<n+1;){
	for(j=tmp[x];j;j=j->next)if(j->v>eps&&d[j->t]+1==d[x])break;
	if(j){
	    s=min(s,j->v);tmp[x]=p[j->t]=j;
	    if((x=j->t)==n){
		for(;x;x=p[x]->rev->t)p[x]->v-=s,p[x]->rev->v+=s;
		flow+=s;s=inf;
	    }
	}else{
	    if(!--cnt[d[x]])break;d[x]=n+1;
	    link(x)if(j->v>eps&&d[x]>d[j->t]+1)tmp[x]=j,d[x]=d[j->t]+1;
	    cnt[d[x]]++;
	    if(p[x])x=p[x]->rev->t;
	}
    }
    return flow;
}


bool check(double t){
    double s;
    mem(e);mem(h);o=e;
    inc(i,1,m)add(0,i,m),add(i,n,m+2*t-c[i]);
    inc(i,1,tot)add(b[i].x,b[i].y,1),add(b[i].y,b[i].x,1);
    s=maxflow();
    return (sqr(m)-s)/2>eps;
}

int main(){
    int _=read();while(_--){
	m=read();tot=0;ans=0;mem(c);
	inc(i,1,m)a[i]=read();
	inc(i,1,m)inc(j,i+1,m)if(a[i]>a[j])b[++tot].x=i,b[tot].y=j,c[i]++,c[j]++;
	n=1+m;
	for(double x=0,y=100;x+eps<y;)
	    if(check(mid))ans=mid,x=mid;else y=mid;
	printf("Case #%d: %.10lf\n",++ca,ans);
    }
    return 0;
}

qkoqhh

关注

0
点赞
踩
0

收藏

觉得还不错? 一键收藏
0
评论
2014西安现场赛C(最大密度子图->分数规划+最大权封闭子图)

题目链接：http://codeforces.com/gym/100548这个题还是比较好的。。首先把能形成逆序对的点连边，发现题目要求的，其实是一个最大的子图密度。。然后很容易联想到封闭子图，然而密度是用除法表示的。。因此需要分数规划一下，令二分答案为g，要求的是ΣV-gΣE=0的解。。(V代表边集，E代表点集)这样可以构造出封闭子图了，把每个逆序对代表的边看成点u，将源点与每个...
复制链接

扫一扫