Sum - 计蒜客

Posted 2021-01-03 zgglj-com

ps：真的好菜，想到莫比乌斯函数，约数又处理不了。模仿 欧拉筛，莫比乌斯筛的递推，也就是当① i % p[ j ] != 0 ② i % p[ j ] == 0，dp[ tp ] 与 dp[ i ]的关系。

const int N = 20000005;

int n, tot;
int dp[N], p[N];

LL sum[N];
bool use[N];

void Inite() {
    dp[1] = sum[1] = 1;
    rep(i, 2, N) {
       if (!use[i]) {
            p[tot++] = i;
            dp[i] = 2;
       }
       for (int j = 0; j < tot && i * p[j] < N; ++j) {
            int tp = i * p[j];
            use[tp] = 1;
            if (i % p[j] == 0) {
                if ((i / p[j]) % p[j] == 0) dp[tp] = 0;
                else dp[tp] = dp[i] / 2;
                break;
            }
            dp[tp] = dp[i] * 2;
       }
    }
    rep(i, 2, N) sum[i] = sum[i - 1] + dp[i];
}

int main()
{
    Inite();
    BEGIN() {
        sc(n);
        pr(sum[n]);
    }
    return 0;
}