1. 程式人生 > >P1128 [HNOI2001]求正整數

P1128 [HNOI2001]求正整數

傳送門

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

//minamoto
#include<bits/stdc++.h>
#define fp(i,a,b) for(register int i=a,I=b+1;i<I;++i)
#define fd(i,a,b) for(register int i=a,I=b-1;i>I;--i)
using namespace std;
const double min(const double &x,const double &y){return x<y?x:y;}
const int p[]={
   2,  3,  5,  7, 11,
  13, 17, 19, 23, 29,
  31, 37, 41, 43, 47,
  53, 59, 61, 67, 71
};double Log[25],f[505][25];int d[505],n,tot,len,A[100005],pos;
void mul(int x){
    int add=0;
    fp(i,1,len)A[i]=A[i]*x+add,add=A[i]/10,A[i]%=10;
    while(add)A[++len]=add%10,add/=10;
}
int main(){
//  freopen("testdata.in","r",stdin);
    scanf("%d",&n);fp(i,1,n)if(n%i==0)d[++tot]=i;
    fp(i,0,19)f[0][i]=0;fp(i,0,19)Log[i]=log(p[i]);
    fp(i,2,tot){
        fp(k,0,19)f[i][k]=1e9;
        fp(j,1,i-1)if(d[i]%d[j]==0){
            int t=d[i]/d[j];
            fp(k,1,19)f[i][k]=min(f[i][k],f[j][k-1]+Log[k-1]*(t-1));
        }
    }A[1]=len=1;
    fp(i,0,19)if(f[tot][i]<f[tot][pos])pos=i;
    for(register int i=tot,nxt;i>1;i=nxt,--pos){
        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);
//      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);
        fp(k,1,d[i]/d[nxt]-1)mul(p[pos-1]);
    }
    while(len--)printf("%d",A[len+1]);return 0;
}