【CF472G】Design Tutorial 壓位
阿新 • • 發佈:2018-03-05
一個 for code signed 位置 blog 子串 系列 inline 壓成一個\(64\)位整數,每次暴力統計。
題目大意
給出兩個\(01\)序列\(A\)和\(B\)
漢明距離定義為兩個長度相同的序列中,有多少個對應位置上的數字不一樣
\(00111\) 和 \(10101\)的距離為\(2\)
\(Q\)次詢問,每次詢問給出\(p_1,p_2,len\)
求\(a_{p_1},a_{p_1+1}\ldots a_{p_1+len?1}\)和\(b_{p_1},b_{p_1+1}\ldots b_{p_1+len?1}\)兩個子串的漢明距離
\(n\leq 2\times{10}^5,q\leq 4\times {10}^5\)
題解
wys【挑戰】弱化版
暴力碾分塊系列
把\(a_x\ldots a_{x+63}\)
時間復雜度:\(O(\frac{nq}{64})\)
代碼
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cstdlib>
#include<ctime>
#include<utility>
using namespace std;
typedef long long ll;
typedef unsigned long long ull;
typedef pair<int,int> pii;
char s1[200010];
char s2[200010];
ull a[200010];
ull b[200010];
int cnt[100010];
int count(ull x)
{
return cnt[x&0xffff]+cnt[(x>>16)&0xffff]+cnt[(x>>32)&0xffff]+cnt[(x>>48)&0xffff];
}
int main()
{
int n,m,q;
scanf("%s%s%d",s1+1,s2+1,&q);
n=strlen(s1+1);
m=strlen(s2+1);
int i;
for(i=n;i>=1;i--)
a[i]=(a[i+1]>>1)|(s1[i]==‘1‘?(1ll<<63):0);
for(i=m;i>=1;i--)
b[i]=(b[i+1]>>1)|(s2[i]==‘1‘?(1ll<<63):0);
cnt[0]=0;
for(i=1;i<(1<<16);i++)
cnt[i]=cnt[i>>1]+(i&1);
int x,y,l;
for(i=1;i<=q;i++)
{
scanf("%d%d%d",&x,&y,&l);
x++;
y++;
int ans=0;
while(l>=64)
{
ans+=count(a[x]^b[y]);
x+=64;
y+=64;
l-=64;
}
if(l)
ans+=count((a[x]^b[y])>>(64-l));
printf("%d\n",ans);
}
return 0;
}
【CF472G】Design Tutorial 壓位