#数论,哈希#洛谷 2429

题目

问在 1 … m 1\dots m 1m内有多少个数是 p 1 , p 2 , … , p n ( 1 ≤ n ≤ 20 ) p_1,p_2,\dots,p_n(1\leq n\leq 20 ) p1,p2,,pn(1n20)的倍数,求这些数的和


分析

首先暴力(20!)应该是过不了的,虽然 O ( 1 0 10 ) O(10^{10}) O(1010)貌似能过,但是还得讲一讲神奇的方法,线性筛,当然同时用哈希维护,然后就没有什么了


代码

#include <cstdio>
#include <algorithm>
#define mod 9973//没错,眼见为实
int n,m,prime[21],h[mod],ans;
void in(int &ans){
    ans=0; char c=getchar();
    while (c<48||c>57) c=getchar();
    while (c>47&&c<58) ans=ans*10+c-48,c=getchar();
}
void print(int ans){
    if (ans>9) print(ans/10);
    putchar(ans%10+48);
}
int locate(int x){//哈希函数
    int pos=x%mod,i=0;
    while (i<mod&&h[(pos+i)%mod]!=x&&h[(pos+i)%mod]) i++;
    return (pos+i)%mod;
}
int main(){
    in(n); in(m);
    for (register int i=1;i<=n;i++) in(prime[i]);
    std::sort(prime+1,prime+1+n);
    for (register int i=1;i<=m/prime[1];i++)
    for (register int j=1;j<=n&&(long long)prime[j]*i<=m;j++){//强制转换(必然)
        register int pos=locate(prime[j]*i);
        if (h[pos]!=prime[j]*i){
            h[pos]=prime[j]*i;
            ans=(ans+prime[j]*i)%376544743;//统计答案
        }
        if (i%prime[j]==0) break;
    }
    print(ans);
    return 0;
}
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值