Codeforces 851D Arpa and a list of numbers【思维+前缀和】

D. Arpa and a list of numbers
time limit per test
2 seconds
memory limit per test
256 megabytes
input
standard input
output
standard output

Arpa has found a list containing n numbers. He calls a list bad if and only if it is not empty and gcd (see notes section for more information) of numbers in the list is 1.

Arpa can perform two types of operations:

  • Choose a number and delete it with cost x.
  • Choose a number and increase it by 1 with cost y.

Arpa can apply these operations to as many numbers as he wishes, and he is allowed to apply the second operation arbitrarily many times on the same number.

Help Arpa to find the minimum possible cost to make the list good.

Input

First line contains three integers nx and y (1 ≤ n ≤ 5·1051 ≤ x, y ≤ 109) — the number of elements in the list and the integers x and y.

Second line contains n integers a1, a2, ..., an (1 ≤ ai ≤ 106) — the elements of the list.

Output

Print a single integer: the minimum possible cost to make the list good.

Examples
input
4 23 17
1 17 17 16
output
40
input
10 6 2
100 49 71 73 66 96 8 60 41 63
output
10
Note

In example, number 1 must be deleted (with cost 23) and number 16 must increased by 1 (with cost 17).

gcd (greatest common divisor) of a set of numbers is the maximum integer that divides all integers in the set. Read more about gcd here.


题目大意:


给出N个数,删除一个数的花费是X,改变一个数从num变成num+1的花费是Y,问将整个序列的Gcd改成不是1的最小花费。


思路:


我们知道,Gcd的结果是素数是最好的,所以我们首先筛出素数,1e6内的素数个数约为8w(8e4)个。


然后考虑对应数字是删除还是增添。很显然我们这里可以Dp去做,但是时间复杂度很爆炸,所以考虑X和Y两种操作之间的关系。

①不难想到,如果X<Y的话,我们就不用进行Y操作了。

②所以我们设定move=X/Y表示我们最多对一个数进行增添数字的次数,如果超过了这个次数,我们不如将其删除来的划算。


接下来考虑如何计算结果。首先我们肯定要O(8e4)去枚举答案,接下来考虑,既然我们已经知道了move这个信息,那么对于当前枚举的素数p,我们有一个区间【p,2p】的话,我们肯定能够通过move这个信息找到区间的一个分界点,使得【p~分界点】的数我们都将其删除,使得【分界点~2p】的数字都增添变成2p.


那么这里我们只要枚举出来一个素数p,然后每一次计算之后将区间移动p个长度即可。整体时间复杂度O(Len*LogLen)这里Len表示素数的个数。


那么过程维护一下即可,具体计算过程参考代码。


Ac代码:

#include<stdio.h>
#include<string.h>
#include<iostream>
using namespace std;
#define ll long long int
ll have[3005000];
ll sum[3005000];
int main()
{
    ll n,x,y;
    while(~scanf("%lld%lld%lld",&n,&x,&y))
    {
        memset(have,0,sizeof(have));
        memset(sum,0,sizeof(sum));
        for(ll i=1;i<=n;i++)
        {
            ll tmp;scanf("%lld",&tmp);
            have[tmp]++;sum[tmp]+=tmp;
        }
        ll ans=1000000000000000000;
        for(ll i=1;i<=2000000;i++)have[i]+=have[i-1],sum[i]+=sum[i-1];
        ll move=(x/y);
        for(ll i=2;i<=1000000;i++)
        {
            ll now=0;
            for(ll j=i;j<=2000000;j+=i)
            {
                ll fen=max(j-i+1,j-move);
                ll Inum=have[j]-have[fen-1];
                ll Isum=sum[j]-sum[fen-1];
                ll Dnum=have[fen-1]-have[j-i];
                now+=Dnum*x;
                now+=y*(j*Inum-Isum);
            }
            ans=min(ans,now);
        }
        printf("%lld\n",ans);
    }
}










评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值