P1500 丘比特的烦恼

题目描述

随着社会的不断发展，人与人之间的感情越来越功利化。最近，爱神丘比特发现，爱情也已不再是完全纯洁的了。这使得丘比特很是苦恼，他越来越难找到合适的男女，并向他们射去丘比特之箭。于是丘比特千里迢迢远赴中国，找到了掌管东方人爱情的神——月下老人，向他求教。

月下老人告诉丘比特，纯洁的爱情并不是不存在，而是他没有找到。在东方，人们讲究的是缘分。月下老人只要做一男一女两个泥人，在他们之间连上一条红线，那么它们所代表的人就会相爱——无论他们身处何地。而丘比特的爱情之箭只能射中两个距离相当近的人，选择的范围自然就小了很多，不能找到真正的有缘人。

丘比特听了月下老人的解释，茅塞顿开，回去之后用了人间的最新科技改造了自己的弓箭，使得丘比特之箭的射程大大增加。这样，射中有缘人的机会也增加了不少。

情人节(Valentine’s day)的午夜零时，丘比特开始了自己的工作。他选择了一组数目相等的男女，感应到他们互相之间的缘分大小，并依此射出了神箭，使他们产生爱意。他希望能选择最好的方法，使被他选择的每一个人被射中一次，且每一对被射中的人之间的缘分的和最大。

当然，无论丘比特怎么改造自己的弓箭，总还是存在缺陷的。首先，弓箭的射程尽管增大了，但毕竟还是有限的，不能像月下老人那样，做到“千里姻缘一线牵”。其次，无论怎么改造，箭的轨迹终归只能是一条直线，也就是说，如果两个人之间的连线段上有别人，那么莫不可向他们射出丘比特之箭，否则，按月下老人的话，就是“乱点鸳鸯谱”了。

作为一个凡人，你的任务是运用先进的计算机为丘比特找到最佳的方案。

输入格式

输入第一行为正整数k，表示丘比特之箭的射程，第二行为正整数n(n<30)，随后有2n行，表示丘比特选中的人的信息，其中前n行为男子，后n行为女子。每个人的信息由两部分组成：他的姓名和他的位置。姓名是长度小于20且仅包含字母的字符串，忽略大小写的区别，位置是由一对整数表示的坐标，它们之间用空格分隔。格式为x y Name。输入文件剩下的部分描述了这些人的缘分。每一行的格式为Name1 Name2 p。Name1和Name2为有缘人的姓名，p是他们之间的缘分值(p为小于等于255的正整数)。以一个End作为文件结束标志。每两个人之间的缘分至多只被描述一次。如果没有被描述，则说明他们缘分值为1。

输出格式

输出仅一个正整数，表示每一对被射中的人之间的缘分的总和。这个和应当是最大的。

---

题解：
比较明显的最大费用最大流，就是要注意的地方比较多

1.名字不区分大小写，需要全部转换成小写或者大写
2.洛谷明确标注了每两个人缘分最多被描述一次，而BZOJ会有这种数据，直接覆盖直接的缘分值即可
3.一定不能判断前三个字母是End就退出，会有些数据的人名前三个字母就是End…
4.两点间是否有一个点满足三点一线随便找个板子就好了
5.没有被描述的两个人缘分值为1

注意到这些建图完了直接跑一遍EK就OK了

---

AC代码：

#pragma GCC optimize(2)
#include<bits/stdc++.h>
#include<ext/rope>
using namespace std;
using namespace __gnu_cxx;
#define LL long long
const int MAXN = 2000+50;
const int MAXM = 2e6+50;
const int MOD = 1e9+7;
const int INF = 0x3f3f3f3f;
int n,d,s,t,tot=1,maxflow,res,head[MAXN],to[MAXM],w[MAXM],nxt[MAXM],co[MAXM];
int dis[MAXN],vis[MAXN],pre[MAXN],flow[MAXN],g[50][50];
struct node{ int x,y; string ss; }a[MAXN];
inline void deal(string &s){
    for(int i=0;i<(int)s.length();i++) s[i]=tolower(s[i]);
}
inline double distance(int x,int y,int xx,int yy){
    return sqrt((x-xx)*(x-xx)+(y-yy)*(y-yy));
}
inline bool check(int x,int y){
    if(distance(a[x].x,a[x].y,a[y].x,a[y].y)>d) return false;
    for(int i=1;i<=2*n;i++){
        if(i==x || i==y) continue;
        double d1=distance(a[x].x,a[x].y,a[i].x,a[i].y);
        double d2=distance(a[y].x,a[y].y,a[i].x,a[i].y);
        double d3=distance(a[x].x,a[x].y,a[y].x,a[y].y);
        if(d1+d2-d3<1e-9) return false;
    }
    return true;
}
inline void ade(int u,int v,int ww,int cost){
    to[++tot]=v; w[tot]=ww; co[tot]=cost; nxt[tot]=head[u]; head[u]=tot;
}
inline void add(int u,int v,int w,int cost){ ade(u,v,w,cost); ade(v,u,0,-cost); }
inline int spfa(){
    for(int i=0;i<=t;i++) dis[i]=-INF;
    queue<int> que; que.push(s); dis[s]=0; vis[s]=1; flow[s]=INF;
    while(!que.empty()){
        int u=que.front(); que.pop(); vis[u]=0;
        for(int i=head[u];i;i=nxt[i]){
            if(w[i] && dis[to[i]]<dis[u]+co[i]){
                dis[to[i]]=dis[u]+co[i]; pre[to[i]]=i;
                flow[to[i]]=min(flow[u],w[i]);
                if(!vis[to[i]]) vis[to[i]]=1,que.push(to[i]);
            }
        }
    }
    return dis[t]!=-INF;
}
inline void EK(){
    while(spfa()){
        res+=flow[t]*dis[t]; int x=t;
        while(x!=s){
            int i=pre[x];
            w[i]-=flow[t]; w[i^1]+=flow[t]; x=to[i^1];
        }
    }
}
signed main(){
#ifndef ONLINE_JUDGE
    freopen("C:\\Users\\Administrator\\Desktop\\in.txt","r",stdin);
#endif // ONLINE_JUDGE
    scanf("%d%d",&d,&n); s=0,t=2*n+1;
    for(int i=1;i<=2*n;i++) cin>>a[i].x>>a[i].y>>a[i].ss,deal(a[i].ss);
    string s1,s2; int k;
    while(cin>>s1){
        if(s1=="End") break;
        cin>>s2>>k; int p1=0,p2=0;
        deal(s1); deal(s2);
        for(int i=1;i<=2*n;i++) if(a[i].ss==s1) p1=i;
        for(int i=1;i<=2*n;i++) if(a[i].ss==s2) p2=i;
        if(p1>p2) swap(p1,p2); g[p1][p2]=k;
    }
    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=n+1;j<=2*n;j++){
            if(check(i,j)){
                if(g[i][j]) add(i,j,1,g[i][j]);
                else add(i,j,1,1);
            }
        }
    }
    EK();
    printf("%d\n",res);
    return 0;
}