P3705 [SDOI2017]新生舞会(01分数规划+费用流)

P3705 [SDOI2017]新生舞会

见识到了大名鼎鼎的$01$分数规划…

$显然这是个二分图,但是这个价值计算放置太头疼了....$

$令\frac{\sum a_i}{\sum b_i}=maxx$

稍微变化一下$maxx\ast \sum b_i=\sum a_i$

也就是$\sum (a_i-maxx\ast b_i)=0$

$那我们二分maxx,把边权看成a_i-maxx\ast b_i$

$求一个最大费用最大流,随着maxx的增大费用肯定会变小!!$

$所以二分maxx,若此时费用仍然大于0,说明maxx可以继续增大.$

很神奇吧?不复杂,却能巧妙地解决问题呢

#include 
using namespace std;
#define int long long
const int maxn=2e5+10;
const double eps=1e-7;
const int inf=1e9;
int n,m,s,t;
double mincost;
int head[maxn<<1],cnt=1,incf[maxn],pre[maxn],vis[maxn];
int a[109][109],b[109][109];
struct edge{
	int to,nxt,flow;double w;//分别代表 
}d[maxn<<1];
void add(int u,int v,int flow,double w)//最大流量,单位费用
{
	d[++cnt]=(edge){v,head[u],flow,w},head[u]=cnt;
	d[++cnt]=(edge){u,head[v],0,-w},head[v]=cnt;
} 
double dis[maxn];
bool spfa()
{
	queueq;
	for(int i=0;i<=t;i++)	dis[i]=-inf;
	memset(vis,0,sizeof(vis));
	q.push(s);
	dis[s]=0,vis[s]=1;
	incf[s] = inf;//初始流量无限大
	while( !q.empty() )
	{
		int u=q.front(); q.pop();
		vis[u]=0;//出队
		for(int i=head[u];i;i=d[i].nxt)
		{
			if( !d[i].flow )	continue;//无流量了	
			int v=d[i].to;
			if( dis[v]=0 )	return true;
	return false;
}
signed main()
{
	cin >> n;
	for(int i=1;i<=n;i++)
	for(int j=1;j<=n;j++)
		cin >> a[i][j];
	for(int i=1;i<=n;i++)
	for(int j=1;j<=n;j++)
		cin >> b[i][j];
	double l=0,r=1e6,mid,ans;
	while( (r-l)>=eps )
	{
		mid=(l+r)/2.0;
		if( isok(mid) )	l=mid,ans=mid;
		else	r=mid;
	}
	printf("%.6f",ans);
}