[spoj CIRU2]关于圆面积并的思考，非simpson！！

最新推荐文章于 2022-04-12 19:07:58 发布
最新推荐文章于 2022-04-12 19:07:58 发布
阅读量512 收藏
点赞数 1
文章标签： 算法
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/u013593293/article/details/54934113
版权

aekdycoin神犇的讲解
picture1
图一是普通求圆并的思路，复杂度O（n^2 logn）
picture2
图二是神犇给的提示，没有题解…

    一开始困扰了好久，我觉得这和hzwer开的坏头有关，他在BZOJ上圆并用了simpson之后一大批同学思路都跟着“奇怪”了起来，这题simpson没试过，不过评论说会挂精度。
    首先我们要理解非simpson的圆并在干什么？这个算法很神的地方在于他把圆并分成了好多个圆弧和一个多边形，然后分别求面积，圆弧部分是一个圆上未被其它圆覆盖的部分的圆弧（可能有多段），而多边形则被分解为每条边单独算贡献。这样我们的工作就简单许多了，枚举一个圆，与其它圆求交点，搞出区间后排序，统计未覆盖的角度算圆弧，统计坐标算叉积求贡献。最后统统加起来就好了。
    注意一些被包含的圆没有作用，标记掉即可！
    现在再看那张提示图，我们发现一个圆未被其它圆覆盖的部分就是圆并的“一次覆盖面积”（见spoj ciru2），我们大胆的猜一个结论：“k次覆盖面积”就是k次覆盖的圆弧加多边形面积。呵呵...结果就是这样...
    注意这时被包含的圆不是去掉，是统计时覆盖次数加一！
    废话就这么多，总之细节很烦，对拍造数据要考虑数据的合法性...

普通圆并

#include <cmath>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <iostream>
#include <algorithm>
using namespace std;
typedef long long LL;
inline int read() {
    int x=0,f=1; char ch=getchar();
    while (ch<'0'||ch>'9') { if (ch=='-') f=-1; ch=getchar(); }
    while (ch>='0'&&ch<='9') { x=x*10+ch-'0'; ch=getchar(); }
    return x*f;
}

const int MAXN=1003;
const double PI=3.141592653589793238462,eps=1e-8;
int n,m;
double ans1,ans2,b[MAXN],R[MAXN];
struct P{
    double x,y;
    P(double v=0,double u=0){ x=v; y=u; }
    double len(){ return sqrt(x*x+y*y); }
    P operator +(const P &E)const{ return P(x+E.x,y+E.y); }
    P operator -(const P &E)const{ return P(x-E.x,y-E.y); }
    P operator *(const double &z)const{ return P(x*z,y*z); }
    P operator /(const double &z)const{ return P(x/z,y/z); }
    double operator *(const P &E)const{ return x*E.x+y*E.y; }
}a[MAXN];
inline double cross(P p,P q){ return p.x*q.y-p.y*q.x; }
inline double dist(P p,P q){ return sqrt((p.x-q.x)*(p.x-q.x)+(p.y-q.y)*(p.y-q.y)); }
struct Data{
    P ax,ay;
    double x,y;
    bool operator <(const Data &E)const{ return (x<E.x)||((x==E.x)&&(y<E.y)); }
}c[MAXN*2];
inline P rot(P p,double a){ return P(cos(a)*p.x-sin(a)*p.y,sin(a)*p.x+cos(a)*p.y); }
inline double deg(P p){
    double tmp=atan2(p.y,p.x);
    return tmp<0?tmp+PI*2:tmp;
}
inline void work(P o1,P o2,double r1,double r2){
    double s,p,h,alpha,tmp=dist(o1,o2);
    p=(r1+r2+tmp)/2; s=sqrt(p*(p-r1)*(p-r2)*(p-tmp));
    h=s*2/tmp; alpha=asin(h/r1);
    if (r1*r1+tmp*tmp<r2*r2) alpha=PI-alpha;
    m++;
    c[m].ax=rot(o2-o1,PI*2-alpha)/tmp*r1+o1;
    c[m].ay=rot(o2-o1,alpha)/tmp*r1+o1;
    c[m].x=deg(c[m].ax-o1);
    c[m].y=deg(c[m].ay-o1);
    if (c[m].x>c[m].y) {
        m++; c[m].y=c[m-1].y; c[m].ay=c[m-1].ay;
        c[m-1].y=PI*2; c[m].x=0;
        c[m-1].ay=c[m].ax=o1+P(r1,0);
    }
}
int main() {
    //freopen("1.in","r",stdin);
    //freopen("2.out","w",stdout);
    int i,j,k,l;
    m=read();
    for (i=1,n=0;i<=m;++i) {
        n++; a[n].x=read(); a[n].y=read(); R[n]=read();
        if (fabs(R[n])<=eps) n--;
    }
    memset(b,0,sizeof(b));
    for (i=1;i<=n;++i)
    if (!b[i])
        for (j=1;j<=n;++j)
        if (!b[j]&&i!=j&&dist(a[i],a[j])<=R[j]-R[i]) { b[i]=1; break; }
    //for (i=1;i<=n;++i) if (!b[i]) printf("%d\n",i);
    for (i=1,ans1=ans2=0;i<=n;++i)
    if (!b[i]) {
        for (j=1,m=0;j<=n;++j)
        if (!b[j]&&i!=j&&dist(a[i],a[j])<R[i]+R[j]) work(a[i],a[j],R[i],R[j]);
        if (m==0) {
            ans2+=PI*R[i]*R[i];
            continue;
        }
        sort(c+1,c+1+m); c[0].y=0; c[m+1].x=PI*2; c[0].ay=c[m+1].ax=a[i]+P(R[i],0);
        /*for (j=1;j<=m;++j) {
            printf("%.8lf %.8lf\n",c[j].x,c[j].y);
            printf("%.8lf %.8lf||%.8lf %.8lf\n",c[j].ax.x,c[j].ax.y,c[j].ay.x,c[j].ay.y);
        }*/
        double tmp=0;
        for (j=0;j<=m;j=k+1) {
            for (k=l=j;k<m&&c[k+1].x<=c[l].y;++k) if (c[k+1].y>c[l].y) l=k+1;
            ans1+=cross(c[l].ay,c[k+1].ax);
            tmp=c[k+1].x-c[l].y;
            ans2+=R[i]*R[i]*(tmp-sin(tmp))/2;
        }
        //printf("%d %d %.8lf %.8lf\n",i,m,ans1,ans2);
    }
    printf("%.5lf\n",fabs(ans1)/2+ans2);

    //fclose(stdin);
    //fclose(stdout);
    return 0;
}

/*
2
4 2 5
3 3 4

*/

spoj ciru2

#include <cmath>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <iostream>
#include <algorithm>
using namespace std;
typedef long long LL;
inline int read() {
    int x=0,f=1; char ch=getchar();
    while (ch<'0'||ch>'9') { if (ch=='-') f=-1; ch=getchar(); }
    while (ch>='0'&&ch<='9') { x=x*10+ch-'0'; ch=getchar(); }
    return x*f;
}

const int MAXN=1003;
const double PI=3.1415926535897932384626433832795028841971693993,eps=1e-8;
int n,m;
double ans[MAXN][2],b[MAXN],R[MAXN];
struct P{
    double x,y;
    P(double v=0,double u=0){ x=v; y=u; }
    double len(){ return sqrt(x*x+y*y); }
    P operator +(const P &E)const{ return P(x+E.x,y+E.y); }
    P operator -(const P &E)const{ return P(x-E.x,y-E.y); }
    P operator *(const double &z)const{ return P(x*z,y*z); }
    P operator /(const double &z)const{ return P(x/z,y/z); }
    double operator *(const P &E)const{ return x*E.x+y*E.y; }
}a[MAXN];
inline double cross(P p,P q){ return p.x*q.y-p.y*q.x; }
inline double dist(P p,P q){ return sqrt((p.x-q.x)*(p.x-q.x)+(p.y-q.y)*(p.y-q.y)); }
struct Data{
    P ax; double x; int y;
    bool operator <(const Data &E)const{ return x<E.x; }
}c[MAXN*10];
inline P rot(P p,double a){ return P(cos(a)*p.x-sin(a)*p.y,sin(a)*p.x+cos(a)*p.y); }
inline double deg(P p){
    double tmp=atan2(p.y,p.x);
    return tmp<0?tmp+PI*2:tmp;
}
inline void work(P o1,P o2,double r1,double r2){
    double s,p,h,alpha,tmp=dist(o1,o2);
    p=(r1+r2+tmp)/2; s=sqrt(p*(p-r1)*(p-r2)*(p-tmp));
    h=s*2/tmp; alpha=asin(h/r1);
    if (r1*r1+tmp*tmp<r2*r2) alpha=PI-alpha;
    c[++m].ax=rot(o2-o1,PI*2-alpha)/tmp*r1+o1; c[m].x=deg(c[m].ax-o1); c[m].y=1;
    c[++m].ax=rot(o2-o1,alpha)/tmp*r1+o1; c[m].x=deg(c[m].ax-o1); c[m].y=-1;
    if (c[m-1].x>c[m].x) {
        c[++m].ax=o1+P(r1,0); c[m].x=0; c[m].y=1;
        c[++m].ax=o1+P(r1,0); c[m].x=PI*2; c[m].y=-1;
    }
}
int main() {
    //freopen("1.in","r",stdin);
    //freopen("2.out","w",stdout);

    int i,j,k,l;
    n=read();
    for (i=1;i<=n;++i) { a[i].x=read(); a[i].y=read(); R[i]=read(); }
    memset(b,0,sizeof(b));
    for (i=1;i<=n;++i)
        for (j=1;j<=n;++j)
        if (i!=j&&dist(a[i],a[j])<=R[j]-R[i]) b[i]++;
    memset(ans,0,sizeof(ans));
    for (i=1;i<=n;++i) {
        for (j=1,m=0;j<=n;++j)
        if (i!=j&&dist(a[i],a[j])<R[i]+R[j]&&dist(a[i],a[j])>fabs(R[j]-R[i]))
            work(a[i],a[j],R[i],R[j]);
        c[++m].ax=a[i]+P(R[i],0); c[m].x=0; c[m].y=b[i]+1;
        c[++m].ax=a[i]+P(R[i],0); c[m].x=PI*2; c[m].y=-b[i]-1;
        sort(c+1,c+1+m);
        /*printf("%d\n",i);
        for (j=1;j<=m;++j) {
            printf("%d %.3lf %d\n",j,c[j].x,c[j].y);
            printf("(%.3lf,%.3lf)\n",c[j].ax.x,c[j].ax.y);
        }*/
        int sa=1,sb=0,sc;
        for (j=1;j<=m;j=k+1) {
            for (k=j;k<m&&fabs(c[k+1].x-c[j].x)<eps;++k);
            for (l=j,sc=sb;l<=k;++l) sc+=c[l].y;
            if (sc==sb&&k!=m) continue;
            if (fabs(c[sa].x-c[j].x)>eps) {
                ans[sb][0]+=cross(c[sa].ax,c[j].ax);
                double tmp=c[j].x-c[sa].x;
                ans[sb][1]+=R[i]*R[i]*(tmp-sin(tmp))/2;
            }
            sb=sc; sa=j;
        }
    }
    //for (i=1;i<=n;++i) printf("%.3lf %.3lf %.3lf\n",ans[i][0],ans[i][1],fabs(ans[i][0])/2+ans[i][1]);
    for (i=1;i<=n;++i) ans[i][0]=fabs(ans[i][0])/2+ans[i][1];
    for (i=1;i<=n;++i) printf("[%d] = %.3lf\n",i,ans[i][0]-ans[i+1][0]);

    //fclose(stdin);
    //fclose(stdout);
    return 0;
}

/*
5
3 0 3
4 3 5
1 2 4
4 3 5
1 4 4

*/

    身为一个蒟蒻非常同情蒟蒻写完程序死WA不止，希望我的程序对蒟蒻们的学习(对拍)有帮助。
llllxxx666888
关注
  • 1
    点赞
  • 踩
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
SPOJ CIRUT(CIRU2-圆的面积并与恰被覆盖k次区域的面积)
nike0good |Oier&ACMer | 熟能生巧
08-11 1643
给n<=1000个圆，对于所有k，求恰被覆盖k次区域的面积。考虑积分。 根据格林公式，对平面的积分可以转换成对边界的积分。 被覆盖k次的边界恰为答案。如何判断每个边界的覆盖次数？一条边界被覆盖的次数=所在圆被覆盖的次数+与该圆相交的圆对圆进行切割后每一段隶属于几个圆里#include<bits/stdc++.h> using namespace std; #define For(i,n) fo
CIRU2（K次圆交板子）
twh233的算法blog
08-12 1198
传送门 圆的面积并 题意：给n个圆（最多1000个），分别求出覆盖1层的面积、覆盖2层的面积、 覆盖3层的面积。。。覆盖n层的面积 //china no.1 #pragma comment(linker, "/STACK:1024000000,1024000000") #include #include #include #include #include #include
SPOJ CIRU SPOJ VCIRCLE 圆的面积并问题
zjy2015302395的博客
06-09 1058
给出若干圆 求面积并 找到所有边界，然后格林公式，就是里面的面积。tip:正方向为逆时针,对于所给的N个圆，我们可以进行去冗杂，表现为: (1) 去掉被包含(内含 or 内切)的小圆 () (2) 去掉相同的圆枚举每个圆，并对于剩下的所有圆和它求交点，对于所求的的交点，可以得到一个角度区间 [A,B], 当然区间如果跨越了(例如 [1.5PI, 0.5PI],注意这里是有方向的) 2PI那么需
[SPOJ8119] CIRUT [自适应Simpson/格林公式][计算几何]
最后一次修改
04-08 326
Link Luogu - https://www.luogu.org/problemnew/show/SP8119 用几何方法瞎搞即可 也可以直接套格林公式 格林公式的具体意义，只要理解了围道积分的定义然后在纸上画一画一个简单情况就懂啦 简单地讲就是，对于我们需要求面积的平面单联通图形，它的轮廓线上的所有点可以被分为两部分 一部分是“面向原点”的，一部分是“背对原点”的， 面向原点的那一部分组成...
格林公式求圆并的面积及重心
RBS的专栏
04-12 4319
ICPC格林公式计算圆并的面积与重心
SPOJ-2.zip_SPOJ_SPOJ-TRIKA.cpp_online judge_sPOJ-Solution_spoj2
09-24
4. **SPOJ-SUMSUMS.cpp**：从名字推测，可能是关于求和或数组处理的题目。这类问题常常需要掌握快速求和技巧，如前缀和、 Fenwick树或Segment Tree等高级数据结构。 5. **SPOJ-TAILS.cpp**：可能与概率统计或数列...
SPOJ.rar_SPOJ
09-19
SPOJ，全称Sphere Online Judge，是一个全球知名的在线编程竞赛平台，旨在提供一个环境让程序员可以测试、提交并评估他们的算法解决方案。这个名为"SPOJ.rar_SPOJ"的压缩包文件，很可能包含了越南版SPOJ（VN.SPOJ）...
SPOJ.zip_SPOJ_sgu_them
09-24
2. SPOJ-MAXSUMSQ.cpp：此代码可能解决的是在一个数组中找到最大平方子数组和的问题，动态规划可以通过维护一个二维表来记录每个子数组的最大和，从而找到最大的平方子数组。 3. SPOJ-MBLAST.cpp：可能涉及的是寻找...
SPOJ-Solutions:SPOJ算法问题的解决方案
03-13
SPOJ（Sphere Online Judge）是一个在线编程竞赛平台，它提供了大量的算法问题供程序员解决，以提高他们的编程技能和算法理解。"SPOJ Solutions"是这个领域的一个资源集合，包含了解决SPOJ上各种问题的代码示例，...
RandomGoCode:算法，SPOJ，涂鸦...但是在Go中！
03-21
2. **算法应用** RandomGoCode项目涵盖了各种算法，如排序（快速排序、归并排序、冒泡排序等）、搜索（二分查找、深度优先搜索、广度优先搜索等）、图论（Dijkstra算法、Floyd-Warshall算法等）、动态规划等。通过...
【BZOJ2178】圆的面积并
cz_xuyixuan的博客
03-27 1183
【题目链接】点击打开链接【思路要点】将X轴根据整点分成若干段，对每一段横坐标用自适应Simpson积分计算面积。时间复杂度\(O(Cnt*NLogN)\)，其中\(Cnt\)为自适应Simpson积分的迭代次数。【代码】#include&lt;bits/stdc++.h&gt; using namespace std; const int MAXN = 1005; const int INF = ...
BNU 4067 求圆并
weixin_30568715的博客
07-11 94
好久没写过单组数据的题目了 QAQ 赤裸裸的模板题 #include <cstdio> #include <cstring> #include <iostream> #include <algorithm> #include <cmath> using namespace std; #define sqr(x) ((x...
圆并模板
zkxwhanybz828的专栏
03-08 392
#include #include #include #include using namespace std; const int N = 1100; const double eps = 1e-8, PI = 3.14159265358979; int n;bool v[N]; inline double sqr (double x) {return x * x; } inline i
SPOJ 8073 The area of the union of circles（计算几何の圆并）（CIRU）
weixin_30553837的博客
11-15 216
Description You are given N circles and expected to calculate the area of the union of the circles ! Input The first line is one integer n indicates the number of the circles. (1 <= n <= 10...
SPOJ CIRU The area of the union of circles
dezhen7015的博客
01-26 93
You are given N circles and expected to calculate the area of the union of the circles ! Input The first line is one integer n indicates the number of the circles. (1 <= n <= 1000) Then ...
判断矩形相交及求面积的交并集
Handsomk
06-20 874
判断矩形相交及求面积的交并集 定义一个平面矩形类，从键盘输入两个矩形，判断这两个矩形是否相交？计算两矩形并集、交集的面积。 #include <iostream> using namespace std; class rectangle{ private: float x1, y1, x2, y2; public: rectangle(){ x...
圆面积并积
yangactive的专栏
04-23 237
http://hi.baidu.com/aekdycoin/blog/item/c1b28e3711246b3f0b55a95e.html  
【SPOJ】AIBOHP - Aibohphobia
Blog
05-03 585
Problem hereSolution答案＝原字串長度-原字串與其反轉後字串的lcm一開始用top down寫，TLE orz…. 之後改成 bottom up 才AC#include <iostream> #include <string> using namespace std; string input, res; int lcs[6101][6101];int main(){
SPOJ - PT07Z （DFS）
CoAAColA的博客
11-22 423
传送门  题面： You are given an unweighted, undirected tree. Write a program to output the length of the longest path (from one node to another) in that tree. The length of a path in this case is number...
洛谷SPOJ如何注册
最新发布
07-27
洛谷的SPOJ需要注册一个SPOJ账号并进行绑定才能进行交题。您可以按照以下步骤进行注册： 1. 打开洛谷网站（https://www.luogu.com.cn/）并登录您的洛谷账号。 2. 在网站顶部导航栏中找到“题库”选项，将鼠标悬停在上面，然后选择“SPOJ”。 3. 在SPOJ页面上，您会看到一个提示，要求您注册SPOJ账号并进行绑定。点击提示中的链接，将会跳转到SPOJ注册页面。 4. 在SPOJ注册页面上，按照要求填写您的用户名、密码和邮箱等信息，并完成注册。 5. 注册完成后，返回洛谷网站，再次进入SPOJ页面。您会看到一个输入框，要求您输入刚刚注册的SPOJ用户名。输入用户名后，点击“绑定”按钮即可完成绑定。 现在您已经成功注册并绑定了SPOJ账号，可以开始在洛谷的SPOJ题库上刷题了。祝您顺利完成编程练习！\[1\]\[2\] #### 引用[.reference_title] - *1* *3* [（洛谷入门系列，适合洛谷新用户）洛谷功能全解](https://blog.csdn.net/rrc12345/article/details/122500057)[target="_blank" data-report-click={"spm":"1018.2226.3001.9630","extra":{"utm_source":"vip_chatgpt_common_search_pc_result","utm_medium":"distribute.pc_search_result.none-task-cask-2~all~insert_cask~default-1-null.142^v91^control_2,239^v3^insert_chatgpt"}} ] [.reference_item] - *2* [luogu p7492 序列](https://blog.csdn.net/zhu_yin233/article/details/122051384)[target="_blank" data-report-click={"spm":"1018.2226.3001.9630","extra":{"utm_source":"vip_chatgpt_common_search_pc_result","utm_medium":"distribute.pc_search_result.none-task-cask-2~all~insert_cask~default-1-null.142^v91^control_2,239^v3^insert_chatgpt"}} ] [.reference_item] [ .reference_list ]

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值