P2429 制杖題

題目描述

求不大于 m 的、 質(zhì)因數(shù)集與給定質(zhì)數(shù)集有交集的自然數(shù)之和。

輸入輸出格式

輸入格式：

第一行二個(gè)整數(shù) n，m。

第二行 n 個(gè)整數(shù)，表示質(zhì)數(shù)集內(nèi)的元素 p[i]。

輸出格式：

一個(gè)整數(shù)，表示答案，對(duì) 376544743 取模。

輸入輸出樣例

輸入樣例#1：

2 15
3 5

輸出樣例#1：

60

說(shuō)明

樣例解釋：所有符合條件的數(shù)為 3,5,6,9,10,12,15 其和為 60。

··· 測(cè)試點(diǎn)編號(hào) 規(guī)模

1 2 3 n*m<=10^7
4 5 n<=2，m<=10^9
6 7 n<=20，m<=10^8
8 9 10 n<=20，m<=10^9
···

?

前三個(gè)點(diǎn)：n * m <= 1e7

不難想到暴力求解

?

后七個(gè)點(diǎn)：n <= 20,m <= 1e9

容斥+等差數(shù)列求和

利用二進(jìn)制枚舉各個(gè)數(shù)的乘積，利用等差數(shù)列求和，容斥原理排除多算的即可

?

#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <cstring> 
#include <algorithm>
#include <vector>
#include <queue>
inline void read(long long &x){x = 0;char ch = getchar();char c = ch;while(ch > '9' || ch < '0')c = ch, ch = getchar();while(ch >= '0' && ch <= '9')x = x * 10 + ch - '0', ch = getchar();}
inline int max(int a, int b){return a > b ? a : b;}
inline int min(int a, int b){return a < b ? a : b;}
inline void swap(int &a, int &b){int tmp = a;a = b;b = tmp;}const int INF = 0x3f3f3f3f;
const int MAXN = 2000 + 5;
const int MOD =  376544743;long long n,m;
long long p[MAXN];
long long ans;int main()
{read(n);read(m);for (register long long i = 1 ; i <= n ; ++ i)read(p[i]);if(n * m <= 100000000){for(int i = 1;i <= m;i ++){for(int j = 1;j <= n;j ++){if(!(i % p[j])){ans += i;ans %= MOD;break;}}}printf("%lld", ans % MOD);return 0;}int S = (1 << n);register long long num,cnt,x;register long long niyuan = (MOD + 1)/ 2;for (register long long i = 1 ; i < S ; ++ i){register long long num = 1,cnt = 0,x = 0;for (register long long j = 1, k = i ; k ; ++j, k >>= 1)if (k & 1)num *= p[j],x ++;cnt = m / num;if (x){if (x & 1)ans += (((num * (1 + cnt))%MOD * cnt)%MOD * niyuan) % MOD;elseans -= (((num * (1 + cnt))%MOD * cnt)%MOD * niyuan) % MOD;ans = ans % MOD;}}printf("%lld", ans % MOD);return 0;
}

?

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

總結(jié)

以上是生活随笔為你收集整理的洛谷P2429 制杖题 [2017年6月计划 数论10]的全部?jī)?nèi)容，希望文章能夠幫你解決所遇到的問(wèn)題。

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