好久不来CSDN,偶尔看到这个帖子:http://community.csdn.net/Expert/topic/4730/4730797.xml?temp=5.944461E-02,楼主说是微软的面试题。我第一次答错了,然后才领会到要点:这题适宜用动态规划算法来解。
下面的代码可能有点笨拙,但自以为清晰,暂且存起来。
#include <cstdlib>
#include <cstdio>
class Calc235
{
unsigned *v;
unsigned minV, va, vb, vc;
int a,b,c;
const int N;
int min(unsigned a, unsigned b, unsigned c)
{
return (a<=b)?((a<=c)?a:c):((b<=c)?b:c);
}
public:
Calc235(int n)
:N(n)
{
v = new unsigned[N];
}
~Calc235()
{
delete[]v;
}
void execute()
{
v[0] = 1;
a = b = c = 0;
va = v[a] * 2;
vb = v[b] * 3;
vc = v[c] * 5;
for (int i=1; i<N; ++i)
{
v[i] = min(va, vb, vc);
if (va == v[i])
{
va = v[++a] * 2;
}
if (vb == v[i])
{
vb = v[++b] * 3;
}
if (vc == v[i])
{
vc = v[++c] * 5;
}
// printf("v[%i]=%10u v[%i]=%10u v[%i]=%10u v[%i]=%10u/n",
// i, v[i], a, v[a]*2, b, v[b]*3, c, v[c]*5);
}
}
unsigned getResult()
{
return v[N-1];
}
};
int main(int argc, char *argv[])
{
int N = 1500;
if (argc > 1)
{
N = atoi(argv[1]);
}
Calc235 calc(N);
calc.execute();
unsigned v = calc.getResult();
printf("result = %u/n", v);
return 0;
}
计算结果:
命令“calc235 2000”的部分输出,这一部分显示,第1500个只有2、3、5质因子的自然数是859963392。
v[1498]= 854296875 v[1363]= 859963392 v[1287]= 859963392 v[1197]= 860934420
v[1499]= 859963392 v[1364]= 860934420 v[1288]= 860934420 v[1197]= 860934420
v[1500]= 860934420 v[1365]= 864000000 v[1289]= 864000000 v[1198]= 864000000
v[1501]= 864000000 v[1366]= 874800000 v[1290]= 874800000 v[1199]= 874800000
下面的部分输出显示,计算第1847个数时发生了溢出:
v[1829]=3955078125 v[1674]=3981312000 v[1589]=3981312000 v[1483]=3981312000
v[1830]=3981312000 v[1675]=3985807500 v[1590]=3985807500 v[1484]=3985807500
v[1831]=3985807500 v[1676]=4000000000 v[1591]=4026531840 v[1485]=4000000000
v[1832]=4000000000 v[1677]=4026531840 v[1591]=4026531840 v[1486]=4026531840
v[1833]=4026531840 v[1678]=4031078400 v[1592]=4031078400 v[1487]=4031078400
v[1834]=4031078400 v[1679]=4050000000 v[1593]=4050000000 v[1488]=4050000000
v[1835]=4050000000 v[1680]=4076863488 v[1594]=4076863488 v[1489]=4081466880
v[1836]=4076863488 v[1681]=4081466880 v[1595]=4081466880 v[1489]=4081466880
v[1837]=4081466880 v[1682]=4096000000 v[1596]=4100625000 v[1490]=4096000000
v[1838]=4096000000 v[1683]=4100625000 v[1596]=4100625000 v[1491]=4100625000
v[1839]=4100625000 v[1684]=4132485216 v[1597]=4132485216 v[1492]=4147200000
v[1840]=4132485216 v[1685]=4147200000 v[1598]=4147200000 v[1492]=4147200000
v[1841]=4147200000 v[1686]=4194304000 v[1599]=4199040000 v[1493]=4194304000
v[1842]=4194304000 v[1687]=4199040000 v[1599]=4199040000 v[1494]=4199040000
v[1843]=4199040000 v[1688]=4218750000 v[1600]=4218750000 v[1495]=4218750000
v[1844]=4218750000 v[1689]=4246732800 v[1601]=4246732800 v[1496]=4246732800
v[1845]=4246732800 v[1690]=4251528000 v[1602]=4251528000 v[1497]=4251528000
v[1846]=4251528000 v[1691]= 0 v[1603]=4271484375 v[1498]=4271484375
v[1847]= 0 v[1692]= 4849664 v[1603]=4271484375 v[1498]=4271484375
v[1848]= 4849664 v[1693]= 9704804 v[1603]=4271484375 v[1498]=4271484375
v[1849]= 9704804 v[1694]= 25032704 v[1603]=4271484375 v[1498]=4271484375
v[1850]= 25032704 v[1695]= 58597376 v[1603]=4271484375 v[1498]=4271484375
v[1851]= 58597376 v[1696]= 79032704 v[1603]=4271484375 v[1498]=4271484375
v[1852]= 79032704 v[1697]= 99563954 v[1603]=4271484375 v[1498]=4271484375
v[1853]= 99563954 v[1698]= 128712704 v[1603]=4271484375 v[1498]=4271484375
v[1854]= 128712704 v[1699]= 133707704 v[1603]=4271484375 v[1498]=4271484375
v[1855]= 133707704 v[1700]= 184008704 v[1603]=4271484375 v[1498]=4271484375