洛谷P2526 小狗散步 二分图匹配

给出 n n n个地点,表示主人的一条路线。然后又给出了 m m m个景点,是小狗想要去的地方。小狗的速度不超过主人的两倍,小狗必须和主人在所有的地点汇合,而且两次汇合之间小狗最多只能去一个地点。 n , m ≤ 1 e 2 n,m\leq1e2 n,m1e2
由于小狗最多只能去一个景点,我们不妨直接预处理出所有从当前点 u u u出发,能到到达景点 x x x,并且还能回到 v v v的所有的点,这样就是一个二分图匹配的问题。
直接匹配并且同时记录下连边的点就可以了。

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int inf=0x3f3f3f3f;
const ll INF=LONG_LONG_MAX;
const int N=202;
int n,m; 
struct Point {
	double x,y;
}r[N],a[N];
vector<int> G[N];
int lk1[N],lk2[N];
int used[N],tim=0;

double sq(double x) {
	return x*x;
}
double dis(Point a,Point b) {
	return sqrt(sq(a.x-b.x)+sq(a.y-b.y));
}
int dfs(int u) {
	for(auto &v:G[u]) {
		if(used[v]==tim) continue;
		used[v]=tim;
		if(lk2[v]==-1||dfs(lk2[v])) {
			lk2[v]=u;
			lk1[u]=v;
			return 1;
		}
	}
	return 0;
}
int solve() {
	memset(lk1,-1,sizeof(lk1));
	memset(lk2,-1,sizeof(lk2));
	tim=0;
	int ans=0;
	for(int i=1;i<n;i++) {
		tim++;
		if(dfs(i)) ++ans;
	} 
	return ans;
}
int main() {
	scanf("%d%d",&n,&m);
	for(int i=1;i<=n;i++)
		scanf("%lf%lf",&r[i].x,&r[i].y);
	for(int i=1;i<=m;i++)
		scanf("%lf%lf",&a[i].x,&a[i].y);
	for(int i=1;i<n;i++) {
		for(int j=1;j<=m;j++) {
			if(dis(r[i],a[j])+dis(a[j],r[i+1])<=2.0*dis(r[i],r[i+1])) {
				G[i].push_back(j);
			}
		}
	}
	int ans=solve();
	ans+=n;
	printf("%d\n",ans);
	for(int i=1;i<=n;i++) {
		if(i<n) printf("%.0lf %.0lf ",r[i].x,r[i].y);
		else printf("%.0lf %.0lf\n",r[n].x,r[n].y);
		if(lk1[i]!=-1) printf("%.0lf %.0lf ",a[lk1[i]].x,a[lk1[i]].y);
	}
	return 0;
}
  • 0
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值