Tishreen-CPC 2018 G. Colors Overflow(分块)
Input
The first line contains a single integer T, denoting the number of test cases.
Each test case start with a line containing an integer Q (1 ≤ Q ≤ 105), denoting the number of queries.
The next line contains Q lines describe the queries. Each line has one query of the form "t L R"(1 ≤ t ≤ 2)
(1 ≤ L ≤ R ≤ 105).
Output
For each test case and for each query of the second type, print a single line, indicating the number of
unique colors in the range [L, R].
思路:
将每种颜色看成一条线段,则每次1操作相当于在区间内画了若干条线段,2操作相当于询问指定区间内线段数。
则对于每次询问,统计[1,L)的右端点数和[1-R]的左端点数,后者减前者,即得到答案(具体可画图理解),而每次1操作增加的左端点数为递减等差数列,增加的右操作数为与之倒序的递增等差数列。
参考了HH学长的实现思路,利用这个性质,用分块的思想来维护区间内左右端点数。
在实现上,对于全部询问Q,第一次,我们只更新右端点数,并在询问操作时把指定区间的右端点数减入到对应ans里;第二次,我们只更新左端点数,并在询问操作时把指定区间内左端点数加入到对应ans里。其中第一次完成后,将整个区间左右旋转,并初始化,使得不同但对称的操作,代码能够复用,很巧妙~~具体实现见代码
#include<iostream>
#include<algorithm>
#define de(x) cout<< #x <<" = "<<x<<endl
using namespace std;
typedef long long ll;
const int bsz=350,N=1e5+7,n=1e5+2;
int bn,bl[N],br[N],sz[N];
int qop[N],ql[N],qr[N];
ll ans[N],bs[N],bb[N],bd[N],a[N];
void init(int n)
{
bn=(n-1)/bsz+1;
for (int i=0;i<bn;++i)
{
bl[i]=i*bsz;
br[i]=min(n,(i+1)*bsz);
sz[i]=br[i]-bl[i];
bs[i]=bd[i]=bb[i]=0;
}
for (int i=0;i<n;++i)
a[i]=0;
}
void build(int b)
{
for (int i=bl[b];i<br[b];++i)
a[i]+=bb[b]+(i-bl[b]+1)*bd[b];
bb[b]=bd[b]=0;
}
void update(int b,int l,int r)
{
int nl=max(l,bl[b]),nr=min(br[b],r);
if (nl>=nr)
return;
if (nr-nl==sz[b])
{
bb[b]+=nl-l; //维护这一块的基数
bd[b]++; //维护这一块的公差
bs[b]+=(1+sz[b])*sz[b]/2+(nl-l)*sz[b]; //用求和公式维护这一块端点数
return;
}
for (int i=nl;i<nr;++i) //如果该块没有被包含,则直接暴力维护
{
bs[b]+=i-l+1;
a[i]+=i-l+1;
}
return;
}
ll query(int b,int l,int r)
{
int nl=max(l,bl[b]),nr=min(br[b],r);
if (nl>=nr)
return 0;
if (nr-nl==sz[b]) //被包含直接返回这一块的值
return bs[b];
//没有被包含时
ll res=0;
build(b);//通过之前维护的公差和基数暴力统计端点数
for (int i=nl;i<nr;++i)
res+=a[i];
return res;
}
int main()
{
std::ios::sync_with_stdio(0);
cin.tie(0);
int T;
cin>>T;
while (T--)
{
int q;
cin>>q;
fill_n(ans,q,0);
for (int i=0;i<q;++i)
{
cin>>qop[i]>>ql[i]>>qr[i];
--ql[i];
--qr[i];
}
init(n);
for (int i=0;i<q;++i)
{
if (qop[i]==1)
for (int j=0;j<bn;++j)
update(j,ql[i],qr[i]+1);
else
for (int j=0;j<bn;++j)
ans[i]-=query(j,0,ql[i]);
//旋转
ql[i]=n-1-ql[i];
qr[i]=n-1-qr[i];
swap(ql[i],qr[i]);
}
init(n);
for (int i=0;i<q;++i)
{
if (qop[i]==1)
for (int j=0;j<bn;++j)
update(j,ql[i],qr[i]+1);
else
{
for (int j=0;j<bn;++j)
ans[i]+=query(j,ql[i],n);
cout<<ans[i]<<endl;
}
}
}
return 0;
}