Cod sursa(job #522097)

Utilizator BlaugranasEnal Gemaledin Blaugranas Data 14 ianuarie 2011 12:08:01
Problema Arbori indexati binar Scor 10
Compilator cpp Status done
Runda Arhiva educationala Marime 1.88 kb
#include<stdio.h>
#include<math.h>
#include<iostream.h>
#include<fstream.h>
#define N 100001
long n,m,y,z,i,j,step,a[N];
int x,poz;
unsigned long long c[N],b[N],p;
       
void modificare(unsigned long long c[N],long y,long z,long n)
{long i=y;
int poz=0;
while(i<=n)
      {c[i]+=z;
      while((i&(long)pow(2,poz))==0)
             poz++;
      i+=(long)pow(2,poz);
      poz++;}}
       
long pozitie(unsigned long long c[N],unsigned long long p,long n)
{long step,i=0,j,k=p,r=p;
for(step=1;step<n;step<<=1);
step>>=1;
while(step!=0&&i<n)
       {j=i+step;
       if(k==c[j])
              return j;
       else
              if(k>c[j])
                      {i=j;
                      k-=c[j];}
       step>>=1;}
if(k!=0||step==0)
       return -1;
else
       return i;}

unsigned long long suma(unsigned long long c[N],long y,long z,long n)
{long j=y-1,i=z;
int poz=0;
unsigned long long s1=0,s2=0;
while(i>0)
       {s1+=c[i];
       while((i&(long)pow(2,poz))==0)
             poz++;
       i-=(long)pow(2,poz);
       poz++;}
poz=0;
while(j>0)
       {s2+=c[j];
       while((j&(long)pow(2,poz))==0)
             poz++;
       j-=(long)pow(2,poz);
       poz++;}
return s1-s2;}

int main()
{fstream f1("aib.in",ios::in);
fstream f2("aib.out",ios::out);
f1>>n>>m;
b[0]=0;
for(i=1;i<=n;i++)
       {f1>>a[i];
       b[i]=b[i-1]+a[i];
       j=i;
       poz=0;
       while((j&(long)pow(2,poz))==0)
              poz++;
       j=j-(long)pow(2,poz);
       c[i]=b[i]-b[j];}
for(i=1;i<=m;i++)
       {f1>>x;
       if(x==2)
              {f1>>p;
              f2<<pozitie(c,p,n)<<endl;}
       else
              if(x==0)
                      {f1>>y>>z;
                      modificare(c,y,z,n);}
              else
                      {f1>>y>>z;
                      f2<<suma(c,y,z,n)<<endl;}}
f1.close();
f2.close();
return 0;}