Gym - 102056C(2018EC final) -Heretical … Möbius ——CRT

题意

给出一个长为200的01序列，判断是否在前1e9个莫比乌斯*值中。（这里的莫比乌斯值加了绝对值）

分析

意到因为4的倍数一定是0，9的倍数一定是0……169的倍数一定是0。那么我们可以对4,9,25,49,121,169这6个200以内这质数平方进行考虑。

我们可以枚举起点位置 $x$ 对这6个数的模数，然后用CRT求出 $x$。对每个起点位置，暴力对比即可。不可能存在所有的mu值，只能单个求。

由于0的个数介于65~95，所以符合条件的起点位置并不多，因此不会超时。

#include <bits/stdc++.h>
using namespace std;

typedef long long ll;
const int inf = 0x3f3f3f3f;
const ll up = 1e9;
const int N = 5e4 + ;

int prime[N];
bool notprime[N];
void getprime()
{
    for (int i = ; i < N; i++) {
        if(!notprime[i]) prime[++prime[]] = i;
        for (int j = ; j <= prime[] && i * prime[j] < N; j++)
        {
            notprime[i * prime[j]] = ;
            if(i%prime[j]==) break;
        }
    }
}

int mu(int x) {    //求单个mu值
    for (int i = ; i <= prime[] && prime[i] * prime[i] <= x; i++) {
        if(x%(prime[i]*prime[i]) == ) return ;
        if(x%prime[i] == ) x/=prime[i];
    }
    return ;
}

ll exgcd(ll a,ll b,ll &x, ll &y) {
    ll g = a;
    if(b==) x=,y=; else g=exgcd(b,a%b,y,x),y-=x*(a/b);
    return g;
}

ll inv(ll a,ll m) {
    ll x, y;
    ll d = exgcd(a, m, x, y);
    return (d == ) ? (x % m + m) % m : -;
}

int ti[], m[] = {, , , , , }, a[];
int M;

void CRT() {
    M = ;
    for(int i = ; i < ; i ++) M = M*m[i];
    for(int i = ; i < ; i++) {
        int Mi = M/m[i];
        ti[i] = 1LL*inv(Mi, m[i])*Mi%M;
    }
}

string s,t;

int check(int v, int r) {
    while(v < ) {
        if(s[v] == '') return ;
        v += r;
    }
    return ;
}

int ok(int x) {
    for (int i = ; i < ; i++)
       if(mu(i+x) != s[i]-'') return ;
   return ;
}

int ans = inf;

void dfs(int x) {
    if(x == ) {
        int v = ;
        for(int i = ; i < ; i++) v = (v + 1LL*ti[i]*a[i]%M)%M; //v为CRT的值
        if(v == ) v = M;
        while(v+ <= up && v < ans) {
            if(ok(v)) ans = v;
            v = v + M;
        }
        return;
    }
    for(int i = ; i < m[x]; i++) {  //枚举余数a[i]
        if(check(i, m[x])) {
            a[x] = ( m[x] - i );
            dfs(x+);
        }
    }
}

int main(){
    getprime();
    CRT();
    for(int i = ; i < ; i++) {cin >> t, s = s + t;}
    int num = ;
    for(int i = ; i < s.size(); i++)
        if(s[i] == '') num++;
    if(num <  || num > )        //0的个数在一个范围内
        cout << - << endl;
    else {
        dfs();
        if(ans == inf)
            cout << - << endl;
        else
            cout << ans << endl;
    }
    return ;
}

参考链接：

1. https://blog.csdn.net/BUAA_Alchemist/article/details/86652706

2. https://blog.csdn.net/a1214034447/article/details/86373308

3. https://ac.nowcoder.com/acm/contest/view-submission?submissionId=40622831