靈淵(seals.cpp/c/pas)

標籤：style   blog   http   color   os   art   

題意：p(m)的值為m的正因數個數（包括1和m本身）。

　　　求滿足p(x)=n的x的最小值。

對於任意正整數n，有n=p1^a1 * p2^a2 * p3^a3 * …… * pn^an;（pi為質數)
n的因數個數(a1+1)*(a2+1)*(a3+1)*……*(an+1);

舉個例子，8=2*2*2；
可以這樣 2^1 * 3^1 * 5^1 =30;
因為 8也可分解為這種形式： 8=2*4；2^3 * 3^1 =24;

 正確答案為 24 。

假設 有一個數 n=b1*b2*b3*……*bi*……*bn；
如果把它分解成 n=b1*b2*b3*……*(bj*bi)*……*b[i-1]*b[i+1]*……*bn，計算出來的結果比上式更小，
那一定滿足這個條件：
p[n-i+1]^bi * p[n-j+1]^(bj)  >  p[n-j+1]^(bj*bi)   
約掉p[n-j+1]^(bj) 得：
p[n-i+1]^bi  >  p[n-j+1]^(bj*(bi - 1)) ;   
只要滿足這個條件 就可一把bj賦為 bj*bi ，把bi賦為0，表示清除；
之後就可以了。

 

代碼如下

#include<iostream>#include<cmath>#include<cstdio>using namespace std;int prime[]={0,2,3,5,7,11,13,17,19,23,29,31,37,41,43,47,53,59,61,67,71,73,79,89,97,101,103,107,109,113,127,131,137,139};int n,i,temp,a[1000]={0},j,ans[10000][2],c[10000];double pro(int x,int y){    double ans=x,t=x;    int c[100],i;c[0]=0;    while (y!=0){c[++c[0]]=y%2;y/=2;}    for (i=c[0]-1;i>=1;i--){        ans*=ans;        if(c[i]==1)ans*=t;    }    return ans;}void clean(){    int i,j;    for(i=1;i<=a[0];i++)        for(j=a[0];j>i;j--)        if(a[i]*a[j]!=0)            if (pro(prime[a[0]-i+1],a[i]-1)>pro(prime[a[0]-j+1],a[j]*(a[i]-1)))            {               a[j]*=a[i]; a[i]=0; break;            }}int main(){    freopen("seals.in","r",stdin);    freopen("seals.out","w",stdout);    scanf("%d",&n);    for    (i=2;i<=sqrt(n);i++)     while (n%i==0) { a[++a[0]]=i; n/=i;}     if (n>1)a[++a[0]]=n;      clean();      int tot=0,t;       for (i=1;i<=a[0];i++)    if (a[a[0]-i+1]>0){         temp=a[a[0]-i+1]-1;           ans[++tot][1]=prime[i];           ans[tot][2]=temp;    }    for (i=1;i<tot;i++) printf("%d^%d*",ans[i][1],ans[i][2]);     printf("%d^%d",ans[tot][1],ans[tot][2]); }

View Code