【题解】CF1320D Reachable Strings

思路

对于 1 并不好考虑，可以针对 0 考虑。

发现每次移动会在一堆 1 之间移动把一个 0 移动 \(2\) 位，因此只需要记录所有 0 位置的奇偶性即可，具体实现可以用 Hash，对开头奇数/偶数分别开 Hash 表记录。

这里 Hash 实现有一个方法：(BKDRHash)

\[h_i = ph_{i-1} + s_i\]

这样可以非常容易地求出区间 Hash 值。（当然这里没写）

不要 unsigned long long 自然溢出啊，会被卡爆的 ......

Code

#include <cstdio>
#include <cctype>
#include <cstring>
#include <algorithm>
using namespace std;
#define File(s) freopen(s".in", "r", stdin), freopen(s".out", "w", stdout)
typedef long long ll;
typedef unsigned long long ull;
template<typename T> inline void gi(T &x){
    char ch; x = 0;
    int f = 1;
    while(isspace(ch = getchar()));
    if(ch == '-') ch = getchar(), f = -1;
    do x = x * 10 + ch - '0'; while(isdigit(ch = getchar()));
    x *= f;
}
template<class T> void upmax(T &x, T y){x = x>y ? x : y;}
template<class T> void upmin(T &x, T y){x = x<y ? x : y;}

const int N = 200005;
const ull p1 = 999983, p2 = 10007, mod1 = 104857601, mod2 = 167772161;
char s[N];
int cnt[N];

template<ull p, ull mod>
struct Hash{
    ull h[N], pw[N];
    int len;
    Hash(){pw[0] = 1; len = 0; h[0] = 0;}
    void hash(ull x){
        ++len;
        pw[len] = pw[len - 1] * p % mod;
        h[len] = (h[len - 1] + x * pw[len]) % mod;
    }
    void space(){
        ++len;
        h[len] = h[len - 1];
        pw[len] = pw[len - 1];
    }
    ull ask(int l, int r){
        return (h[r] - h[l - 1] + mod) % mod;
    }
};
Hash<p1, mod1> h01, h11;
Hash<p2, mod2> h02, h12;

int main(){
    // File("cf1320d");
    int n;
    gi(n);
    scanf("%s", s + 1);
    for(int i=1; i<=n; i++){
        if(s[i] == '1'){
            h01.space(); h11.space(); h02.space(); h12.space();
            cnt[i] = cnt[i-1];
        }
        else{
            h01.hash((i & 1) + 1); h02.hash((i & 1) + 1);
            h11.hash(((i & 1) ^ 1) + 1); h12.hash(((i & 1) ^ 1) + 1);
            cnt[i] = cnt[i-1] + 1;
        }
    }
    int q; gi(q);
    while(q--){
        int _l, _r, _len;
        gi(_l); gi(_r); gi(_len);
        int l1, r1, l2, r2;
        l1 = _l; r1 = _l  + _len - 1; l2 = _r; r2 = _r + _len - 1;
        if((cnt[r1] - cnt[l1 - 1]) != (cnt[r2] - cnt[l2 - 1])){
            puts("No");
            continue;
        }
        if((l1 ^ l2) & 1){
            ull h1 = h01.ask(l1, r1) * h11.pw[l2 - 1] % mod1;
            ull h2 = h11.ask(l2, r2) * h01.pw[l1 - 1] % mod1;
            if(h1 != h2){
                puts("No");
                continue;
            }
            h1 = h02.ask(l1, r1) * h12.pw[l2 - 1] % mod2;
            h2 = h12.ask(l2, r2) * h02.pw[l1 - 1] % mod2;
            if(h1 != h2){
                puts("No");
                continue;
            }
        }
        else{
            ull h1 = h01.ask(l1, r1) * h01.pw[l2 - 1] % mod1;
            ull h2 = h01.ask(l2, r2) * h01.pw[l1 - 1] % mod1;
            if(h1 != h2){
                puts("No");
                continue;
            }
            h1 = h02.ask(l1, r1) * h02.pw[l2 - 1] % mod2;
            h2 = h02.ask(l2, r2) * h02.pw[l1 - 1] % mod2;
            if(h1 != h2){
                puts("No");
                continue;
            }
        }
        puts("Yes");
    }
    return 0;
}