洛谷P2429 制杖题 [2017年6月计划 数论10]

P2429 制杖题

题目描述

求不大于 m 的、 质因数集与给定质数集有交集的自然数之和。

输入输出格式

输入格式：

第一行二个整数 n，m。

第二行 n 个整数，表示质数集内的元素 p[i]。

输出格式：

一个整数，表示答案，对 376544743 取模。

输入输出样例

输入样例#1：

2 15
3 5

输出样例#1：

60

说明

样例解释：所有符合条件的数为 3,5,6,9,10,12,15 其和为 60。

··· 测试点编号 规模

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
···

前三个点：n * m <= 1e7

不难想到暴力求解

后七个点：n <= 20,m <= 1e9

容斥+等差数列求和

利用二进制枚举各个数的乘积，利用等差数列求和，容斥原理排除多算的即可

#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <algorithm>
#include <vector>
#include <queue>
inline void read(long long &x){x = ;char ch = getchar();char c = ch;while(ch > '' || ch < '')c = ch, ch = getchar();while(ch >= '' && ch <= '')x = x *  + ch - '', 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 =  + ;
const int MOD =  ;

long long n,m;
long long p[MAXN];
long long ans;

int main()
{
    read(n);read(m);
    for (register long long i =  ; i <= n ; ++ i)
        read(p[i]);
    if(n * m <= )
    {
        for(int i = ;i <= m;i ++)
        {
            for(int j = ;j <= n;j ++)
            {
                if(!(i % p[j]))
                {
                    ans += i;
                    ans %= MOD;
                    break;
                }
            }
        }
        printf("%lld", ans % MOD);
        return ;
    }

    int S = ( << n);
    register long long num,cnt,x;
    register long long niyuan = (MOD + )/ ;

    for (register long long i =  ; i < S ; ++ i)
    {
        register long long num = ,cnt = ,x = ;
        for (register long long j = , k = i ; k ; ++j, k >>= )
            if (k & )num *= p[j],x ++;
        cnt = m / num;
        if (x)
        {
            if (x & )
                ans += (((num * ( + cnt))%MOD * cnt)%MOD * niyuan) % MOD;
            else
                ans -= (((num * ( + cnt))%MOD * cnt)%MOD * niyuan) % MOD;
            ans = ans % MOD;
        }
    }
    printf("%lld", ans % MOD);
    return ;
}