2020icpc 上海站 （计算几何）

题目链接：I-Sky Garden_第 45 届国际大学生程序设计竞赛（ICPC）亚洲区域赛（上海） (nowcoder.com)

题意：有n个圆他们的圆心都是（0,0），并且第i个圆的半径是i。然后有m条直线将这些圆等分成了2*m分。并且产生了一些交点。让你求任意两个不同的交点之间的最短距离之和是多少。

首先我们要特判的一种情况就是所有点中是否包含原点，显然如果直线条数大于等于2，那么一定是要包含原点的。我们下面就着重讨论直线条数大于等于2的情况。

任意两个不同的交点之间的最短距离之和，我们可以把原点单独讨论，原点到其他点之间的最短距离之和很好求，我们可以单独求出原点到其他点的最短距离之和，除了原点之外，其他点之间只有两种关系，一种是在同一个圆上，还有一种就是不在同一个圆上，我们可以暴力求出半径为1的圆上不同点之间的最短距离之和，不妨假设值为ans，那么半径为r的圆上不同点之间的最短距离之和等于多少呢？显然是等于r*ans，是与半径成正比的。现在我们已经解决了同一个圆上不同的点之间的最小距离和，剩下的问题就是求不在同一个圆上的点之间的最短距离和了，半径为r的圆上的所有点？到半径为r+1的圆上的所有点的最短距离和是多少呢是ans*2+cnt*cnt*r,cnt是圆上点的个数，从半径为r的圆上的一个点到半径为r+1的圆上一个点的最短路线是什么样的呢?不妨我们假设从半径为r的圆上的第i个点移动到半径为r+1的圆上的第j个点，那么最短路径就是先从半径为r的圆上的第i个点移动到半径为r的圆上的第j个点，再移动1就可以了，由于在两个圆上选取的点是相互独立的，所以一共会多走cnt*cnt次半径，由于我们每次移动都是在当前圆上进行的，所以我们要加上ans*2，半径差别不是1的两个圆分析也是一样的，我说的可能不是很好理解，建议在图上画一画，这样会比较清楚，下面是代码：

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<iostream>
#include<queue>
#include<cmath>
using namespace std;
#define pi acos(-1.0)
int main()
{
	int n,m;
	cin>>n>>m;
	int cnt=2*m;
	double o=pi/m;
	double ans=0;
	for(int i=1;i<cnt;i++)
		ans+=min(2.0,o*min(i,cnt-i));//计算半径最小的圆上1点到其他点的最小距离和 
	ans=ans*cnt/2;//计算半径最小的圆上的点之间的最小距离和 
	double tans=0;//记录答案 
	for(int i=1;i<=n;i++)
		tans+=i*ans;//计算每个圆上的点到该圆上其他点的最小距离和 
	if(m!=1)//原点存在贡献 
		tans+=cnt*(n+1)*n/2;
	//计算每个圆上的点到半径比他大的圆上各点的最小距离和 
	for(int i=1;i<=n;i++)
	for(int j=i+1;j<=n;j++)
		tans+=i*ans*2+(j-i)*cnt*cnt;
	printf("%.10lf",tans);
	return 0;
 }