省选专练[SDOI2017]新生舞会

我发现只要是带舞会二字的必然和网络流有关，原因是舞会会出现匹配问题。

然后证明了，当前各省省选的DAY2T1的难度大体就是分数规划套一个省选知识点的check函数，HAOI，SCOI，SDOI都有过。

然后这是标准分数规划模型

sigmaAi/sigmaBi=C

乘过来于是Ai-C*Bi=0；

当C较小的时候便大于0；

于是用费用流或者KMcheck一下就好了

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<queue>
using namespace std;
#define INF 1e18+120
inline void read(int &x){
	x=0;
	int f=1;
	char ch=getchar();
	while(ch<'0'||ch>'9'){
		if(ch=='-')
			f=-1;
		ch=getchar();
	}
	while(ch>='0'&&ch<='9'){
		x=x*10+ch-'0';
		ch=getchar();
	}
	x*=f;
}
struct Front_star{
	int u,v,w,nxt;
	double c;
}e[40000*4+120];
double a[200][200]={0};
double b[200][200]={0};
int n;
int first[500]={0};
int cnt=1;
void addedge(int u,int v,int w,double c){
	cnt++;
	e[cnt].u=u;
	e[cnt].v=v;
	e[cnt].w=w;
	e[cnt].c=c;
	e[cnt].nxt=first[u];
	first[u]=cnt;
}
void add(int u,int v,int w,double c){
	addedge(u,v,w,c);
	addedge(v,u,0,-c);
}
int S=0;
int T=251;
void build(double ans){
	cnt=1;
	memset(first,0,sizeof(first));
	for(int i=1;i<=n;i++){
		for(int j=1;j<=n;j++){
			add(i,j+n,1,a[i][j]-ans*b[i][j]);
		}
	}
	for(int i=1;i<=n;i++){
		add(S,i,1,0);
		add(i+n,T,1,0);
	}
}
queue<int> q;
double dis[500]={0};
int inqueue[500]={0};
int pre[500]={0};
bool SPFA(){
	for(int i=S;i<=T;i++){
		dis[i]=-INF;
		pre[i]=0;
	}
	dis[S]=0;
	q.push(S);
	while(!q.empty()){
		int x=q.front();
		q.pop();
		inqueue[x]=0;
		for(int i=first[x];i;i=e[i].nxt){
			int v=e[i].v;
			if(e[i].w&&dis[x]+e[i].c>dis[v]){
				pre[v]=i;
				dis[v]=dis[x]+e[i].c;
				if(!inqueue[v]){
					q.push(v);
					inqueue[v]=1;
				}
			}
		}
	}
	return dis[T]!=-INF;
}
bool check(){
	double ans=0;
	while(SPFA()){
		int s=int(1e9+7);
		for(int i=pre[T];i;i=pre[e[i^1].v]){
			s=min(s,e[i].w);
		}
		for(int i=pre[T];i;i=pre[e[i^1].v]){
			e[i].w-=s;
			e[i^1].w+=s;
		}
		ans+=s*dis[T];
//		cout<<"ans= "<<ans<<endl;
	}
	if(ans>0){
		return true;
	} 
	else{
		return false;
	}
}
int main(){
	read(n);
	for(int i=1;i<=n;i++){
		for(int j=1;j<=n;j++){
			int x;
			read(x);
			a[i][j]=(double)(x);
		}
	}
	for(int i=1;i<=n;i++){
		for(int j=1;j<=n;j++){
			int x;
			read(x);
			b[i][j]=(double)x;
		}
	}
	double l=0;
	double r=1e4+7;
	double eps=0.0000001;
	double ans;
	while(r-l>eps){
		double mid=(l+r)/2;
		build(mid);
		if(check()){
			l=mid;
		}
		else{
			r=mid;
		}
//		printf("l=%.10lf r=%.10lf\n",l,r);
	}
	printf("%.6lf",l);
}