首页 > 代码库 > hdu6069[素数筛法] 2017多校3

hdu6069[素数筛法] 2017多校3

/*hdu6069[素数筛法] 2017多校3*/
#include <bits/stdc++.h>
using namespace std;
typedef long long LL;
LL l, r, k;
const LL MOD = 998244353LL;
int T, n, prime[1100000], primesize;
bool isprime[11000000];
void getlist(int listsize)
{
    memset(isprime, 1, sizeof(isprime));
    isprime[1] = false;
    for (int i = 2; i <= listsize; i++)
    {
        if (isprime[i])prime[++primesize] = i;
        for (int j = 1; j <= primesize && i * prime[j] <= listsize; j++)
        {
            isprime[i * prime[j]] = false;
            if (i % prime[j] == 0)break;
        }
    }
}
LL num[1000005], ans[1000005];
void solve() {
    LL n = r - l + 1;
    for (int i = 0; i < n; i++) {
        num[i] = i + l;
        ans[i] = 1;
    }
    for (int i = 1; (LL)prime[i]*prime[i] <= r; i++) {
        for (LL j = prime[i] * (l / prime[i]); j <= r; j += prime[i]) {
            if (j < l) continue;
            LL cnt = 0;
            while (num[j - l] % prime[i] == 0) {
                cnt++;
                num[j - l] /= prime[i];
            }
            ans[j - l] = (ans[j - l] * (1LL + cnt * k)) % MOD;
        }
    }
    LL res = 0;
    for (int i = 0; i < n; i++) {
        if (num[i] > 1) {
            ans[i] = (ans[i] * (1LL + k)) % MOD;
        }
        res = (res + ans[i]) % MOD;
    }
    printf("%lld\n", res);
}
int main() {
    getlist(1000005);
    scanf("%d", &T);
    while (T--) {
        scanf("%lld%lld%lld", &l, &r, &k);
        solve();
    }
    return 0;
}

 

hdu6069[素数筛法] 2017多校3