版权声明：本文为博主原创文章，未经博主允许不得转载。 https://blog.csdn.net/hanyuweining/article/details/84674403
想出的一道题竟然是原题QAQ
非常有趣的一个题
根据三角形两边之和大于第三边 所以相交的线段一定是比不相交的线段要长的
所以直接二分图构图 最小费用最大流即可
（我不管我不管我要把这个出到NOIP膜你赛）
代码如下。
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<queue>
#define inf 20021225
#define ll long long
#define db double
#define mxn 400
#define eps 1e-6
using namespace std;
struct edge{int to,lt,f,fr;db c;}e[mxn*mxn];
struct point{int x,y;}p[mxn];
int in[mxn],from[mxn],cnt=1,nn,s,t,n,m;
db dis[mxn],ans;bool vis[mxn];
void add(int x,int y,int f,db c)
{
e[++cnt].to=y;e[cnt].lt=in[x];e[cnt].f=f;e[cnt].fr=x;e[cnt].c=c;in[x]=cnt;
e[++cnt].to=x;e[cnt].lt=in[y];e[cnt].f=0;e[cnt].fr=y;e[cnt].c=-c;in[y]=cnt;
}
queue<int> que;
bool spfa()
{
for(int i=1;i<=nn;i++) dis[i]=inf,vis[i]=0;
vis[s]=1;dis[s]=0.0;que.push(s);
while(!que.empty())
{
int x=que.front();que.pop();vis[x]=0;
//printf("%d ",x);
for(int i=in[x];i;i=e[i].lt)
{
int y=e[i].to;
if(e[i].f&&dis[y]>dis[x]+e[i].c)
{
dis[y]=dis[x]+e[i].c;
from[y]=i;
if(!vis[y]) que.push(y),vis[y]=1;
}
}
}
//printf("%lf %lf\n",dis[t],inf+eps);
return dis[t]<inf-eps;
}
db flow()
{
int flow=inf;
for(int i=t;i!=s;i=e[from[i]].fr) flow=min(flow,e[from[i]].f);
for(int i=t;i!=s;i=e[from[i]].fr) e[from[i]].f-=flow,e[from[i]^1].f+=flow;
return flow*dis[t];
}
void dinic()
{
while(spfa()) ans+=flow();
}
db dist(point a,point b)
{
return sqrt((db)(a.x-b.x)*(a.x-b.x)+(db)(a.y-b.y)*(a.y-b.y));
}
int main()
{
while(scanf("%d",&n)!=EOF)
{
cnt=1;memset(in,0,sizeof(in));
s=n*2+1;t=s+1;nn=t;
for(int i=1;i<=2*n;i++) scanf("%d%d",&p[i].x,&p[i].y);
for(int i=1;i<=n;i++) add(s,i,1,0),add(i+n,t,1,0);
for(int i=1;i<=n;i++)
for(int j=1;j<=n;j++)
add(i,j+n,1,dist(p[i],p[j+n]));
dinic();
for(int i=1;i<=n;i++)
for(int j=in[i];j;j=e[j].lt)
if(!e[j].f)
{
printf("%d\n",e[j].to-n);
break;
}
}
return 0;
}