BZOJ3233 [Ahoi2013]找硬币(线性筛+dp)
2015-06-26 03:04
337 查看
【题解】
本蒻一直在想二维dp,看了题解才发现竟然一维就可以
设f[i]为最大面值为i时,买下所有兔纸花费的最小硬币数
f[i] = min{ f[j] - sigma(a[k]/i*(i/j-1)) } , j|i,其中,j为次大面值,这个方程考虑的是选了i能减小多少j的使用
注意,如果硬币种类很多,是不影响最优答案的(不用就行了) -----------> 重要的性质
所以,可以在枚举j的这一步下手优化:
规定i/j必须为质数,否则可以在它们之间加一个面值k,使 i/k,k/j都是质数,不影响答案
可以在线性筛的过程中求出数组min[i]:i的最小质因子,用它实现:仅枚举i的质因子作为i/j的值
n的质因子个数是log(n)级别的,复杂度:O( n*Max*log(Max) )
【代码】
本蒻一直在想二维dp,看了题解才发现竟然一维就可以
设f[i]为最大面值为i时,买下所有兔纸花费的最小硬币数
f[i] = min{ f[j] - sigma(a[k]/i*(i/j-1)) } , j|i,其中,j为次大面值,这个方程考虑的是选了i能减小多少j的使用
注意,如果硬币种类很多,是不影响最优答案的(不用就行了) -----------> 重要的性质
所以,可以在枚举j的这一步下手优化:
规定i/j必须为质数,否则可以在它们之间加一个面值k,使 i/k,k/j都是质数,不影响答案
可以在线性筛的过程中求出数组min[i]:i的最小质因子,用它实现:仅枚举i的质因子作为i/j的值
n的质因子个数是log(n)级别的,复杂度:O( n*Max*log(Max) )
【代码】
#include<stdio.h> #include<stdlib.h> int a[100005],f[100005],pri[100005],min[100005]; int main() { int n,i,j,k,Ma=0,cnt=0,t,ans; scanf("%d",&n); for(i=1;i<=n;i++) { scanf("%d",&a[i]); f[1]+=a[i]; if(Ma<a[i]) Ma=a[i]; } for(i=2;i<=Ma;i++) { if(min[i]==0) { pri[++cnt]=i; min[i]=i; } for(j=1;j<=cnt&&pri[j]*i<=Ma;j++) { min[pri[j]*i]=pri[j]; if(i%pri[j]==0) break; } } for(i=2;i<=Ma;i++) f[i]=f[1]; ans=f[1]; for(i=2;i<=Ma;i++) { j=i; while(j>1) { t=f[i/min[j]]; for(k=1;k<=n;k++) t-=a[k]/i*(min[j]-1); if(f[i]>t) f[i]=t; while(min[j]==min[j/min[j]]) j/=min[j]; j/=min[j]; } if(ans>f[i]) ans=f[i]; } printf("%d",ans); return 0; }
相关文章推荐
- how convert large HEX string to binary array ?
- js判断是否为数组的函数: isArray()
- redis info参数详解
- 【温故知新】c#异步编程模型(APM)--使用委托进行异步编程
- [windows通用应用开发] win8.1中制作类似win10菜单栏效果
- CentOS安装部署Cacti
- 并发性测试AB(apache benchmark)
- More Effective C++ 条款4 非必要不提供 default constructor
- 如何自己编写一个easyui插件
- 个人applicationthread解惑
- linux kbuild module
- shell 中如何判断前一个命令是否执行成功
- shell 中如何判断前一个命令是否执行成功
- 【笔记】iOS开发基础笔记<二>
- python3.4学习笔记(七) 学习网站博客推荐
- Windows Phone 二十一、联系人存储
- 查看当前发行版可以使用的shell
- 查看当前发行版可以使用的shell
- bestcoder #45 分类: bestcoder ...
- bestcoder #45