题目描述
众所周知,模数的hash会产生冲突。例如,如果模的数p=7
,那么4
和11
便冲突了。
B君对hash冲突很感兴趣。他会给出一个正整数序列value[]
。
自然,B君会把这些数据存进hash池。第value[k]
会被存进(k%p)
这个池。这样就能造成很多冲突。
B君会给定许多个p
和x
,询问在模p
时,x
这个池内**数的总和
**。
另外,B君会随时更改value[k]
。每次更改立即生效。
保证1<=p<n1<=p<n.
输入格式
第一行,两个正整数n,m
,其中n
代表序列长度,m
代表B君的操作次数。
第一行,n
个正整数,代表初始序列。
接下来m
行,首先是一个字符cmd
,然后是两个整数x,y
。
-
若
cmd='A'
,则询问在模x
时,y
池内数的总和。 -
若
cmd='C'
,则将value[x]
修改为y
。
输出格式
对于每个询问输出一个正整数,进行回答。
输入输出样例
输入 #1复制
10 5 1 2 3 4 5 6 7 8 9 10 A 2 1 C 1 20 A 3 1 C 5 1 A 5 0
输出 #1复制
25 41 11
说明/提示
样例解释
A 2 1
的答案是1+3+5+7+9=25
.
A 3 1
的答案是20+4+7+10=41
.
A 5 0
的答案是1+10=11
.
数据规模
对于10%
的数据,有n<=1000,m<=1000
.
对于60%
的数据,有n<=100000.m<=100000
.
对于100%
的数据,有n<=150000,m<=150000
.
保证所有数据合法,且1<=value[i]<=1000
.
。。。。。数据有点弱,直接暴力只会T一个点。那么先给个暴力的代码:
#include <bits/stdc++.h>
using namespace std;
const int mac=1e6+10;
int a[mac];
int main()
{
ios::sync_with_stdio(false);
cin.tie(0);cout.tie(0);
int n,m;
cin>>n>>m;
for (int i=1; i<=n; i++)
cin>>a[i];
for (int i=1; i<=m; i++){
char s[2];
int x,y;
cin>>s>>x>>y;
if (s[0]=='A') {
int ans=0;
for (int j=y; j<=n; j+=x) ans+=a[j];
cout<<ans<<endl;
}
else {
a[x]=y;
}
}
return 0;
}
T了一个的点是死活都过不去,除了打表。。。。那么我们就要优化我们的算法了。按照分块一贯的尿性就是将n处理成sqrt(n)。有什么用呢?简单来讲这一题可以预处理,暴力预处理的复杂度是n^2,我们将模降成sqrt(n)那么复杂度就会下降:
for (int i=1; i<=sqrt(n); i++)//枚举模数
for (int j=1; j<=n; j++)
ans[i][j%i]+=a[j];
接下来每次询问就是O(1)复杂度了,修改的时候也是sqrt(n)的复杂度:
for (int j=1; j<=sqrt(n); j++)
ans[j][x%j]-=a[x],ans[j][x%j]+=y;
a[x]=y;
那么还有一个就是当模大于sqrt(n)的时候,这个我们可以直接暴力求,和上面的T了代码是一样的。但由于这个时候模数比较大,大于sqrt(n),所以暴力求的时候的复杂度小于sqrt(n)。那么我们总的复杂度就是O((N+M)*(sqrt(n))。
以下是AC代码:
#include <bits/stdc++.h>
using namespace std;
const int mac=1.5e5+10;
const int N=sqrt(mac);
int a[mac];
int ans[mac][N];
int main()
{
ios::sync_with_stdio(false);
cin.tie(0);cout.tie(0);
int n,m;
cin>>n>>m;
for (int i=1; i<=n; i++)
cin>>a[i];
for (int i=1; i<=sqrt(n); i++)
for (int j=1; j<=n; j++)
ans[i][j%i]+=a[j];
for (int i=1; i<=m; i++){
char s[2];
int x,y;
cin>>s>>x>>y;
if (s[0]=='A') {
if (x<=sqrt(n)) cout<<ans[x][y]<<endl;
else {
int sum=0;
for (int j=y; j<=n; j+=x)
sum+=a[j];
cout<<sum<<endl;
}
}
else {
for (int j=1; j<=sqrt(n); j++)
ans[j][x%j]-=a[x],ans[j][x%j]+=y;
a[x]=y;
}
}
return 0;
}