题目背景
Grant喜欢带着他的小狗Pandog散步。Grant以一定的速度沿着固定路线走，该路线可能自交。Pandog喜欢游览沿途的景点，不过会在给定的N个点和主人相遇。小狗和主人同时从(X1，Y1)点出发，并同时在(Xn，Yn)点汇合。小狗的速度最快是Grant的两倍。当主人从一个点以直线走向另一个点时，Pandog跑向一个它感兴趣的景点。Pandog每次与主人相遇之前最多只去一个景点。
题目描述
你现在的任务是：为Pandog寻找一条路线(有可能与主人的路线部分相同)，使它能够游览最多的景点，并能够准时与主人在给定地点相遇或者汇合。
输入格式
输入文件第一行是两个整数N和M( 1≤N，M≤100 )；
输入文件第二行的N个坐标给出了Grant的散步路线，即Pandog和主人相遇地点；
输入文件第三行的M个坐标给出了所有Pandog感兴趣的景点。
所有输入的坐标均不相同，且绝对值不超过1000。
输出格式
输出小狗的移动路线。
第一行是经过的点数，第二行依次为经过的点的坐标(直角坐标系)
输入输出样例
输入 #1复制
4 5
1 4 5 7 5 2 -2 4
-4 -2 3 9 1 2 -1 3 8 -3
输出 #1复制
6
1 4 3 9 5 7 5 2 1 2 -2 4
说明/提示
"The way is wrong!"表示输出方案错误（可能是坐标不存在输入文件中，两个相遇点间存在多个景点，或距离超出范围）
因为两个景点之间的距离已知，所以我们可以求出这两个景点之间可以有哪些景点是小狗可以到达的。
所以其实是一张二分图。
扫描二维码关注公众号，回复：
8515276 查看本文章
怎么输出方案呢？从汇点T开始跑一下，判断流量即可。
AC代码：
#pragma GCC optimize("-Ofast","-funroll-all-loops")
#include<bits/stdc++.h>
//#define int long long
using namespace std;
const int inf=0x3f3f3f3f;
const int N=210,M=1e5+10;
int n,m,h[N],s,t;
int head[N],nex[M],to[M],w[M],tot=1;
struct node{int x,y;};
inline void ade(int a,int b,int c){
to[++tot]=b; nex[tot]=head[a]; w[tot]=c; head[a]=tot;
}
inline void add(int a,int b,int c){ade(a,b,c); ade(b,a,0);}
inline double dis(node a,node b){
return sqrt((a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y));
}
inline int bfs(){
queue<int> q; q.push(s); memset(h,0,sizeof h); h[s]=1;
while(q.size()){
int u=q.front(); q.pop();
for(int i=head[u];i;i=nex[i]){
if(w[i]&&!h[to[i]]){
h[to[i]]=h[u]+1; q.push(to[i]);
}
}
}
return h[t];
}
int dfs(int x,int f){
if(x==t) return f; int fl=0;
for(int i=head[x];i&&f;i=nex[i]){
if(w[i]&&h[to[i]]==h[x]+1){
int mi=dfs(to[i],min(w[i],f));
w[i]-=mi,w[i^1]+=mi,fl+=mi,f-=mi;
}
}
if(!fl) h[x]=-1;
return fl;
}
inline int dinic(){
int res=0;
while(bfs()) res+=dfs(s,inf);
return res;
}
signed main(){
cin>>n>>m; t=n+m+1; node p[N],g[N];
for(int i=1;i<=n;i++) cin>>g[i].x>>g[i].y;
for(int i=1;i<=m;i++) cin>>p[i].x>>p[i].y,add(i+n,t,1);
for(int i=2;i<=n;i++){
add(s,i,1);
for(int j=1;j<=m;j++){
if(dis(g[i-1],p[j])+dis(g[i],p[j])<=dis(g[i-1],g[i])*2.0)
add(i,j+n,1);
}
}
cout<<n+dinic()<<'\n';
for(int i=1;i<=n;i++) cout<<g[i].x<<' '<<g[i].y<<' ';
for(int i=head[t];i;i=nex[i]){
if(w[i]) cout<<p[to[i]-n].x<<' '<<p[to[i]-n].y<<' ';
}
return 0;
}
