luogu P4515 [COCI2009-2010#6] XOR 容斥
LINK:XOR
一个不常见的容斥套路题。
以往是只求三角形面积的交 现在需要求被奇数次覆盖的区域的面积。
打住 求三角形面积的交我也不会写 不过这道题的三角形非常特殊 等腰直角 且直角点都在左下方 这就有很多的性质了。
容易发现最后交出的三角形为等腰直角三角形。
考虑如何求若干个三角形交出的面积 不太会证明 题解区的一个神仙给出了一个式子。
设 \(c_i=x_i+y_i+z_i\)最终交出的三角形的直角边边长为 \(MAX(0,min(c_i)-max(x_i)-max(y_i))\)
数据范围这么小 显然可以子集容斥 不过对于枚举到的三角形 需要配上一定的容斥系数满足 偶消奇不消。
对于一个集合s来说 容斥系数为\(2^{|S|-1}(-1)^{|S|-1}\)
怎么说 这是 对于这种容斥的常用套路(系数。
证明:\(\sum_{k=1}^nC(n,k)2^{k-1}(-1)^{k-1}=[![2|n]]\)
\(\sum_{k=1}^nC(n,k)(-2)^{k-1}=\frac{\sum_{k=1}^nC(n,k)(-2)^{k}}{-2}=\frac{-1+\sum_{k=0}^nC(n,k)(-2)^{k}}{-2}\)
二项式定理合并起来 可得\(\frac{1-(-2+1)^n}{2}=\frac{1-(-1)^n}{2}=[![2|n]]\)
const int MAXN=12;
int n;
struct wy
{
int x,y,r,w;
}t[MAXN];
db ans;
inline void dfs(int v,int sz,int z,int x,int y,int op)
{
if(v==n+1)
{
if(!sz)return;
ans=ans+(1ll<<sz-1)*op*((z-x-y)<0?0:(ll)(z-x-y)*(z-x-y));
return;
}
dfs(v+1,sz+1,min(z,t[v].w),max(x,t[v].x),max(y,t[v].y),-op);
dfs(v+1,sz,z,x,y,op);
}
int main()
{
freopen("1.in","r",stdin);
get(n);
rep(1,n,i)
{
int x,y,z;
get(x);get(y);get(z);
t[i]=(wy){x,y,z};
t[i].w=x+y+z;
}
dfs(1,0,INF,0,0,-1);
printf("%.1lf",ans/2);
return 0;
}