POJ 2352 Stars

Description

Astronomers often examine star maps where stars are represented by points on a plane and each star has Cartesian coordinates. Let the level of a star be an amount of the stars that are not higher and not to the right of the given star. Astronomers want to know the distribution of the levels of the stars. 

For example, look at the map shown on the figure above. Level of the star number 5 is equal to 3 (it's formed by three stars with a numbers 1, 2 and 4). And the levels of the stars numbered by 2 and 4 are 1. At this map there are only one star of the level 0, two stars of the level 1, one star of the level 2, and one star of the level 3. 

You are to write a program that will count the amounts of the stars of each level on a given map.

Input

The first line of the input file contains a number of stars N (1<=N<=15000). The following N lines describe coordinates of stars (two integers X and Y per line separated by a space, 0<=X,Y<=32000). There can be only one star at one point of the plane. Stars are listed in ascending order of Y coordinate. Stars with equal Y coordinates are listed in ascending order of X coordinate. 

Output

The output should contain N lines, one number per line. The first line contains amount of stars of the level 0, the second does amount of stars of the level 1 and so on, the last line contains amount of stars of the level N-1.

Sample Input

5
1 1
5 1
7 1
3 3
5 5

Sample Output

1
2
1
1
0
题目大意: 有n个星星, 现在分别给出它们按y递增的坐标,每个星星有一个等级(该星星的等级是x坐标和y坐标都不大于该星的星星数),先要求出每个等级的星星有多少个
这个是线段树第一次写,详解请关注大神博客:
http://blog.youkuaiyun.com/ulquiorra0cifer/article/details/7769675
LANGUAGE:C
CODE:

  
#include<stdio.h>
#include<string.h>
#define max 33000

struct
{
    int l,r,sum;
}tree[max<<2];

void build(int l,int r,int idx)
{
    int mid=(l+r)>>1;//get the midle
    tree[idx].sum=0;//init sum
    if(l==r)//at same point
    {
        tree[idx].l=tree[idx].r=r;
        return ;
    }
    tree[idx].l=l;tree[idx].r=r;
    build(l,mid,idx<<1);//build left subtree
    build(mid+1,r,idx<<1|1);//build(mid+1,r,idx*2+1);build right subtree
}

void add(int l,int r,int idx)
{
    int mid=(tree[idx].l+tree[idx].r)>>1;
    if(tree[idx].l==l&&tree[idx].r==r)//if find the segment
    {
        tree[idx].sum++;//current segement
        return ;
    }
    if(l>mid)
	    add(l,r,idx<<1|1);//if the segement at root's right:add to the right subtree
    else if(r<=mid)add(l,r,idx<<1);//if the segement at root's left:add to the left subtree
    else
    {
        add(l,mid,idx<<1);//add to the left subtree
        add(mid+1,r,idx<<1|1);//add to the right subtree
    }
    tree[idx].sum=tree[idx<<1|1].sum+tree[idx<<1].sum;//get sum
}

int query(int l,int r,int idx)
{
    int mid=(tree[idx].l+tree[idx].r)>>1;
    if(tree[idx].l==l&&tree[idx].r==r)//if find the segement,return sum;
        return tree[idx].sum;
    if(l>mid)
        return query(l,r,idx<<1|1);//return right subtree's sum;
    else if(r<=mid)
        return query(l,r,idx<<1);//return left subtree's sum;
    else
        return query(l,mid,idx<<1)+query(mid+1,r,idx<<1|1);//return the sum of subtree's sum;
}

int main()
{
    int i,n,x,y,lev[max];
    while(scanf("%d",&n)!=EOF)
    {
        memset(lev,0,sizeof(lev));
        build(0,max,1);
        for(i=0;i<n;i++)
        {
            scanf("%d%d",&x,&y);
            lev[query(0,x,1)]++;
            add(x,x,1);
        }
        for(i=0;i<n;i++)
            printf("%d\n",lev[i]);
    }
    return 0;
}



 
 
基于可靠性评估序贯蒙特卡洛模拟法的配电网可靠性评估研究(Matlab代码实现)内容概要:本文围绕“基于可靠性评估序贯蒙特卡洛模拟法的配电网可靠性评估研究”,介绍了利用Matlab代码实现配电网可靠性的仿真分析方法。重点采用序贯蒙特卡洛模拟法对配电网进行长时间段的状态抽样与统计,通过模拟系统元件的故障与修复过程,评估配电网的关键可靠性指标,如系统停电频率、停电持续时间、负荷点可靠性等。该方法能够有效处理复杂网络结构与设备时序特性,提升评估精度,适用于含分布式电源、电动汽车等新型负荷接入的现代配电网。文中提供了完整的Matlab实现代码与案例分析,便于复现和扩展应用。; 适合人群:具备电力系统基础知识和Matlab编程能力的高校研究生、科研人员及电力行业技术人员,尤其适合从事配电网规划、运行与可靠性分析相关工作的人员; 使用场景及目标:①掌握序贯蒙特卡洛模拟法在电力系统可靠性评估中的基本原理与实现流程;②学习如何通过Matlab构建配电网仿真模型并进行状态转移模拟;③应用于含新能源接入的复杂配电网可靠性定量评估与优化设计; 阅读建议:建议结合文中提供的Matlab代码逐段调试运行,理解状态抽样、故障判断、修复逻辑及指标统计的具体实现方式,同时可扩展至不同网络结构或加入更多不确定性因素进行深化研究。
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值