传送门

rqy是我们的红太阳没有它我们就会死

可以考虑dp,设\(dp[i][j]\)表示只包含前\(j\)个质数的数中,因子个数为\(i\)的数的最小值是多少,那么有转移方程

\[f[i][j]=min(f[i/k][j-1]\times p_j^{k-1})
\]

这玩意儿肯定是不能高精dp的……于是看到乘法就想到对数……根据对数的基本定理,有

\[log\ n=\sum c_i\ log\ p_i
\]

那么我们可以改一下转移

\[f[i][j]=min(f[i/k][j-1]+ (k-1)log\ p_j)
\]

然后算出最后的答案之后倒着找一下转移的方向高精乘回去就是了

因为没有写eps于是调了一个小时都找不出错我再也不偷懒不写eps了QAQ

  1. //minamoto
  2. #include<bits/stdc++.h>
  3. #define fp(i,a,b) for(register int i=a,I=b+1;i<I;++i)
  4. #define fd(i,a,b) for(register int i=a,I=b-1;i>I;--i)
  5. using namespace std;
  6. const double min(const double &x,const double &y){return x<y?x:y;}
  7. const int p[]={
  8. 2, 3, 5, 7, 11,
  9. 13, 17, 19, 23, 29,
  10. 31, 37, 41, 43, 47,
  11. 53, 59, 61, 67, 71
  12. };double Log[25],f[505][25];int d[505],n,tot,len,A[100005],pos;
  13. void mul(int x){
  14. int add=0;
  15. fp(i,1,len)A[i]=A[i]*x+add,add=A[i]/10,A[i]%=10;
  16. while(add)A[++len]=add%10,add/=10;
  17. }
  18. int main(){
  19. // freopen("testdata.in","r",stdin);
  20. scanf("%d",&n);fp(i,1,n)if(n%i==0)d[++tot]=i;
  21. fp(i,0,19)f[0][i]=0;fp(i,0,19)Log[i]=log(p[i]);
  22. fp(i,2,tot){
  23. fp(k,0,19)f[i][k]=1e9;
  24. fp(j,1,i-1)if(d[i]%d[j]==0){
  25. int t=d[i]/d[j];
  26. fp(k,1,19)f[i][k]=min(f[i][k],f[j][k-1]+Log[k-1]*(t-1));
  27. }
  28. }A[1]=len=1;
  29. fp(i,0,19)if(f[tot][i]<f[tot][pos])pos=i;
  30. for(register int i=tot,nxt;i>1;i=nxt,--pos){
  31. for(nxt=1;d[i]%d[nxt]||f[i][pos]<f[nxt][pos-1]+Log[pos-1]*(d[i]/d[nxt]-1)-1e-5;++nxt);
  32. // for(nxt=1;d[i]%d[nxt]||f[i][pos]<f[nxt][pos-1]+Log[pos-1]*(d[i]/d[nxt]-1)-1e-5;++nxt);
  33. fp(k,1,d[i]/d[nxt]-1)mul(p[pos-1]);
  34. }
  35. while(len--)printf("%d",A[len+1]);return 0;
  36. }

P1128 [HNOI2001]求正整数的更多相关文章

  1. luogu P1128 [HNOI2001]求正整数 dp 高精度

    LINK:求正整数 比较难的高精度. 容易想到贪心不过这个贪心的策略大多都能找到反例. 考虑dp. f[i][j]表示前i个质数此时n的值为j的最小的答案. 利用高精度dp不太现实.就算上FFT也会T ...

  2. 高精度+搜索+质数 BZOJ1225 [HNOI2001] 求正整数

    // 高精度+搜索+质数 BZOJ1225 [HNOI2001] 求正整数 // 思路: // http://blog.csdn.net/huzecong/article/details/847868 ...

  3. BZOJ 1225: [HNOI2001] 求正整数( dfs + 高精度 )

    15 < log250000 < 16, 所以不会选超过16个质数, 然后暴力去跑dfs, 高精度计算最后答案.. ------------------------------------ ...

  4. bzoj1225 [HNOI2001] 求正整数

    1225: [HNOI2001] 求正整数 Time Limit: 10 Sec  Memory Limit: 162 MBSubmit: 762  Solved: 313[Submit][Statu ...

  5. 【BZOJ】1225: [HNOI2001] 求正整数

    http://www.lydsy.com/JudgeOnline/problem.php?id=1225 题意:给一个数n,求一个最小的有n个约数的正整数.(n<=50000) #include ...

  6. [HNOI2001]求正整数

    题目描述 对于任意输入的正整数n,请编程求出具有n个不同因子的最小正整数m. 例如:n=4,则m=6,因为6有4个不同整数因子1,2,3,6:而且是最小的有4个因子的整数. 输入输出格式 输入格式: ...

  7. BZOJ 1225: [HNOI2001] 求正整数 高精度+搜索+质数

    题意:给定n求,有n个因子的最小正整数. 题解:水题,zcr都会,我就不说什么了. 因数个数球求法应该知道,将m分解质因数,然后发现 a1^p1*a2^p2....an^pn这样一个式子, (1+p1 ...

  8. [HNOI2001] 求正整数 - 背包dp,数论

    对于任意输入的正整数n,请编程求出具有n个不同因子的最小正整数m. Solution (乍一看很简单却搞了好久?我真是太菜了) 根据因子个数计算公式 若 \(m = \prod p_i^{q_i}\) ...

  9. 求正整数n所有可能的和式的组合(如;4=1+1+1+1、1+1+2、1+3、2+1+1、2+2

    作者:张小二 nyoj90 ,可以使用递归的方式直接计算个数,也可以通过把满足的个数求出来计数,因为在juLy博客上看到整数划分,所以重写了这个代码,就是列出所m的可能性,提交后正确.acmer的入门 ...

随机推荐

  1. 动态规划之最长公共子序列(LCS)

             在字符串S中按照其先后顺序依次取出若干个字符,并讲它们排列成一个新的字符串,这个字符串就被称为原字符串的子串          有两个字符串S1和S2,求一个最长公共子串,即求字符串 ...

  2. uva1366/LA3530

    简单的矩阵上的dp 从左上角扫到右下角,dp[i][j][0] 代表在i,j处选A ,dp[i][j][1] 代表在i,j处选B. dp[i][j][0]=max(dp[i-1][j][1],dp[i ...

  3. 一起去打CS(codevs 5059)

    题目描述 Description 早就和lyk约好了去打cs,一直没找着时间,终于今天我家没人,他家也没人,总算可以出去了.但是偏偏天公不作美,某某人非要留那么多题要做.没办法只能尽快做完然后抓紧时间 ...

  4. 【BZOJ3697】采药人的路径(点分治)

    题意:采药人的药田是一个树状结构,每条路径上都种植着同种药材.采药人以自己对药材独到的见解,对每种药材进行了分类.大致分为两类,一种是阴性的,一种是阳性的.采药人每天都要进行采药活动.他选择的路径是很 ...

  5. codevs3728 联合权值

    题目描述 Description 输入描述 Input Description 输出描述 Output Description 样例输入 Sample Input 样例输出 Sample Output ...

  6. Nginx配置文件nginx.conf 详解

    #定义Nginx运行的用户和用户组 user www www;   #nginx进程数,建议设置为等于CPU总核心数. worker_processes 8;   #全局错误日志定义类型,[ debu ...

  7. P1160 队列安排 洛谷

    https://www.luogu.org/problem/show?pid=1160 题目描述 一个学校里老师要将班上N个同学排成一列,同学被编号为1-N,他采取如下的方法: 1.先将1号同学安排进 ...

  8. Swap Nodes in Pairs(链表操作)

    Given a linked list, swap every two adjacent nodes and return its head. For example,Given 1->2-&g ...

  9. Ubuntu 16.04设置rc.local开机启动命令/脚本的方法(通过update-rc.d管理Ubuntu开机启动程序/服务)

    注意:rc.local脚本里面启动的用户默认为root权限. 一.rc.local脚本 rc.local脚本是一个Ubuntu开机后会自动执行的脚本,我们可以在该脚本内添加命令行指令.该脚本位于/et ...

  10. js 合并对象

    对象的合并 需求:设有对象 o1 ,o2,需要得到对象 o3 var o1 = { a:'a' }, o2 = { b:'b' }; // 则 var o3 = { a:'a', b:'b' } 方法 ...