2020ICPC上海站 I.Sky Garden（计算几何+思维）

题目链接：https://ac.nowcoder.com/acm/contest/9925/I

题意：

        画n个圆和m条直线，圆的中心点为(0,0)，圆的半径分别从1到n，而直线都必经过(0,0)点，并且所有直线会把每个圆平均分成2m个面积相等的区域，直线会和圆形成交点，求所有交点两两经过线的最短路径的和。

思路：

1. 结论一：同一圆上两点的最短路径是两点通过圆弧的路径或者是两点与圆心连线的距离和，即图中A到B的最短路径为红色路径或棕色路径。

2. 结论二：两个不在一个圆上的点之间的最短路径是外圈圆的点通过直线到达内圈圆的点所在的圆上，然后再通过到达的点与圆上另一个点求最短路径。即如下图B到A的最短路径是B到C后，再加上C到A的最短路径。

1. 首先要特判，如果m=1，则不存在中心交点。如果m!=1，那么就可以先计算中心交点到各点的路径和，公式为(1+n)*n*m。

2. 接着从最内部的圆往最外部的圆进行计算，每执行一次计算用sum记录此圆内某个点到达圆内其他2m-1个点的距离和（在此过程中要利用结论一进行判断），为了计算当前圆内的两两各点之间的最短距离和，距离和为sum*2*m-sub，所以sub记录的是重复计算的数据，需要一直更新。

3. 完成了圆内两两各点的最短距离和后，则需要计算圆内各点到达外层圆上各点的最短距离和，根据结论二则可得到圆内某个点到达外层圆上的所有点的距离和为m*(1+(n-i))*(n-i)+sum*(n-i)，其中m*(1+(n-i))*(n-i)计算的是当前第i层圆上某点到达外层圆上各点时走到直线上的路程，sum*(n-i)则计算的是圆内某点到达外层圆上各点投影到当前圆上各点的距离和，最后要乘上2m，因为圆上有2m个点。

代码：

#include<bits/stdc++.h>
using namespace std;
const double pi = acos(-1);

int main()
{
	int n,m;
	scanf("%d %d",&n,&m);
	double ans=0;
	if(m!=1) ans=ans+(1+n)*n*m;
	for(int i=1;i<=n;i++){
		double sub=0;
		double sum=0;
		for(int j=1;j<=m;j++){
			if(j==m){
				sub=sub+m*2*i;
				sum=sum*2+2*i;
			}
			else{
				if((pi*i*j)/m<2*i){
					sum=sum+(pi*i*j)/m;
					sub=sub+2*m*((pi*i*j)/m);
				}
				else{
					sum=sum+2*i;
					sub=sub+2*m*2*i;
				}
			}
		}
		ans+=sum*2*m-sub;
		ans=ans+(m*(1+(n-i))*(n-i)+sum*(n-i))*2*m;
	}
	printf("%.10f\n",ans);
	return 0;
}