luogu P4515 [COCI2009-2010#6] XOR

描述

坐标系下有若干个等腰直角三角形,且每个等腰直角三角形的直角顶点都在左下方,两腰与坐标轴平行。被奇数个三角形覆盖的面

积部分为灰色,被偶数个三角形覆盖的面积部分为白色,如下图所示。

已知 NN个等腰直角三角形的顶点坐标及腰长,求灰色部分面积。

输入输出格式

输入格式:

输入第一行包含一个整数 NN,表示等腰直角三角形数量。

接下来 NN行,每行三个整数 X, Y, RX,Y,R,分别表示等腰直角三角形的顶点坐标 (X, Y)(X,Y)与腰长 RR。

输入输出样例

输入样例#1: 复制

3
1 1 2
7 1 6
5 3 4

输出样例#1:

24.0

这是自己做出的第一道容斥题(除了一些SB容斥),虽然这道题也不算太难,而且我做了一个晚上。总之就是自己在容斥上还是太菜了。

还是来说题吧。首先要会求多个三角形的交。显然这道题中两个等腰直角的交还是一个等腰直角三角形。我的做法是分类讨论,不如洛谷上题解那么简洁,于是就不说了。

重点是算出容斥系数。我们设个三角形的交的容斥系数为。显然。对于1" class="mathcode" src="https://private.codecogs.com/gif.latex?i%3E1">,的情况,我们先设初值,显然如果为奇数,那么初值为1,否则为0。然后我们要容斥去重。。计算个三角形的交的时候,个三角形的交会被计算次。

然后通过观察证明可以知道

代码:

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstring>
#include<cmath>
#include<queue>
#include<set>
#include<map>
#include<vector>
#include<ctime>
#include<iomanip>
#define ll long long
#define N 11 using namespace std;
inline int Get() {int x=0,f=1;char ch=getchar();while(ch<'0'||ch>'9') {if(ch=='-') f=-1;ch=getchar();}while('0'<=ch&&ch<='9') {x=(x<<1)+(x<<3)+ch-'0';ch=getchar();}return x*f;} int n;
ll ans;
ll f[N];
ll c[N][N];
struct tri {ll x,y,r;}s[N],tem,g;
void work(tri a,tri b) {
if(a.y>b.y) swap(a,b);
if(a.x<=b.x) {
int r=min(b.r,a.r-(b.x+b.y-a.x-a.y));
if(r<0) return g.r=-1,void();
g=b;
g.r=r;
} else {
int r=min(a.y+a.r-b.y,b.x+b.r-a.x);
if(r<0) return g.r=-1,void();
g.x=a.x,g.y=b.y;
g.r=r;
}
}
void dfs(int v,ll x,ll y,ll r,int tot) {
if(r<=0&&tot) return ;
if(v>n) {
if(!tot) return ;
ans+=r*r*f[tot];
return ;
}
dfs(v+1,x,y,r,tot);
if(!tot) {
dfs(v+1,s[v].x,s[v].y,s[v].r,tot+1);
} else {
tem.x=x,tem.y=y,tem.r=r;
work(s[v],tem);
dfs(v+1,g.x,g.y,g.r,tot+1);
}
}
int main() {
n=Get();
c[0][0]=1;
for(int i=1;i<=n;i++) {
for(int j=0;j<=i;j++) {
c[i][j]=(!j||i==j)?1:c[i-1][j-1]+c[i-1][j];
}
}
for(int i=1;i<=n;i++) {
f[i]=(i&1)?1:0;
for(int j=1;j<i;j++) {
f[i]-=f[j]*c[i][j];
}
}
//f[i]=(-1)^(i+1)*2^(i-1)
for(int i=1;i<=n;i++) {
s[i].x=Get(),s[i].y=Get(),s[i].r=Get();
}
dfs(1,0,0,0,0);
cout<<fixed<<setprecision(1)<<1.0*ans/2;
return 0;
}

luogu P4515 [COCI2009-2010#6] XOR的更多相关文章

  1. luogu P4515 [COCI2009-2010#6] XOR 容斥

    LINK:XOR 一个不常见的容斥套路题. 以往是只求三角形面积的交 现在需要求被奇数次覆盖的区域的面积. 打住 求三角形面积的交我也不会写 不过这道题的三角形非常特殊 等腰直角 且直角点都在左下方 ...

  2. [Luogu P4180][BJWC 2010]严格次小生成树

    严格次小生成树,关键是“严格”,如果是不严格的其实只需要枚举每条不在最小生成树的边,如果得到边权和大于等于最小生成树的结束就行.原理就是因为Kruskal非常贪心,只要随便改一条边就能得到一个非严格的 ...

  3. luogu P4385 [COCI2009]Dvapravca

    传送门 我真的弱,正解都不会还打了个错的暴力 考虑平行线与x轴平行,那么可以按照y为第一关键字升序,x为第二关键字升序排序,然后合法的一段红点就是连续的一段,答案也就是最大的连续红色段 推广到一般情况 ...

  4. POJ 1703 Find them, Catch them(并查集高级应用)

    手动博客搬家:本文发表于20170805 21:25:49, 原地址https://blog.csdn.net/suncongbo/article/details/76735893 URL: http ...

  5. luogu P2574 XOR的艺术 (线段树)

    luogu P2574 XOR的艺术 (线段树) 算是比较简单的线段树. 当区间修改时.\(1 xor 1 = 0,0 xor 1 = 1\)所以就是区间元素个数减去以前的\(1\)的个数就是现在\( ...

  6. [luogu]P1800 software_NOI导刊2010提高(06)[DP][二分答案]

    [luogu]P1800 software_NOI导刊2010提高(06) 题目描述 一个软件开发公司同时要开发两个软件,并且要同时交付给用户,现在公司为了尽快完成这一任务,将每个软件划分成m个模块, ...

  7. Luogu P1801 黑匣子_NOI导刊2010提高(06)

    P1801 黑匣子_NOI导刊2010提高(06) 题目描述 Black Box是一种原始的数据库.它可以储存一个整数数组,还有一个特别的变量i.最开始的时候Black Box是空的.而i等于0.这个 ...

  8. 【luogu P1774 最接近神的人_NOI导刊2010提高(02)】 题解

    题目链接:https://www.luogu.org/problemnew/show/P1774 归并排序求逆序对. #include <cstdio> #define livelove ...

  9. 【luogu P1801 黑匣子_NOI导刊2010提高(06)】 题解

    题目链接:https://www.luogu.org/problemnew/show/P1801 替罪羊树吼啊! #include <cstdio> #include <cstrin ...

随机推荐

  1. [译]ABP vNext介绍

    译者注 ASP.NET Boilerplate是.Net平台非常优秀的一个开源Web应用程序框架,在国内也有大量的粉丝. 近日, 本人在github上闲逛, 发现ASP.NET Boilerplate ...

  2. mysql真的不能做搜索引擎吗?

    大家都对电商的商品查询并不陌生,比如我们想根据商品名称查询所有商品信息. 有些技术的童鞋第一念头是搜索引擎:有些技术的童鞋第一念头是模糊查询,如like?(如果商品信息存放到mysql里,我们一般使用 ...

  3. #8 Python数学方法

    前言 前几节了解了Python的不同数据类型,有小伙伴会问,不同的数据类型之间是否可以相互转换?肯定是可以的,本篇博文主要记录数字类型的转换,其他类型的相互转换会在下几节记录,Here we go! ...

  4. ajax+ashx:实现文件的批量导出

    背景: 最近公司有一个需求,就是实现excle的批量导出(一次性导出多个excle). 实现方式: 想到的实现方式: 1.发起一个导出请求,然后批量生产需要导出的excle文件,最后将文件生成一个压缩 ...

  5. 浏览器能正常访问的url,superagent不能正常访问

    在写音乐播放器的过程中,我需要获取qq音乐排行榜的信息,于是我向以前一样,在后台的MusicController中添加一个getTopList方法 然后写下以下代码 // 获取排行 async get ...

  6. 简单说一下UWP中的JumpList

    在Windows10的10856这个版本中,微软为桌面版提供了一组新的应用交互方式,磁贴和Toast通知的个性化都有了一定的改善.针对磁贴方面,微软为我们提供了一组新的API来扩充我们对应用的交互方式 ...

  7. Linux常用基本命令:三剑客命令之-awk内置函数用法

    awk的内置函数大致可以分类为算数函数.字符串函数.时间函数.其他函数等 算数函数 最常用的算数函数有rand函数.srand函数.int函数. 可以使用rand函数生成随机数,但是使用rand函数时 ...

  8. JavaScript中的递归

    译者按: 程序员应该知道递归,但是你真的知道是怎么回事么? 原文: All About Recursion, PTC, TCO and STC in JavaScript 译者: Fundebug 为 ...

  9. linux学习笔记-lrmi源码包的编译安装方法

    我的邮箱地址:zytrenren@163.com欢迎大家交流学习纠错! 官方的lrmi包没有人更新了,如果碰到需要这个编译安装这个包,可以参考我的解决思路,如下: https://pkgs.org/这 ...

  10. js识别设备

    console.log(window.navigator); Navigator 对象属性         appCodeName     返回浏览器的代码名.         appMinorVer ...