test20181018 B君的第三题
题意
B 君的第三题(shenyang)
题目描述
客似云来,万里无云
B 君得到了一个数组\(\{a_1,a_2,\dots,a_n\}\)。
B 君想通过修改让数组中个每对数都互质。
每次使一个数+1 或者-1 的代价是1。
不能将\(a_i\) 修改为0 或者负数。
问至少多少代价才可以让所有数两两互质。
输入格式
第一行一个整数n。
第二行n 个整数\(a_i\),表示数组初始值。
输出格式
一行一个数表示答案。
样例输入
5
2 4 6 8 10
样例输出
4
样例解释
修改为2, 3, 5, 7, 11。
数据规模与约定
对于100% 的数据,满足\(1 \leq n \leq 100, 1 \leq a_i \leq 30\)。
对于30% 的数据,满足\(1 \leq n \leq 4\)。
对于另30% 的数据,满足\(1 \leq a_i \leq 10\)。
分析
因为\(a_i \leq 30\),1跟所有数互质,所以往下最多改到1,往上最多改到59,所以每个数所含质因数是固定的。
考虑dp,用\(f(i,s)\)表示前i个数改动至只含s集合中的质因数且两两互质所需最小代价,转移方程为:
枚举从原来的含k的集合转移过来,其中k与s没有交集 \\
f(i,k∨s)=\min_{x=1}^{59}f(i-1,k)+|a_i-x|
\]
时间复杂度\(O(n \cdot 59 \cdot 2^{17})\),上界十分不紧。
代码
#include<cstdlib>
#include<cstdio>
#include<cmath>
#include<cstring>
#include<iostream>
#include<string>
#include<vector>
#include<list>
#include<deque>
#include<stack>
#include<queue>
#include<map>
#include<set>
#include<bitset>
#include<algorithm>
#include<complex>
#define rg register
#define il inline
#define co const
#pragma GCC optimize ("O0")
using namespace std;
template<class T> il T read()
{
T data=0;
int w=1;
char ch=getchar();
while(!isdigit(ch))
{
if(ch=='-')
w=-1;
ch=getchar();
}
while(isdigit(ch))
data=10*data+ch-'0',ch=getchar();
return data*w;
}
template<class T> il T read(T&x)
{
return x=read<T>();
}
typedef long long ll;
const int INF=0x7fffffff;
int p[17]={2,3,5,7,11,13,17,19,23,29,31,37,41,43,47,53,59};
il int bit(rg int x)
{
rg int s=0;
for(rg int i=0;i<17;++i)
if(x%p[i]==0)
s|=(1<<i);
return s;
}
int b[60];
const int MAXN=107;
int f[MAXN][1<<17];
int main()
{
freopen("shenyang.in","r",stdin);
freopen("shenyang.out","w",stdout);
memset(f,0x3f,sizeof f);
f[0][0]=0;
rg int n=read<int>();
for(rg int i=1;i<=59;++i)
b[i]=bit(i);
for(rg int i=1;i<=n;++i)
{
rg int a=read<int>();
for(rg int j=1;j<=59;++j)
{
rg int ub=((1<<17)-1)^b[j];
for(rg int k=ub;;--k&=ub)
{
f[i][k|b[j]]=min(f[i][k|b[j]],f[i-1][k]+abs(a-j));
if(k==0)
break;
}
}
}
rg int ans=INF,up=(1<<17)-1;
for(rg int i=0;i<=up;++i)
ans=min(ans,f[n][i]);
printf("%d\n",ans);
// fclose(stdin);
// fclose(stdout);
return 0;
}
test20181018 B君的第三题的更多相关文章
- test20181016 B君的第三题
题意 B 君的第三题(haskell) 题目描述 大学四年,我为什么,为什么不好好读书,没找到和你一样的工作. B 君某天看到了这样一个题,勾起了无穷的回忆. 输入\(n, k\) 和一棵\(n\) ...
- test20181015 B 君的第三题
题意 B 君的第三题(zhengzhou) 题目描述 让你在战争和耻辱中做一块选择,你选择耻辱,可你将来还得进行战争. 在平面上有n 个整点(横纵坐标都是整数) B 君想找到一个整点,使得这个点,到所 ...
- test20181019 B君的第三题
题意 B 君的第三题(urumqi) 题目描述 风雨如晦,鸡鸣不已. B 君最近在研究自己的学长都在做什么工作,每个学长属于一个公司. B 君会获得一些信息,比如x 和y 在相同公司,x 和y 在不同 ...
- noiac132 B君的第三题 (树形dp)
传送门 本来想用点分治做,结果root又求不对 算的时候还算错了 我好菜啊 结果szr大佬告诉我是树形dp 我好菜啊!! 我们有$\lceil \frac{x}{k} \rceil = \frac{x ...
- shanquan2的两年三题系列
好像只有2个月就退役啦 不管了,先说一下哪三题:多点求值.lcm.替罪羊树(bzoj3065) [upd0]2016.3.29 多点求值A掉啦,myy卡常数sxbk(不是说好的是shanquan2出的 ...
- Java-集合-第三题 有如下Student 对象, private String name; private int age; private int score; private String classNum; 其中,classNum 表示学生的班号,例如“class05”。 有如下List List list = new ArrayList(); l
第三题 有如下Student 对象, private String name; private int age; private int score; private String classNum; ...
- NOIP 2008提高组第三题题解by rLq
啊啊啊啊啊啊今天已经星期三了吗 那么,来一波题解吧 本题地址http://www.luogu.org/problem/show?pid=1006 传纸条 题目描述 小渊和小轩是好朋友也是同班同学,他们 ...
- NOIP2008提高组(前三题) -SilverN
此处为前三题,第四题将单独发布 火柴棒等式 题目描述 给你n根火柴棍,你可以拼出多少个形如“A+B=C”的等式?等式中的A.B.C是用火柴棍拼出的整数(若该数非零,则最高位不能是0).用火柴棍拼数字0 ...
- 《学习OpenCV》练习题第四章第三题b
#include <highgui.h> #include <cv.h> #include "opencv_libs.h" /* *<学习OpenCV ...
随机推荐
- A题:Common Substrings(KMP应用)
原题链接 注意:2号和3号get_next()函数中next[i]赋值时的区别,一个是0,一个是1,且不能互换 #include<cstdio> #include<cstring&g ...
- python实现指定目录下批量文件的单词计数:并发版本
在 文章 <python实现指定目录下批量文件的单词计数:串行版本>中, 总体思路是: A. 一次性获取指定目录下的所有符合条件的文件 -> B. 一次性获取所有文件的所有文件行 - ...
- 使用 Laravel 数据填充功能生成中文测试数据
今晚……不对,是昨晚,折腾一个的小项目,发现自动填充的中文数据显示起来总不太美观,于是开始琢磨如何填充中文数据进行测试. 然而一番搜索后惊奇的发现,官方.以及一些非官方的文档均未提及这一功能.期间看到 ...
- Spring整合Quartz定时发送邮件
功能描述:刚开始接触Quartz,试着用Quartz整合spring实现每隔一分钟发送一封邮件连续发送10次 核心jar: 邮件发送:commons-email-1.2.jar mail.jar(必须 ...
- MAC nginx代理设置
问题: 10.154.156.83:10081私服不存在了.但是不能改.用nginx代理至maven.xx.cn 增加换回地址: sudo ifconfig lo0 add 10.154.156.83 ...
- 《Effective Java 2nd》第2章 创建和销毁对象
目录 第1条:考虑使用静态工厂方法代替构造器 第2条:遇到多个构造器参数时考虑用构建器 第3条:用私有构造器或者枚举类型强化Singleton属性 第4条:通过私有构造器强化不可实例化的能力 第5条: ...
- BootStrap iCheck插件全选与获取value值的解决方法
这篇文章主要介绍了BootStrap iCheck插件全选与获取value值的解决方法,解决方法其实很简单,下面小编给大家分享下这方面的知识 在使用jQuery iCheck 插件的时候遇到了一个问题 ...
- .Net t图片生成水印
借鉴于博客园园友的方法,大神神风(https://www.cnblogs.com/tandyshen/archive/2012/04/14/picwater.html) ,很实用的一种 自定义水印方法 ...
- POJ 2185 Milking Grid(KMP最小循环节)
http://poj.org/problem?id=2185 题意: 给出一个r行c列的字符矩阵,求最小的覆盖矩阵可以将原矩阵覆盖,覆盖矩阵不必全用完. 思路: 我对于字符串的最小循环节是这么理解的: ...
- c++ 列表删除元素(erase)
#include <list> #include <iostream> #include <iterator> using namespace std; int m ...