关闭

bzoj2960: 跨平面

标签: bzoj平面图转对偶朱刘算法
486人阅读 评论(2) 收藏 举报

题意:给定平面图,求对偶图的最小树形图。

平面转对偶可以将每条边分成两条,然后暴力找一圈。

最小树形图使用朱刘算法,然而我也不太懂。。

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cmath>
#include<set>
#define N 3005
#define M 30005
#define ll long long
using namespace std;
int n,m,st[M],belong[M],cnt,l;
int id[N],fl[N],fa[N],f[N];
bool ok[M];
struct D{int x,y;}a[N];
struct E{int x,y,p;double k;}b[M];
struct cmp{bool operator()(int x,int y){return b[x].k<b[y].k;}};
set<int,cmp>S[N];

ll cross(D a,D b){return (ll)a.x*b.y-(ll)a.y*b.x;}
void link(int x,int y,int z)
{
	b[++l].x=x;b[l].y=y;b[l].p=z;
}
int solve(int n)
{
	int Ans=0,cnt,rt=n;
	while(1)
	{
		for (int i=0;i<=n;i++)f[i]=1000000000;
		for (int i=1;i<=l;i++)
			if (b[i].p<f[b[i].y]&&b[i].x!=b[i].y)
				f[b[i].y]=b[i].p,fa[b[i].y]=b[i].x;
		f[rt]=0;cnt=-1;
		for (int i=0;i<=n;i++) id[i]=fl[i]=-1;
		for (int i=0;i<=n;i++)
		{
			int j;Ans+=f[i];
			for (j=i;j!=rt&&fl[j]!=i&&id[j]==-1;j=fa[j])fl[j]=i;
			if (j!=rt&&id[j]==-1)
			{
				id[j]=++cnt;
				for (int k=fa[j];k!=j;k=fa[k])
					id[k]=cnt;
			}
		}
		if (cnt==-1) break;
		for (int i=0;i<=n;i++)
			if (id[i]==-1) id[i]=++cnt;
		for (int i=1;i<=l;i++)
		{
			int t=f[b[i].y];
			b[i].x=id[b[i].x];
			b[i].y=id[b[i].y];
			if (b[i].x!=b[i].y) b[i].p-=t;
		}
		n=cnt;rt=id[rt];
	}
	return Ans;
}
int main()
{
	scanf("%d%d",&n,&m);
	for (int i=1;i<=n;i++)
		scanf("%d%d",&a[i].x,&a[i].y);
	for (int i=1;i<=m;i++)
	{
		int t=i<<1;
		scanf("%d%d%d%d",&b[t].x,&b[t].y,&b[t].p,&b[t|1].p);
		b[t|1].x=b[t].y;b[t|1].y=b[t].x;
	}
	m=m*2+1;
	for (int i=2;i<=m;i++)
	{
		b[i].k=-atan2(a[b[i].y].y-a[b[i].x].y,a[b[i].y].x-a[b[i].x].x);
		S[b[i].x].insert(i);
	}
	set<int,cmp>::iterator i1;
	for (int i=2;i<=m;i++)
		if (!ok[i])
		{
			int top=1;ll sum=0;st[1]=i;
			while(1)
			{
				ok[st[top]]=1;
				i1=S[b[st[top]].y].find(st[top]^1);i1++;
				if (i1==S[b[st[top]].y].end())
					i1=S[b[st[top]].y].begin();
				if (*i1==i)break;
				st[++top]=*i1;
			}
			for (int j=1;j<=top;j++)
				sum+=cross(a[b[st[j]].x],a[b[st[j]].y]);
			if (sum>=0)
			{
				cnt++;
				for (int j=1;j<=top;j++) belong[st[j]]=cnt;
			}
		}
	n=cnt+1;int x=1;
	for (int i=2;i<=m;i++)
		if (b[i].p) x+=b[i].p,link(belong[i],belong[i^1],b[i].p);
	for (int i=0;i<n;i++) link(n,i,x);
	printf("%d\n",solve(n)-x);
	return 0;
}


0
0
查看评论
发表评论
* 以上用户言论只代表其个人观点,不代表CSDN网站的观点或立场

bzoj 2960: 跨平面 最小树形图 朱-刘算法

首先用set找出所有的区域,然后对于一条边在相邻的区域之间连有向边,然后新建一个根连向所有点跑朱-刘算法即可。 AC代码如下: #include #include #include #include ...
  • lych_cys
  • lych_cys
  • 2016-04-14 16:13
  • 594

bzoj 2960 跨平面 平面图转对偶图 最小树形图

平面图转对偶图就是枚举每一条没选过的边,每次找反边顺时针方向最近的一条边,这样就可以找到每一个简单多边形(包括最大的多边形) 这个可以用set维护。注意找完后一条边后不要删!!!!!!!!!!!然后...
  • make_it_for_good
  • make_it_for_good
  • 2016-10-24 18:14
  • 221

【BZOJ2960】跨平面(对偶图+朱刘算法)

DescriptionInput第一行两个整数n和m,表示点与线段的数目。 接下来n行,每行两个整数x和y,表示第i个点的坐标,点从1到n编号。 接下来m行,每行四个整数p,q,V1和V2,表示存...
  • mousemjq
  • mousemjq
  • 2017-03-29 17:38
  • 221

分割平面的几种方法

(1) n条直线最多分平面问题        题目大致如:n条直线,最多可以把平面分为多少个区域。        析:可能你以前就见过这题目,这充其量是一道...
  • long71751380
  • long71751380
  • 2013-07-30 18:21
  • 2349

PCL系列——平面模型分割

PCL系列 PCL系列——读入PCD格式文件操作 PCL系列——将点云数据写入PCD格式文件 PCL系列——拼接两个点云 PCL系列——从深度图像(RangeImage)中提取NARF关键点 PCL系...
  • xuezhisdc
  • xuezhisdc
  • 2016-04-20 11:12
  • 2793

PCL中RANSAC使用,点云平面检测,显示,存储

#include #include #include #include #include #include #include #include #include #include ...
  • fandq1223
  • fandq1223
  • 2016-11-24 16:35
  • 3231

AR demo 简单的平面图像识别

1. 首先到Vuforia的官网下载SDK,安装好最新版本的unity 以及下载好 android SDK。 我这里用到的是5-5-9版本。 2. 在官网的develop分栏里配置...
  • NeoAdrian
  • NeoAdrian
  • 2016-07-27 21:43
  • 3320

LintCode 22 平面列表

题目:flatten要求:给定一个列表,该列表中的每个要素要么是个列表,要么是整数。将其变成一个只包含整数的简单列表。 注意事项如果给定的列表中的要素本身也是一个列表,那么它也可以包含列表。样例:给定...
  • linglian0522
  • linglian0522
  • 2017-02-22 16:27
  • 720

利用 ransac 算法拟合平面

1.前言最近项目中遇到一个问题, 老板给了一组数据然后要求获取其中处于同一个平面上的数据点的信息, 很明显就是使用ransac 算法进行处理。2. ransac算法思想这里我们使用自己的理解来说明下这...
  • zhyh1435589631
  • zhyh1435589631
  • 2016-11-01 10:16
  • 2266

空间离散点拟合成空间平面

1.原理 空间中的离散点得到拟合平面,其实这就是一个最优化的过程。即求这些点到某个平面距离和最小的问题。我们知道一个先验消息,那就是该平面一定会过众散点的平均值。接着我们需要做的工作就是求这个平面的...
  • shenziheng1
  • shenziheng1
  • 2016-04-17 21:11
  • 4398
    个人资料
    • 访问:40712次
    • 积分:1346
    • 等级:
    • 排名:千里之外
    • 原创:93篇
    • 转载:0篇
    • 译文:0篇
    • 评论:19条
    友情链接
    最新评论