1. 程式人生 > >【CF472G】Design Tutorial 壓位

【CF472G】Design Tutorial 壓位

一個 for code signed 位置 blog 子串 系列 inline

題目大意

  給出兩個\(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}\)

壓成一個\(64\)位整數,每次暴力統計。

  時間復雜度:\(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 壓位