题目描述：

　　给定 n 个同心的扇形，求有多少面积，被至少k 个扇形所覆盖。

输入输出格式

输入格式：

第一行是三个整数 n，m，k。n 代表同心扇形个数，m代表将(−π ,π ]的角度区间平均分成2m 份。

从第二行开始的 n 行，每行三个整数r，a1，a2。描述了一个圆心在原点的扇形，半径为r，圆心角是从弧度

　　输出格式：

　　　　　　输出一个整数 ans ，数据保证答案在 $2^{63} - 1263−1范围内。$

var
    next_insert,val_insert,next_delete,val_delete:array[1..4000000]of longint;
    head_insert,head_delete:array[-1000000..1000000]of longint;
    cnt:array[-1000000..1000000]of longint;
    a:array[1..2000000]of longint;
    tot,tot_insert,tot_delete,v:longint;
    ans,x,y,n,m,need,s,t,r:int64;
    i,j:longint;
procedure add_insert(x,v:longint);
begin
    inc(tot_insert);
    next_insert[tot_insert]:=head_insert[x];
    head_insert[x]:=tot_insert;
    val_insert[tot_insert]:=v;
end;
procedure add_delete(x,v:longint);
begin
    inc(tot_delete);
    next_delete[tot_delete]:=head_delete[x];
    head_delete[x]:=tot_delete;
    val_delete[tot_delete]:=v;
end;
procedure add(s,t,r:int64);
begin          
    if s=t then exit;
    inc(cnt[s]); add_insert(s,r);
    dec(cnt[t]); add_delete(t,r);
end;
procedure update(k,l,r,x,z:longint);
var
    mid:longint;
begin
    a[k]:=a[k]+z;
    if l=r then exit;
    mid:=(l+r)>>1;
    if x<=mid then update(k*2,l,mid,x,z) else update(k*2+1,mid+1,r,x,z);
end;
function query(k,l,r,need:longint):longint;
var
    mid:longint;
begin
    if l=r then exit(l);
    mid:=(l+r)>>1;
    if a[k*2+1]>=need then exit(query(k*2+1,mid+1,r,need))
    else exit(query(k*2,l,mid,need-a[k*2+1]));
end;
begin
    read(n,m,need);
    for i:=1 to n do
    begin
        read(r,s,t);　　
　　　　　　　　//是不是觉得这个读入处理很恶心，嗯，我也这么觉得。。。。。。（当我看到a1不一定小于a2时，我真想提着西瓜刀去找出题人拼命）
        if (s>=0)and(t>=0) then
        begin
            if s<t then add(s,t,r)
            else
            begin
                add(s,m,r);
                add(-m,0,r);
                add(0,t,r);
            end;
        end else
        if (s<0)and(t<0) then
        begin
            if s<t then add(s,t,r)
            else
            begin
                add(-m,t,r);
                add(0,m,r);
                add(s,0,r);
            end;
        end else
        if (s>=0)and(t<0) then
        begin
            add(s,m,r);
            add(-m,t,r);
        end else
        if (s<0)and(t>=0) then
        begin
            add(0,t,r);
            add(s,0,r);
        end;
    end;
    for i:=-m to m-1 do
    begin
        j:=head_insert[i];
        while j>0 do
        begin
            v:=val_insert[j];
            update(1,0,100000,v,1);
            j:=next_insert[j];
        end;
        j:=head_delete[i];
        while j>0 do
        begin
            v:=val_delete[j];
            update(1,0,100000,v,-1);
            j:=next_delete[j];
        end;
        x:=x+cnt[i];
        if x>=need then
        begin
            y:=query(1,0,100000,need);
            ans:=ans+sqr(y);
        end;
    end;
    writeln(ans);
end.

SHOI 2013 【扇形面积并】

题目描述：

思路分析：

代码实现：