LOJ 6278

#include<cstdio>
#include<cmath>
#include<vector>
#include<algorithm>
#define ll long long
using namespace std;
int n,blo,bl[50005];
ll a[50005],atag[505];
vector<int>ve[505];
void reset(int x)
{
ve[x].clear();
for(int i=(x-1)*blo+1;i<=min(x*blo,n);i++)
{
ve[x].push_back(a[i]);
}
sort(ve[x].begin(),ve[x].end());
}
void add(int l,int r,int c)
{
for(int i=l;i<=min(bl[l]*blo,r);i++)
{
a[i]+=c;
}
reset(bl[l]);
if(bl[l]!=bl[r])
{
for(int i=(bl[r]-1)*blo+1;i<=r;i++)
{
a[i]+=c;
}
reset(bl[r]);
}
for(int i=bl[l]+1;i<=bl[r]-1;i++)
{
atag[i]+=c;
}
}
int qu(int left,int right,int val)
{
int ans=0;
for(int i=left;i<=min(bl[left]*blo,right);i++)
{
if(a[i]+atag[bl[i]]<val)
{
ans++;
}
}
if(bl[left]!=bl[right])
{
for(int i=(bl[right]-1)*blo+1;i<=right;i++)
{
if(a[i]+atag[bl[i]]<val)
{
ans++;
}
}
}
for(int i=bl[left]+1;i<=bl[right]-1;i++)
{
ans+=lower_bound(ve[i].begin(),ve[i].end(),val-atag[i])-ve[i].begin();
}
return ans;
}
int main()
{
scanf("%d",&n);
blo=sqrt(n);
for(int i=1;i<=n;i++)
{
scanf("%lld",&a[i]);
bl[i]=(i-1)/blo+1;
ve[bl[i]].push_back(a[i]);
}
for(int i=1;i<=bl[n];i++)
{
sort(ve[i].begin(),ve[i].end());
}
for(int e=1;e<=n;e++)
{
int mod,l,r,c;
scanf("%d%d%d%d",&mod,&l,&r,&c);
if(mod==0)
{
add(l,r,c);
}
else
{
if(mod==1)
{
printf("%d\n",qu(l,r,c*c));
}
}
}
return 0;
}