UVA - 10856 Recover Factorial（二分查找）

题目：（Sample Output是错的！要不是因为题目里面有描述输出，估计没人做的对吧

）



这个题目，和这个OJ里面的另外一个题目有着紧密的联系。点击打开我的博客

（这一点，其实从标题也看的出来）

把这个题目的代码改改就能得到这个题目的代码。

在编代码的过程中就已经用到了二分了，可以确定，需要的数组大小的大致范围，大于2700000，小与2750000

代码：

#include<iostream>
using namespace std;

int sum[2750000];

bool isprime(int n)
{
if (n == 2)return true;
if (n % 2 == 0)return false;
for (int i = 3; i*i <= n; i += 2)if (n%i == 0)return false;
return true;
}

int main()
{
for (int i = 0; i < 2750000; i++)sum[i] = 0;
for (int i = 2; i < 2750000; i++)
{
if (isprime(i))
{
long long m = i;
while (m < 2750000)
{
for (int j = m; j < 2750000; j += m)sum[j]++;
m *= i;
}
}
}
for (int i = 1; i < 2750000; i++)sum[i] += sum[i - 1];
int n, low, high, mid;
int cas = 1;
while (cin >> n)
{
if (n < 0)break;
cout << "Case " << cas++ << ": ";
low = 0, high = 2749999;
while (low < high - 1)
{
mid = (low + high) / 2;
if (sum[mid] < n)low = mid;
else if (sum[mid]>n)high = mid;
else break;
}
if (n == 0)mid = 0;
if (sum[mid] == n)cout << mid << "!\n";
else cout << "Not possible.\n";
}
return 0;
}

在查找的时候，还是二分，因为数组已经是升序的了。

不过，有个特殊情况，sum[0]=sum[1]=0，这是唯一的不是严格递增的地方，需要特判。