HOJ 1867 经理的烦恼 树状数组

/*

题目:

   给出n个商店,每个商店有初始化的商品数量,现在有两种指令:

   0 x y 连锁店x的商品数量变化值为y,y > 0商品数量增加, y < 0减少

   1 i j 输出编号在[i,j]区间内的连锁店中商品数量为素数的有多少家

   当指令为1时,求给出 i到j的商品数量为素数的商店的数目

 

分析:

   用树状数组做,只不过在update()处改动一下,并且注意到有可能改动完后和改动前

   均为素数的情况,所以用visit[]数组记录此前是否为素数,另外需要注意的是当初始

   化时即为奇数的话,当时就要初始化数组c[]以及visit[]数组

*/

#include <iostream>

#include <cstdio>

#include <cstring>

using namespace std;

const int X = 1000010;

int a[X],c[X],n,m,k;

bool visit[X];

bool check(int x)  //判断是否是素数

{

   if(x<=1)    //小于等于一直接返回false

      return false;

   for(int i=2;i*i<=x;i++)

      if(x%i==0)

         return false;

   return true;

}

int lowbit(int x)

{

   return x & -x;

}

int sum(int x)        //求和

{

   int ans = 0;

   while(x>0)

   {

      ans += c[x];

      x -= lowbit(x);

   }

   return ans;

}

void update(int x)

{

   if(check(a[x]))       //如果这个商店的商品数为素数

   {

      if(!visit[x])      //如果之前不是素数的话,现在为素数

      {

         visit[x] = true;

         while(x<=n)

         {

            c[x]++;

            x += lowbit(x);

         }

      }

   }

   else               //如果这个商店的商品数为偶数

   {

      if(visit[x])    //如果之前是素数,但现在不是素数

      {

         visit[x] = false;

         while(x<=n)

         {

            c[x]--;

            x += lowbit(x);

         }

      }

   }

}

int main()

{

   freopen("sum.in","r",stdin);

   freopen("sum.out","w",stdout);

   int x,y,z,cnt = 0;

   while(scanf("%d%d%d",&n,&m,&k),n||m||k)

   {

      printf("CASE #%d:\n",++cnt);

      for(int i=1;i<=n;i++)    //初始化

         a[i] = k;

      if(check(k))          //如果刚开始即为素数

      {

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

         {

            c[i] = lowbit(i);

            visit[i] = true;   //visit数组为true时表示之前的数为素数

         }

      }

      else

      {

         memset(c,0,sizeof(c));

         memset(visit,false,sizeof(visit));

      }

      for(int i=1;i<=m;i++)

      {

         scanf("%d%d%d",&z,&x,&y);

         if(z)

            printf("%d\n",sum(y)-sum(x-1));

         //计算的是[i,j]的区间内的连锁店中商品数量为素数的有多少家

 

         else

         {

            a[x] += y;         //x商店的变化

            update(x);         //更新x商店的变化给其他商店带来的变化

         }

      }

      cout<<endl;

   }

   return 0;

}

 

转载于:https://www.cnblogs.com/yejinru/archive/2012/03/22/2411671.html

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

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值