有n朵花,每朵花有三个属性:花形(s)、颜色(c)、气味(m),用三个整数表示。
现在要对每朵花评级,一朵花的级别是它拥有的美丽能超过的花的数量。
定义一朵花A比另一朵花B要美丽,当且仅Sa>=Sb,Ca>=Cb,Ma>=Mb。显然,两朵花可能有同样的属性。需要统计出评出每个等级的花的数量。

维护一个W*W的矩阵，初始值均为S.每次操作可以增加某格子的权值,或询问某子矩阵的总权值.修改操作数M<=160000,询问数Q<=10000,W<=2000000.

第一行两个整数,S,W;其中S为矩阵初始值;W为矩阵大小
接下来每行为一下三种输入之一(不包含引号):
"1 x y a"
"2 x1 y1 x2 y2"
"3"
输入1:你需要把(x,y)(第x行第y列)的格子权值增加a
输入2:你需要求出以左下角为(x1,y1),右上角为(x2,y2)的矩阵内所有格子的权值和,并输出
输入3:表示输入结束

某国为了防御敌国的导弹袭击，发展出一种导弹拦截系统。但是这种导弹拦截系统有一个缺陷：虽然它的第一发炮弹能够到达任意的高度、并且能够拦截任意速度的导弹，但是以后每一发炮弹都不能高于前一发的高度，其拦截的导弹的飞行速度也不能大于前一发。某天，雷达捕捉到敌国的导弹来袭。由于该系统还在试用阶段，所以只有一套系统，因此有可能不能拦截所有的导弹。

在不能拦截所有的导弹的情况下，我们当然要选择使国家损失最小、也就是拦截导弹的数量最多的方案。但是拦截导弹数量的最多的方案有可能有多个，如果有多个最优方案，那么我们会随机选取一个作为最终的拦截导弹行动蓝图。

我方间谍已经获取了所有敌军导弹的高度和速度，你的任务是计算出在执行上述决策时，每枚导弹被拦截掉的概率。

第一行包含一个正整数

，分别表示第 枚导弹的高度和速度。

输出包含两行。

第一行为一个正整数，表示最多能拦截掉的导弹数量；

第二行包含n个0到1之间的实数，第i个数字表示第i枚导弹被拦截掉的概率（你可以保留任意多位有效数字）。

 #include<algorithm>

 #include<iostream>

 #include<cstring>

 #include<cstdio>

 #include<cmath>

 #include<vector>

 using namespace std;

 const int maxn=;

 inline int read()

 {

     int x=,f=;char cc;cc=getchar();

     while(cc>''||cc<''){if(cc=='-')f=-;cc=getchar();}

     while(cc>=''&&cc<=''){x=(x<<)+(x<<)+(cc^);cc=getchar();}

     return x;

 }

 struct tree

 {

     int f;double w;

     tree(){f=,w=;}

 }t[maxn];

 int n;

 int st[maxn],top=,th,tv;

 inline int lowbit(int x){return x&(-x);}

 inline void add(int p,int f,double w)

 {

     while(p<n)

     {

         if(t[p].f<f)

         {

             if(t[p].f==)st[++top]=p;

             t[p].f=f;t[p].w=w;

         }

         else if(t[p].f==f)t[p].w+=w;

         p+=lowbit(p);

     }

     return;

 }

 tree ask(int p)

 {

     tree res;

     while(p)

     {

         if(t[p].f>res.f) res=t[p];

         else if(t[p].f==res.f) res.w+=t[p].w;

         p-=lowbit(p);

     }

     return res;

 }

 struct Dan

 {

     int h,v,f[],id,t;

     double g[];

 }a[maxn],q[maxn];

 int wh[maxn],wv[maxn],id[maxn];

 int rk[maxn];

 inline bool cmp(int x,int y){return a[x].h<a[y].h||(a[x].h==a[y].h&&a[x].id<a[y].id);}

 int cmpid(Dan a,Dan b){return a.id<b.id;}

 int cnt=;

 void cdq(int l,int r,int mode)

 {

     if(l==r)

     {

         if(a[l].f[mode]<){a[l].f[mode]=;a[l].g[mode]=;}

         return;

     }

     int mid=(l+r)>>;

     memcpy(q+l,a+l,sizeof(Dan)*(r-l+));

     int q1=l,q2=mid+;

     for(int i=l;i<=r;i++)

         if(q[i].t<=mid)a[q1++]=q[i];

         else a[q2++]=q[i];

     cdq(l,mid,mode);

     q1=l;

     for(int i=mid+;i<=r;i++)

     {

         while(q1<=mid && a[q1].id<a[i].id)

             add(a[q1].v,a[q1].f[mode],a[q1].g[mode]),q1++;

         tree res=ask(a[i].v);

         if(!res.f)continue;

         if(res.f+>a[i].f[mode])

         {

             a[i].f[mode]=res.f+;

             a[i].g[mode]=res.w;

         }

         else if(res.f+==a[i].f[mode]) a[i].g[mode]+=res.w;

     }

     while(top){t[st[top]].w=;t[st[top--]].f=;}

     cdq(mid+,r,mode);

     merge(a+l,a+mid+,a+mid+,a+r+,q+l,cmpid);

     memcpy(a+l,q+l,sizeof(Dan)*(r-l+));

     return;

 }

 int main()

 {

     //freopen("cnm.txt","r",stdin);

     n=read();

     for(int i=;i<=n;i++)

     {

         a[i].h=read();a[i].v=read();a[i].id=i;

         wh[i]=a[i].h;wv[i]=a[i].v;

         rk[i]=i;

     }

     sort(wh+,wh+n+);

     sort(wv+,wv+n+);

     th=unique(wh+,wh+n+)-wh-;

     tv=unique(wv+,wv+n+)-wv-;

     for(int i=;i<=n;i++)

     {

         a[i].h=th-(lower_bound(wh+,wh+th+,a[i].h)-wh)+;

         a[i].v=tv-(lower_bound(wv+,wv+tv+,a[i].v)-wv)+;

     }

     sort(rk+,rk+n+,cmp);

     for(int i=;i<=n;i++)a[rk[i]].t=i;

     cdq(,n,);

     for(int i=;i<=n;i++)

     {

         a[i].h=th-a[i].h+;

         a[i].v=tv-a[i].v+;

         a[i].id=n-a[i].id+;

         a[i].t=n-a[i].t+;

     }

     reverse(a+,a+n+);

     cdq(,n,);

     reverse(a+,a+n+);

     double smm=;

     int ans=;

     for(register int i=;i<=n;i++)

         ans=max(ans,a[i].f[]+a[i].f[]-);

     printf("%d\n",ans);

     for(int i=;i<=n;i++)

         if(a[i].f[]==ans)

             smm+=a[i].g[]*a[i].g[]*1ll;

     for(int i=;i<=n;i++)

     {

         double res=a[i].g[]*a[i].g[];

         if(a[i].f[]+a[i].f[]-!=ans)printf("%.5lf ",0.0);

         else printf("%.5f ",res/smm);

     }

     return ;

 }