題意： 給你 N 個數(shù) 和 一個 M；

　　　　對于 每一個 Ni ， 乘以 K 取摸 M 都有一個 集合， 把所有集合合并， 求和

?????????? Σ ai ( ai → K * Ni % M )

思路 ： 　　最開始 直接求一邊gcd ， 然后容斥。。。。 結(jié)果狀態(tài)有 2 ^? （1e4）....

　　　　　　反著求 M 的約數(shù)， 然后記錄要用到的約數(shù)， 對于這些進(jìn)行容斥就好了（不能狀壓）

　　　　　　

#include<bits/stdc++.h> using namespace std; typedef long long LL; const int maxn = 1e5 + 131; LL gcd(LL a, LL b){return b == 0 ? a : gcd(b, a % b); }LL P[maxn], cnt; void GetP(LL m){cnt = 0;for(LL i = 1; i <= sqrt(m); ++i){if(m % i == 0){P[cnt++] = i;if(i * i != m) P[cnt++] = m / i;}}sort(P,P+cnt); }LL Vis[maxn], Num[maxn];int main() {int T;scanf("%d", &T);for(int kase = 1; kase <= T; ++kase){LL n, m;scanf("%lld %lld", &n, &m);GetP(m);LL u;memset(Vis,0, sizeof(Vis));memset(Num,0, sizeof(Num));for(int i = 0; i < n; ++i){scanf("%lld", &u);LL tmp = gcd(u, m);for(LL j = 0; j < cnt; ++j){if(P[j] % tmp == 0)Vis[j] = 1;}}//cout << P[cnt - 1 ] << endl;Vis[cnt -1] = 0;LL Ans = 0;for(LL i = 0; i < cnt; ++i) // 容斥 {if(Vis[i] != Num[i]){LL tmp = m / P[i], D = Vis[i] - Num[i];Ans += (tmp + 1) * tmp / 2 * P[i] * D;for(LL j = i; j < cnt; ++j){if(P[j] % P[i] == 0)Num[j] += D;} }}printf("Case #%d: %lld\n", kase, Ans);} }

?

轉(zhuǎn)載于:https://www.cnblogs.com/aoxuets/p/4935421.html

《新程序員》：云原生和全面數(shù)字化實踐50位技術(shù)專家共同創(chuàng)作，文字、視頻、音頻交互閱讀

總結(jié)

以上是生活随笔為你收集整理的HDU 5514的全部內(nèi)容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網(wǎng)站內(nèi)容還不錯，歡迎將生活随笔推薦給好友。