最近又去研究暴力美学。。
深深的感到暴力之美。
题目1:给出一个长为 n 的数列,以及 n个操作,操作涉及区间加法,单点查值。
思路:把数列分成根号n块 对于一个区间l r 先暴力处理掉连边不包含整个块的,然后对包含整个块的直接打上标记
#include<iostream>
#include<cstdio>
#include<vector>
#include<queue>
#include<set>
#include<cmath>
#include<algorithm>
#include<cstring>
using namespace std;
#define rep(i,j,k) for(int i=j;i<=k;i++)
#define sca(x) scanf("%d",&x)
#define per(i,j,k) for(int i=j;i>=k;i--)
#define inf 0x3f3f3f3f
#define LL long long
#define N 550005
#define inf 0x3f3f3f3f
int len;
int a[N],b[N],tag[N];
int read()
{
char c;
int x=0,f=1;
c=getchar();
while(c<'0'||c>'9')
{
if(c=='-')f=-1;
c=getchar();
}
while(c>='0'&&c<='9')
{
x=x*10+c-'0';
c=getchar();
}
return x*f;
}
void add(int l,int r,int k)
{
for(int i=l;i<=min(b[l]*len,r);i++)a[i]+=k;
if(b[l]!=b[r])
{
for(int i=(b[r]-1)*len+1;i<=r;i++)a[i]+=k;
}
for(int i=b[l]+1;i<=b[r]-1;i++)tag[i]+=k;
}
int main()
{
int n;
int op,l,r,k;
n=read();
len=sqrt(n);
rep(i,1,n)a[i]=read();
rep(i,1,n)b[i]=(i-1)/len+1;
rep(i,1,n)
{
op=read(),l=read(),r=read(),k=read();
if(op==0)
{
add(l,r,k);
}
else
{
printf("%d\n",a[r]+tag[b[r]]);
}
}
return 0;
}
题目2:给出一个长为 n的数列,以及 n 个操作,操作涉及区间加法,询问区间内小于某个值 x 的元素个数
思路:依旧是分块对于每个块内的数先排好序,不在整个块内的依然是暴力,对于在块内的就二分查询
注意不完整的块更新完后要重新排序。
#include<iostream>
#include<cstdio>
#include<vector>
#include<queue>
#include<set>
#include<cmath>
#include<algorithm>
#include<cstring>
using namespace std;
#define rep(i,j,k) for(int i=j;i<=k;i++)
#define sca(x) scanf("%d",&x)
#define per(i,j,k) for(int i=j;i>=k;i--)
#define inf 0x3f3f3f3f
#define LL long long
#define N 550005
#define pb(x) push_back(x)
#define inf 0x3f3f3f3f
int len;
int a[N],b[N],tag[N];
vector<int>v[25000];
int read()
{
char c;int x=0,f=1;
c=getchar();
while(c<'0'||c>'9'){if(c=='-')f=-1;c=getchar();}
while(c>='0'&&c<='9'){x=x*10+c-'0';c=getchar();}
return x*f;
}
void modify(int x)
{
v[x].clear();
for(int i=(x-1)*len+1;i<=x*len;i++)v[x].pb(a[i]);
sort(v[x].begin(),v[x].end());
}
void add(int l,int r,int k)
{
for(int i=l;i<=min(b[l]*len,r);i++)a[i]+=k;
modify(b[l]);//l所在的整个块重新排序
if(b[l]!=b[r])
{
for(int i=(b[r]-1)*len+1;i<=r;i++)a[i]+=k;
modify(b[r]);//r所在的块重新排序
}
for(int i=b[l]+1;i<=b[r]-1;i++)tag[i]+=k;
}
int query(int l,int r,int k)
{
int ans=0;
for(int i=l;i<=min(b[l]*len,r);i++)
{
if(a[i]+tag[b[l]]<k)ans++;
}
if(b[l]!=b[r])
{
for(int i=(b[r]-1)*len+1;i<=r;i++)
{
if(a[i]+tag[b[r]]<k)ans++;
}
}
for(int i=b[l]+1;i<=b[r]-1;i++)
{
int tmp=k-tag[i];
ans+=lower_bound(v[i].begin(),v[i].end(),tmp)-v[i].begin();
}
return ans;
}
int main()
{
int n;
n=read();
int op,l,r,c;
len=sqrt(n);
rep(i,1,n)a[i]=read();
rep(i,1,n)b[i]=(i-1)/len+1;
rep(i,1,n)v[b[i]].pb(a[i]);
rep(i,1,b[n])sort(v[i].begin(),v[i].end());//先对每个块排序
rep(i,1,n)
{
op=read(),l=read(),r=read(),c=read();
if(op==0)add(l,r,c);
else printf("%d\n",query(l,r,c*c));
}
}
题目3:给出一个长为 n 的数列,以及 n 个操作,操作涉及区间加法,询问区间内小于某个值 x 的前驱(比其小的最大元素)。
思路:与上一题基本上相似,这个要熟悉lower bound 的用法,测试数据里边应该是没有负数的 所以ans的值直接设位-1.
#include<iostream>
#include<cstdio>
#include<vector>
#include<queue>
#include<set>
#include<cmath>
#include<algorithm>
#include<cstring>
using namespace std;
#define rep(i,j,k) for(int i=j;i<=k;i++)
#define sca(x) scanf("%d",&x)
#define per(i,j,k) for(int i=j;i>=k;i--)
#define inf 0x3f3f3f3f
#define LL long long
#define N 550005
#define pb(x) push_back(x)
#define inf 0x3f3f3f3f
int len;
int a[N],b[N],tag[N];
vector<int>v[25000];
int read()
{
char c;int x=0,f=1;
c=getchar();
while(c<'0'||c>'9'){if(c=='-')f=-1;c=getchar();}
while(c>='0'&&c<='9'){x=x*10+c-'0';c=getchar();}
return x*f;
}
void modify(int x)
{
v[x].clear();
for(int i=(x-1)*len+1;i<=x*len;i++)v[x].pb(a[i]);
sort(v[x].begin(),v[x].end());
}
void add(int l,int r,int k)
{
for(int i=l;i<=min(b[l]*len,r);i++)a[i]+=k;
modify(b[l]);
if(b[l]!=b[r])
{
for(int i=(b[r]-1)*len+1;i<=r;i++)a[i]+=k;
modify(b[r]);
}
for(int i=b[l]+1;i<=b[r]-1;i++)tag[i]+=k;
}
int query(int l,int r,int k)
{
int ans=-1;
for(int i=l;i<=min(b[l]*len,r);i++)
{
if(a[i]+tag[b[l]]<k)ans=max(ans,a[i]+tag[b[l]]);
}
if(b[l]!=b[r])
{
for(int i=(b[r]-1)*len+1;i<=r;i++)
{
if(a[i]+tag[b[r]]<k)ans=max(ans,a[i]+tag[b[r]]);
}
}
for(int i=b[l]+1;i<=b[r]-1;i++)
{
int tmp=k-tag[i];
int pos=lower_bound(v[i].begin(),v[i].end(),tmp)-v[i].begin();
if(pos>=1)ans=max(ans,v[i][pos-1]+tag[i]);
}
return ans;
}
int main()
{
int n;
n=read();
int op,l,r,c;
len=sqrt(n);
rep(i,1,n)a[i]=read();
rep(i,1,n)b[i]=(i-1)/len+1;
rep(i,1,n)v[b[i]].pb(a[i]);
rep(i,1,b[n])sort(v[i].begin(),v[i].end());
rep(i,1,n)
{
op=read(),l=read(),r=read(),c=read();
if(op==0)add(l,r,c);
else printf("%d\n",query(l,r,c));
}
}