首页 > 代码库 > 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
声明:以上内容来自用户投稿及互联网公开渠道收集整理发布,本网站不拥有所有权,未作人工编辑处理,也不承担相关法律责任,若内容有误或涉及侵权可进行投诉: 投诉/举报 工作人员会在5个工作日内联系你,一经查实,本站将立刻删除涉嫌侵权内容。