POJ 2546 Circular Area

POJ 2546 Circular Area

题目链接：POJ 2546

题意：

求两个圆相交的面积。

分析：

将 两个圆的面积 减去 两个圆并的面积 就是 答案了。（主要是为了应用前面提及的模板）

实现1：

#include<cstdio>
#include<cmath>
#include<cstring>
#include<iostream>
#include<vector>
#include<algorithm>
using namespace std;

const double eps = 5 * 1e-13;
int dcmp(double x) 
{
  if(fabs(x) < eps) 
	  return 0; 
  else 
	  return x < 0 ? -1 : 1;
}

const double PI = acos(-1.0);
const double TWO_PI = PI * 2;

double NormalizeAngle(double rad, double center = PI) 
{
  return rad - TWO_PI * floor((rad + PI - center) / TWO_PI);
}

struct Point 
{
  double x, y;
  Point(double x=0, double y=0):x(x),y(y) { }
};

typedef Point Vector;

Vector operator + (Vector A, Vector B)
{ 
	return Vector(A.x+B.x, A.y+B.y); 
}
Vector operator - (Point A, Point B) 
{ 
	return Vector(A.x-B.x, A.y-B.y); 
}
Vector operator * (Vector A, double p) 
{ 
	return Vector(A.x*p, A.y*p); 
}
Vector operator / (Vector A, double p) 
{ 
	return Vector(A.x/p, A.y/p); 
}

// 理论上这个“小于”运算符是错的，因为可能有三个点a, b, c, a和b很接近（即a<b好b<a都不成立），b和c很接近，但a和c不接近
// 所以使用这种“小于”运算符的前提是能排除上述情况
bool operator < (const Point& a, const Point& b) 
{
  return dcmp(a.x - b.x) < 0 || (dcmp(a.x - b.x) == 0 && dcmp(a.y - b.y) < 0);
}

bool operator == (Point A, Point B) 
{
  return dcmp(A.x - B.x) == 0 && dcmp(A.y - B.y) == 0;
}

double Dot(Vector A, Vector B) 
{ 
	return A.x*B.x + A.y*B.y; 
}
double Length(Vector A) 
{ 
	return sqrt(Dot(A, A)); 
}
double Cross(Vector A, Vector B) 
{ 
	return A.x*B.y - A.y*B.x; 
}

double angle(Vector v) 
{
  return atan2(v.y, v.x);
}

// 交点相对于圆1的极角保存在rad中
void getCircleCircleIntersection(Point c1, double r1, Point c2, double r2, vector<double>& rad) 
{
  double d = Length(c1 - c2);
  if(dcmp(d) == 0) return; // 不管是内含还是重合，都不相交
  if(dcmp(r1 + r2 - d) < 0) return;
  if(dcmp(fabs(r1-r2) - d) > 0) return;

  double a = angle(c2 - c1);
  double da = acos((r1*r1 + d*d - r2*r2) / (2*r1*d));
  rad.push_back(NormalizeAngle(a-da));
  rad.push_back(NormalizeAngle(a+da));
}

Point GetLineProjection(Point P, Point A, Point B) 
{
  Vector v = B-A;
  return A+v*(Dot(v, P-A) / Dot(v, v));
}
/// 题目相关
const int maxn = 200 + 5;
int n, N, M; // n是圆的总数，N和M是场地长宽
Point P[maxn];
double R[maxn];

// 取圆no弧度为rad的点
Point getPoint(int no, double rad) 
{
  return Point(P[no].x + cos(rad)*R[no], P[no].y + sin(rad)*R[no]);
}

// 场地边界上的点p是否可见
bool visible(Point p) 
{
  for(int i = 0; i < n; i++)
  {
    if(dcmp(Length(p - P[i]) - R[i]) < 0) return false;
  }
  return true;
}

// 求圆的并在(0,0)-(N,M)内的面积
// 使用一般曲边图形的面积算法。下文中，“所求图形”指的是不能种菜的区域，它的边界由圆弧和直线段构成。
// 算法：对于所求图形边界上的每一段（可以是曲线）a~>b，累加Cross(a, b)和它在直线段a->b右边部分的面积（左边部分算负）
// 边界计算：
// 1. 每个圆被其他圆和场地边界分成了若干条圆弧，中点不被其他圆覆盖且在场地内的圆弧在所求图形边界上
// 2. 场地的四条边界被圆分成了若干条线段。中点在某个圆内部的线段在所求图形边界上
double getArea() 
{
  double area = 0;
  // 圆弧部分
  for(int i = 0; i < n; i++)
  {
    vector<double> rad;
    rad.push_back(0);
    rad.push_back(PI*2);

    // 圆和圆的交点
    for(int j = 0; j < n; j++)
      getCircleCircleIntersection(P[i], R[i], P[j], R[j], rad);
    
    sort(rad.begin(), rad.end());
    for(int j = 0; j < rad.size()-1; j++) 
		if(rad[j+1] - rad[j] > eps) 
		{
		  double mid = (rad[j] + rad[j+1]) / 2.0; // 圆弧中点相对于圆i圆心的极角
		  if(visible(getPoint(i, mid)))
		  { // 弧中点可见，因此弧在图形边界上
			area += Cross(getPoint(i, rad[j]), getPoint(i, rad[j+1])) / 2.0;
			double a = rad[j+1] - rad[j];
			area += R[i] * R[i] * (a - sin(a)) / 2.0;
		  }
		}
   }

  // 直线段部分
  for(int i = 0; i < 4; i++) 
  {
    vector<double> dist;
    dist.push_back(0);
    for(int j = 0; j < n; j++) 
	{
      vector<double> rad;
      for(int k = 0; k < rad.size(); k++) 
	  {
        Point p = getPoint(j, rad[k]);
      }
    }

    sort(dist.begin(), dist.end()); // 必须按照到起点的距离排序而不是按照点的字典序排序，否则向量方向可能会反
    vector<Point> points;

    for(int j = 0; j < dist.size()-1; j++) 
	{
      Point midp = (points[j] + points[j+1]) / 2.0;
      if(!visible(midp)) 
		  area += Cross(points[j], points[j+1]) / 2.0; // 线段中点不可见，因此线段在图形边界上
    }
  }

  return area;
}

int main()
{
	// freopen("in.txt","r",stdin);
  while(scanf("%lf %lf %lf %lf %lf %lf", &P[0].x, &P[0].y, &R[0], &P[1].x, &P[1].y, &R[1]) !=EOF)
  {
    sort(P, P+2);
    n = unique(P, P+2) - P;
    printf("%.3lf\n", R[0]*R[0]*PI+R[1]*R[1]*PI-getArea());
  }
  return 0;
}

实现2（针对两个圆的）

#include <stdio.h>
#include <string.h>
#include <iostream>
#include <algorithm>
#include <vector>
#include <queue>
#include <set>
#include <map>
#include <string>
#include <math.h>
#include <stdlib.h>
#include <time.h>
using namespace std;

double Area(double x,double y,double r,double xx,double yy,double rr)  
{  

	double d, s, t, t1;
	d=sqrt((x-xx)*(x-xx)+(y-yy)*(y-yy));  
	if(d>=r+rr) s=0;  
	else if(d<=fabs(r-rr))
		s=min(acos(-1.0)*r*r,acos(-1.0)*rr*rr);  
	else    //（里面有两种情况可以统一） 
	{
		t=(r*r+d*d-rr*rr)/2.0/d;  
		t1=sqrt(r*r-t*t);  
		s=-d*t1+r*r*acos(t/r)+rr*rr*acos((d-t)/rr);  
	}
	return s;
	//cout<<fixed<<setprecision(3)<<s<<endl;  
} 
 
int main()
{
    double x,y,r,xx,yy,rr;
	while(scanf("%lf %lf %lf %lf %lf %lf",&x,&y,&r,&xx,&yy,&rr)!=EOF)
	{
		printf("%.3lf\n",Area(x,y,r,xx,yy,rr));
	}
    return 0;
}