Cod sursa(job #3243870)

Utilizator maryyMaria Ciutea maryy Data 21 septembrie 2024 21:24:47
Problema Arbori indexati binar Scor 100
Compilator cpp-64 Status done
Runda Arhiva educationala Marime 1.62 kb
#include <bits/stdc++.h>

using namespace std;
ifstream in("aib.in");
ofstream out("aib.out");
const int NMAX=1e6;
long long aib[NMAX+1], n, m;
int max_p2;
void fenwick_build()
{
    max_p2 = n;
    while (max_p2 & (max_p2 - 1))
    {
        max_p2 &= (max_p2 - 1);
    }
    for (int i = 1; i <= n; i++)
    {
        int j = i + (i & -i);
        if (j <= n) {
          aib[j] += aib[i];
        }
    }
}
void fenwick_add(int pos, int val)
{
    while(pos<=n)
    {
        aib[pos] += val;
        pos += (pos & -pos);
    }
}
long long fenwick_sum(int pos)
{
    int s = 0;
    while (pos) {
    s += aib[pos];
    pos &= pos - 1;
    }
    return s;
}
long long fenwick_range_sum(int from, int to)
{
    return fenwick_sum(to) - fenwick_sum(from - 1);
}
int fenwick_bin_search(long long sum)
{
    int ans=0,curr_sum=0;
    for(int step=max_p2; step>0; step>>=1)
        if(ans+step<=n && curr_sum+aib[ans+step]<sum)
        {
            ans+=step;
            curr_sum+=aib[ans];
        }
    if(ans+1>n || fenwick_sum(ans+1)!=sum)
        return -1;
    return ans+1;
}

int main()
{
    int op, x, y;
    in>>n>>m;
    for(int i=1; i<=n; i++)
    {
        in>>aib[i];
    }
    fenwick_build();

    for(int i=1; i<=m; i++)
    {
        in>>op;
        if(op==0)
        {
            in>>x>>y;
            fenwick_add(x, y);
        }
        if(op==1)
        {
            in>>x>>y;
            out<<fenwick_range_sum(x, y)<<'\n';
        }
        if(op==2)
        {
            in>>x;
            out<<fenwick_bin_search(x)<<'\n';
        }
    }
}