POJ1338

POJ1338题 丑数

题目描述：素数因子仅为2,3,5的数为丑数。前10个丑数分别是：1,2,3,4,5,6,8,9,10,12。

现在给出n，求第n个丑数。

输入：n<=1500          //第1500个丑数大概为8亿多

输出：第n个丑数

         由分析可得第n个丑数一定是由前n-1个丑数或*2或*3或*5得来的.

设A[P2]*2>A[n-1]且A[P2-1]*2<=A[n-1]。设A[P3]*3>A[n-1]且A[P3-1]*3<=A[n-1]。设A[P5]*5>A[n-1]且A[P5-1]*5<=A[n-1]。则P2所指的丑数正好是前n-1个丑数中乘以2之后刚刚大于A[n-1]的那个。P3,P5所指含义类似。如下例子：

1, 2, 3, 4, 5, 6, 8, 9, 10, 12 当求出前10个丑数（为12）了，要求第11个丑数，此时按语义P2指向8 ，P3指向 5，p5指向3.则第11个丑数要不就是8*2=16，要不就是5*3=15，要不就是3*5=15。

看这三个数哪个最小了，可以看出15最小，则第11个丑数为15.且p3++，p5++；因为当求出15之后p3所指的5*3<=15了，所以p3要加1.

         维护指针p2，p3，p5可以求得第1500个丑数。

代码如下：

#include<cstdio>
using namespacestd;
 
const int maxn =1500+10;
int a[maxn];
 
int main()
{
    a[1]=1;
    int p2=1,p3=1,p5=1;
    for(int i=2; i<=1500; i++)
    {
        int t2=a[p2]*2,t3=a[p3]*3,t5=a[p5]*5;
        if(t2<t3&&t2<t5)
        {
            a[i]=t2;
            p2++;
        }
        else if(t3<t2&&t3<t5)
        {
            a[i]=t3;
            p3++;
        }
        else if(t5<t2&&t5<t3)
        {
            a[i]=t5;
            p5++;
        }
        else if(t2==t3&&t2<t5)
        {
            a[i]=t2;
            p2++;
            p3++;
        }
        else if(t5==t3&&t5<t2)
        {
            a[i]=t5;
            p5++;
            p3++;
        }
        else if(t2==t5&&t2<t3)
        {
            a[i]=t2;
            p2++;
            p5++;
        }
        else
        {
            a[i]=t2;
            p2++;
            p3++;
            p5++;
        }
    }
    int n ;
   while(scanf("%d",&n)==1&&n)
    {
        printf("%d\n",a[n]);
    }
    return 0;
}