栏目分类:
子分类:
返回
名师互学网用户登录
快速导航关闭
当前搜索
当前分类
子分类
实用工具
热门搜索
名师互学网 > IT > 面试经验 > 面试问答

zoj 3868 GCD Expectation

面试问答 更新时间: 发布时间: IT归档 最新发布 模块sitemap 名妆网 法律咨询 聚返吧 英语巴士网 伯小乐 网商动力

zoj 3868 GCD Expectation

#include <iostream>#include <cstdio>#include <cstring>#include <cmath>using namespace std;typedef long long ll;const int maxn = 1e6+5;const ll mod = 998244353;int sum[maxn];ll dp[maxn];ll cal(ll x, int n){    ll s = 1;    while(n > 0) {        if(n & 1) s = (s*x) % mod;        x = (x*x)%mod;        n >>= 1;    }    return s;}int main(){    int T;    scanf("%d", &T);    while(T--) {        int n, k;        scanf("%d %d", &n, &k);        memset(sum, 0, sizeof(sum));        memset(dp, 0, sizeof(dp));        int MAX = 0;        for(int i = 0; i < n; ++i) { int a; scanf("%d", &a); MAX = max(MAX, a); sum[a]++;        }        ll ans = 0;        for(int i = MAX; i >= 1; --i) { int cnt = 0; for(int j = i; j <= MAX; j+=i) {     cnt += sum[j];     dp[i] = (dp[i] - dp[j] + mod) % mod; } dp[i] = ((dp[i]+cal(2, cnt)-1)%mod + mod) % mod; ans = (ans + dp[i]*cal(i, k)) % mod;        }        printf("%lldn", ans);    }    return 0;}
转载请注明:文章转载自 www.mshxw.com
本文地址:https://www.mshxw.com/it/379217.html
我们一直用心在做
关于我们 文章归档 网站地图 联系我们

版权所有 (c)2021-2022 MSHXW.COM

ICP备案号:晋ICP备2021003244-6号