D. Arpa and a list of numbers

time limit per test   2 seconds

memory limit per test     256 megabytes

Arpa has found alist containingn numbers. Hecalls 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 performtwo 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 applythese operations to as many numbers as he wishes, and he is allowed to applythe second operation arbitrarily many times on the same number.

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

Input

First linecontains three integersn,x and y (1 ≤ n ≤ 5·105,1 ≤ x, y ≤ 109) — the number of elements in the list and the integers x and y.

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

Output

Print a single integer: the minimumpossible 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,number1 must be deleted (with cost 23) and number 16 must increased by 1 (with cost 17).

Agcd (greatest common divisor) of a set of numbers is the maximum integer thatdivides all integers in the set. Read more about gcdhere.

【题意】

给出n个数,以及x,y。现在你可以对这n个数进行两种操作。

  1. 把任意一个数删除,花费为x。
  2. 把任意一个数加一,花费为y。

现在要求操作后所有数gcd不为1(不互质),求最小花费。

【思路】

显然可以想到去枚举gcd,考虑到时间上的优化,我们可以枚举素数,因为每个不为1的数一定是一个素数或者一个素数的倍数。(事实上直接暴力枚举gcd也能卡过)

素数用素筛预处理一下即可。

假设我们枚举的素数为prime。那么我们需要把所有数变为这个素数或其倍数或者将其中一个(或几个)删除。

显然如果我们要改变某个数的值,我们一定是增加到离它最近的是prime的倍数的那个数。所以我们应该提前预处理一下前缀和sum[i]、num[i],sum[i]表示小于等于i的数的和,num[i]表示小于等于i的数的个数,然后在区间内进行操作。

假设现在枚举的区间为【 k * prime[i] , (k+1) * prime[i]),那么我们先应该在区间内找一个分界点limit,一个数小于等于limit时删除它花费更少,否则把它一直加到k * prime[i]的花费更少。

limit可由x,y的大小得到。

limit=max( (j+1) * prime[i]-rate-1 , j * prime[i]);

然后就对区间里的数根据limit分成的两部分分别计算。

对于删除的那部分,计算方法为删除个数*x,即:

num[ limit ]-num[ j * prime[i]] * x

而对于加上的那部分,计算方法为(数的个数*加上后所得到的那个素数的倍数-原来的数的和) * y,即:

((num[(j+1) * prime[i]]-num[limit]) * ((j+1) * prime[i])-(sum[(j+1) * prime[i]]-sum[limit])) * y

然后不断更新最小值即可。

#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;
#define mst(a,b) memset((a),(b),sizeof(a))
#define rush() int T;scanf("%d",&T);while(T--) typedef long long ll;
const int maxn = *;
const ll mod = 1e9+;
const int INF = 0x3f3f3f3f;
const double eps = 1e-; int prime[maxn]= {};
int num_prime=;
bool isprime[maxn]= {,};
ll num[maxn],sum[maxn]; void init()
{
for(int i=;i<maxn;i++)
{
if(!isprime[i]) prime[num_prime++]=i;
for(int j=;j<num_prime&&i*prime[j]<maxn;j++)
{
isprime[i*prime[j]]=;
if(i%prime[j]==) break;
}
}
} int main()
{
int o;
int n,x,y;
init();
while(~scanf("%d%d%d",&n,&x,&y))
{
mst(num,);
mst(sum,);
int Max=;
for(int i=;i<n;i++)
{
scanf("%d",&o);
num[o]++;
sum[o]+=o;
Max=max(Max,o);
}
for(int i=;i<=Max*;i++) //预处理前缀和
{
num[i]+=num[i-];
sum[i]+=sum[i-];
}
int rate=x/y;
ll ans=1e18;
for(int i=;i<num_prime&&prime[i-]<=Max;i++)
{
ll cnt=;
for(int j=;j*prime[i]<=Max;j++)
{
int limit=max((j+)*prime[i]-rate-,j*prime[i]);
ll num1=num[limit]-num[j*prime[i]]; //删除
ll sum2=sum[(j+)*prime[i]]-sum[limit]; //加上
ll num2=num[(j+)*prime[i]]-num[limit];
cnt+=num1*x;
cnt+=(num2*((j+)*prime[i])-sum2)*y;
if(cnt>ans) break; //优化
}
ans=min(ans,cnt);
}
printf("%I64d\n",ans);
}
return ;
}

Codeforces Codeforces Round #432 (Div. 2 D ) Arpa and a list of numbers的更多相关文章

  1. Codeforces Round #432 (Div. 1) B. Arpa and a list of numbers

    qtmd的复习pat,老子不想看了,还不如练几道cf 这题首先可以很容易想到讨论最后的共因子为素数 这个素数太多了,1-1e6之间的素数 复杂度爆炸 所以使用了前缀和,对于每个素数k的每个小区间 (k ...

  2. Codeforces Round #432 Div. 1 C. Arpa and a game with Mojtaba

    首先容易想到,每种素数是独立的,相互sg就行了 对于一种素数来说,按照的朴素的mex没法做... 所以题解的简化就是数位化 多个数同时含有的满参数因子由于在博弈中一同变化的,让他们等于相当于,那么这样 ...

  3. Codeforces Beta Round #80 (Div. 2 Only)【ABCD】

    Codeforces Beta Round #80 (Div. 2 Only) A Blackjack1 题意 一共52张扑克,A代表1或者11,2-10表示自己的数字,其他都表示10 现在你已经有一 ...

  4. Codeforces Beta Round #83 (Div. 1 Only)题解【ABCD】

    Codeforces Beta Round #83 (Div. 1 Only) A. Dorm Water Supply 题意 给你一个n点m边的图,保证每个点的入度和出度最多为1 如果这个点入度为0 ...

  5. Codeforces Beta Round #79 (Div. 2 Only)

    Codeforces Beta Round #79 (Div. 2 Only) http://codeforces.com/contest/102 A #include<bits/stdc++. ...

  6. Codeforces Beta Round #77 (Div. 2 Only)

    Codeforces Beta Round #77 (Div. 2 Only) http://codeforces.com/contest/96 A #include<bits/stdc++.h ...

  7. Codeforces Beta Round #76 (Div. 2 Only)

    Codeforces Beta Round #76 (Div. 2 Only) http://codeforces.com/contest/94 A #include<bits/stdc++.h ...

  8. Codeforces Beta Round #75 (Div. 2 Only)

    Codeforces Beta Round #75 (Div. 2 Only) http://codeforces.com/contest/92 A #include<iostream> ...

  9. Codeforces Beta Round #74 (Div. 2 Only)

    Codeforces Beta Round #74 (Div. 2 Only) http://codeforces.com/contest/90 A #include<iostream> ...

随机推荐

  1. linux常用命令(15)whereis命令

    whereis命令只能用于程序名的搜索,而且只搜索二进制文件(参数-b).man说明文件(参数-m)和源代码文件(参数-s).如果省略参数,则返回所有信息.和find相比,whereis查找的速度非常 ...

  2. JSP 传统标签extends TagSupport

    1.控制JSP页面某一部分内容是否执行 public int doStartTag() 返回EVAL_BODY_INCLUDE,执行 返回SKIP_BODY,不执行 2.控制整个JSP页面是否执行 p ...

  3. BOM Summary P268-P269

    The Browser Object Model(BOM) is based on the window object, which represents the browser window and ...

  4. SqlServer:SqlServer(存储过程动态表查询(取消返回值),事务处理,批量还原sqlserver备份,强制删除被占用的数据库)

    1.存储过程动态表查询 USE [NETWORKING_AUDIT] GO /****** Object: StoredProcedure [dbo].[impConfigInfo] Script D ...

  5. CentOS7使用阿里云源安装Docker

    安装步骤 1.删除已安装的Docker # Uninstall installed docker sudo yum remove docker \ docker-client \ docker-cli ...

  6. C学习笔记-内存管理

    作用域 一个C语言变量的作用域可以是代码块 作用域,函数作用域或者文件作用域 代码块是{}之间的一段代码 同一个代码块不可以有重名变量 auto自动变量 一般情况下代码块内部定义的变量都是自动变量 也 ...

  7. 单例模式(一)static、final和单例模式

    static 那天我朋友问了我个问题,static和单例模式有什么区别,所以我觉得static可以讲一下 他的问题是,把对象弄成static是不是就不变了 显然,这是还没弄清楚引用和对象的区别 其实存 ...

  8. PTA(Basic Level)1029.旧键盘

    旧键盘上坏了几个键,于是在敲一段文字的时候,对应的字符就不会出现.现在给出应该输入的一段文字.以及实际被输入的文字,请你列出肯定坏掉的那些键. 输入格式: 输入在 2 行中分别给出应该输入的文字.以及 ...

  9. 小菜鸟之Cisco

    Switch>enable// Switch#configure// Switch#show vlan//展示vlan接口 Switch(config)#enable password 1234 ...

  10. 线性基求交(线段树)--牛客第四场(xor)

    题意: 给你n个基,q个询问,每个询问问你能不能 l~r 的所有基都能表示 x . 思路: 建一颗线性基的线段树,up就是求交的过程,按照线段树区间查询的方法进行check就可以了. #define ...