計算幾何模板

多圓面積交

typedef long long ll;  

typedef unsigned long long ull; 
typedef vector vi; 
const int inf = 0x3f3f3f3f; 
const double eps = 1e-10; 
const int mod = 100000007; 
const int maxn = 1000010; 
const double pi = acos(-1.0); 
#define sqr(x) ((x)*(x)) 
const int n = 1010; 
double area[n]; 
int n; 
int dcmp(double x) 
struct cp 
cp(double xx, double yy, double ang = 0, int t = 0) 
void get() 
} cir[n], tp[n * 2]; 
double dis(cp a, cp b) 
double cross(cp p0, cp p1, cp p2) 
int circrosscir(cp p1, double r1, cp p2, double r2, cp &cp1, cp &cp2) 
bool circmp(const cp& u, const cp& v) 
bool cmp(const cp& u, const cp& v) 
double calc(cp cir, cp cp1, cp cp2) 
void cirunion(cp cir, int n) 
tp[tn++] = cp(cir[i].x - cir[i].r, cir[i].y, pi, -cnt); 
tp[tn++] = cp(cir[i].x - cir[i].r, cir[i].y, -pi, cnt); 
sort(tp, tp + tn, cmp); 
int p, s = cir[i].d + tp[0].d; 
for (int j = 1; j < tn; ++j) 
} } 
void solve() 
//area[i]為重疊了i次的面積 
//tot 為總面積 
double tot = 0; 
for(int i=1; i<=n; i++) tot += area[i]; 
printf("%f\n", tot); 
} 
int main() 

多圓面積並

#include#include#include#include#include#define ld double

#define eps 1e-13
using namespace std;
int n,top,st,ed;
ld xl[1001],xr[1001];
ld ans;
bool del[1001];
struct datat[1001],sk[1001];
struct linep[1001];
ld dis(data a,data b)
bool cmp1(data a,data b)
sort(p+1,p+sz+1,c***);
for(i=1;i<=sz;i++)
st=i;ed=j-1;i=j-1;
mid=(l+r)/2;
fl=getf(l);fr=getf(r);fmid=getf(mid);
ans+=simpson(l,mid,r,fl,fmid,fr,cal(r-l,fl,fmid,fr));
}}int main()

給你一堆點，找銳角三角形。

two pointer 思想。

統計出所有銳角和直=鈍角的數目。

做法是這樣的：對每個點對所有點極角排序，然後two pointer計算每乙個銳角（乙個邊上有好幾個點也會被統計好幾次），直角鈍角。然後ans=(銳角個數-直角鈍角個數*2)/3;因為每乙個角度可能也只可能出現在乙個三角形中。

但是一上來沒弄懂，所以隊友想了乙個思路也能過，先統計直角和鈍角的個數。然後c（n，3）統計所有情況，不能構成三角形的是三點共線的情況，所以n^2logn統計三點共線（two pointer）。然後再減去直角鈍角個數即可。兩種方法均可過。

#include #include #include #include #include #include #include #include #include #include using namespace std;

typedef long long ll;
struct point
point operator -(const point &a) const
}p[2010];
int n;
ll xmul(const point &a,const point &b)
ll dot(const point &a,const point &b)
bool cmp(const point &a,const point &b)
return xmul(a,b) > 0;
}vectorvec;
int main()
}ans=(rui-2*zhidun)/3;
printf ("%i64d\n",ans);
}return 0;
}

計算幾何模板

sgn返回x經過eps處理的符號，負數返回 1，正數返回1，x的絕對值如果足夠小，就返回0。const double eps 1e 8 int sgn double x double mysqrt double x pt是point的縮寫 int版 struct pt pt int x,int y ...

計算幾何模板

include define vct point using namespace std const double pi atan2 0,1 const double eps 1e 8 int sgn double d struct point bool operator point b const...

計算幾何 模板

include include include include typedef double db const db eps 1e 9 const db pi acos 1 判斷符號 inline intsign db ps 比較大小 inline intcmp db ps,db b 點 struc...