标题:Log大侠

atm参加了速算训练班,经过刻苦修炼,对以2为底的对数算得飞快,人称Log大侠。

一天,Log大侠的好友 drd 有一些整数序列需要变换,Log大侠正好施展法力...

变换的规则是: 对其某个子序列的每个整数变为: [log_2 (x) + 1] 其中 [] 表示向下取整,就是对每个数字求以2为底的对数,然后取下整。
例如对序列 3 4 2 操作一次后,这个序列会变成 2 3 2。

drd需要知道,每次这样操作后,序列的和是多少。

【输入格式】
第一行两个正整数 n m 。
第二行 n 个数,表示整数序列,都是正数。
接下来 m 行,每行两个数 L R 表示 atm 这次操作的是区间 [L, R],数列序号从1开始。

【输出格式】
输出 m 行,依次表示 atm 每做完一个操作后,整个序列的和。

  这题暴力肯定可以得一部分的,区间大小最差的情况就是每次都给L=1,R=N,这时暴力肯定不行的,会写线段树的话,可以看出这题也就是单点更新以及总区间查询。而它没说数据范围,但就算是1e18的话,取log对于每个位置的数来说,它最多也就更新个60次左右,所以我们可以加个标记,代表这个区间内还有没有位置需要更新,然后再线段树维护,时间复杂度就是60*nlog(n)这样。

 #include<cstdio>
#include<cmath>
#define L(x) (x<<1)
#define R(x) (x<<1|1)
#define M(x) ((T[x].l+T[x].r)>>1)
typedef long long ll;
const int N=;
struct Tree{
bool flag;
int l,r;
ll sum;
}T[N<<];
ll a[N];
void built(int id,int l,int r)
{
T[id].l=l;
T[id].r=r;
T[id].sum=;
T[id].flag=false;
if(l==r)
{
T[id].sum=a[l];
T[id].flag=(a[l]>2ll);
return ;
}
built(L(id),l,M(id));
built(R(id),M(id)+,r);
T[id].flag=T[L(id)].flag|T[R(id)].flag;
T[id].sum=T[L(id)].sum+T[R(id)].sum;
}
void modify(int id,int l,int r)
{
if(!T[id].flag)
return ;
if(T[id].l==T[id].r)
{
T[id].sum=(ll)floor(log2(1.0*T[id].sum)+1.0);
T[id].flag=(T[id].sum>2ll);
return ;
}
if(l<=M(id))
modify(L(id),l,r);
if(r>M(id))
modify(R(id),l,r);
T[id].flag=T[L(id)].flag|T[R(id)].flag;
T[id].sum=T[L(id)].sum+T[R(id)].sum;
}
int main()
{
int n,m,l,r;
scanf("%d%d",&n,&m);
for(int i=;i<=n;i++)
scanf("%lld",&a[i]);
built(,,n);
while(m--)
{
scanf("%d%d",&l,&r);
modify(,l,r);
printf("%lld\n",T[].sum);
}
return ;
}

线段树下线段果

  不懂线段树的话,知道stl的map的话,还有种map的写法,思路一样,当某个位置的值<=2时就把它从map删去。

 #include<cstdio>
#include<cmath>
#include<map>
using namespace std;
typedef long long ll;
map<int,ll> mmp;
map<int,ll>::iterator b,e,temp;
int main()
{
int n,m,l,r;
ll sum=;
scanf("%d%d",&n,&m);
for(int i=;i<=n;i++)
{
scanf("%lld",&mmp[i]);
sum+=mmp[i];
}
while(m--)
{
scanf("%d%d",&l,&r);
b=mmp.lower_bound(l);
e=mmp.upper_bound(r);
while(b!=e)
{
temp=b;
b++;
sum-=temp->second;
temp->second=(ll)floor(1.0*log2(temp->second)+1.0);
sum+=temp->second;
if(temp->second<=)
mmp.erase(temp);
}
printf("%lld\n",sum);
}
return ;
}

STLwd

05-20 09:27