题目来源: 鲸歌的博客
题解
这是一道简单的题目,但是朴素的算法是$ ext{O}left( n^2 ight)$的,只能过$70\%$的点.如何优化呢?
考虑在什么情况下会使一段区间内出现过的A,B,C数量相等.显然的,设$totAleft[ i ight]$为前i个字符中A出现的次数(以此类推,$totB$和$totC$)
设一段相等的区间起始于$i$,终止于$j$,显然,当且仅当$totA[j]-totA[i]=totB[j]-totB[i]=totC[j]-totC[i]$时,i与j间A,B,C数相等
移项易得$totA[j]-totB[j]=totA[i]-totB[i],totB[j]-totC[j]=totB[i]-totC[i]$,易证该式成立为上式成立充分必要条件.
这两个值拥有相当有趣的推论,即记$ab[i]=totA[i]-totB[i],bc[i]=totB[i]-totC[i]$,
只要任何$i<j$且$ab[i]=ab[j], bc[i]=bc[j]$成立,i~j 间就有相同的A,B,C数目,反之,是不能有相同的A,B,C数目的.
这时算法仍然是$ ext{O}left( n^2
ight)$,但是通过组合学的方法,我们可以将所有$ab[i],bc[i]$相同的点分为一类,每类中可行段的总和即是答案.而每类中只要是两两元素间的段
都是可行段,即段数为$frac{nleft( n-1
ight)}{2}$.n只需要扫描到此类节点时自增一即可.
最后问题,如何储存一类的和?使用hash表,平均复杂度接近常数,不过hash函数要写的好些.(有冲突怎么办?记得前向星么?就那么处理)
最后附上渣代码,不过效率比std高多了,std评测415ms,我的代码180ms.(有个神犇同学写的AC代码总共1500 ms...)
#include <cstdio>
#define nxp 1110007
#define par 111029
int tot[3],hash[nxp],t,t2,t3,t4,ln;
char c;
struct nnode{
int ab,bc,cn,sz;
} nn[2000000];
inline int abs(int n){
return (n<0?-n:n);
}
int hashF(int ab,int bc){
return abs((ab*par+bc)%nxp);
}
int find(int ab,int bc){
t4=hashF(ab,bc);
t=hash[t4];
while(nn[t].ab!=ab||nn[t].bc!=bc){
if(t==0) return ~t4;//-t4-1
t=nn[t].cn;
}
return t;
}
long long sum;
int main(){
ln=1;
nn[ln].ab=0;
nn[ln].bc=0;
t3=~t3;
nn[ln].cn=hash[0];
nn[ln].sz=1;
hash[0]=ln;
++ln;
while(c=getchar(),c!=EOF&&c!='
'&&c!=' '){
t2=c-'A';
if(t2<3){
++tot[t2];
}
t3=find(tot[1]-tot[0],tot[2]-tot[1]);
if(t3<0){//if has significate bit true
nn[ln].ab=tot[1]-tot[0];
nn[ln].bc=tot[2]-tot[1];
t3=~t3;
nn[ln].cn=hash[t3];
nn[ln].sz=1;
hash[t3]=ln;
++ln;
}else{
sum+=nn[t3].sz;
++nn[t3].sz;
}
}
printf("%d
",sum);
return 0;
}
最后说一句,千万注意$ab[0]=bc[0]=0$,不然很可能错.(即是我代码main最前面如此伪和的来历).唉,解释了这么多不如代码实在啊...我这个代码用多次加法(每次sum加上此点前符合要求的起点总数)代替乘法,附上乘法版(不过这个优化根本不明显啊...)
#include <cstdio>
#define nxp 1110007
#define par 111029
int tot[3],hash[nxp],t,t2,t3,t4,ln,i;
char c;
struct nnode{
int ab,bc,cn,sz;
} nn[2000000];
inline int abs(int n){
return (n<0?-n:n);
}
int hashF(int ab,int bc){
return abs((ab*par+bc)%nxp);
}
int find(int ab,int bc){
t4=hashF(ab,bc);
t=hash[t4];
while(nn[t].ab!=ab||nn[t].bc!=bc){
if(t==0) return ~t4;//-t4-1
t=nn[t].cn;
}
return t;
}
long long sum;
int main(){
ln=1;
nn[ln].ab=0;
nn[ln].bc=0;
t3=~t3;
nn[ln].cn=hash[0];
nn[ln].sz=1;
hash[0]=ln;
++ln;
while(c=getchar(),c!=EOF&&c!='
'&&c!=' '){
t2=c-'A';
if(t2<3){
++tot[t2];
}
t3=find(tot[1]-tot[0],tot[2]-tot[1]);
if(t3<0){//if has significate bit true
nn[ln].ab=tot[1]-tot[0];
nn[ln].bc=tot[2]-tot[1];
t3=~t3;
nn[ln].cn=hash[t3];
nn[ln].sz=1;
hash[t3]=ln;
++ln;
}else{
++nn[t3].sz;
}
}
for(i=1;i<ln;++i) sum+=(nn[i].sz)*(nn[i].sz-1)>>1;
printf("%lld
",sum);
return 0;
}
代码中nxp和par即为hash函数常数了...这里选用了两个大小刚刚好的素数.~x是-x-1的一个很方便的代替方法(而且非常快).为什么-x-1主要是考虑到-0=0的情况不好分辨.(当然~x比-x快是显然的,无需取反后加一)