题目链接:http://tz.ahalei.com/problems/view/48
思路:我的方法是暴力穷举(也就是没有方法,虽然通过简单的推理可以知道个位数只能是0,1,5,6,但没有大用处)。但这题要注意,穷举的i最大达到1e8,平方i*i达到1e16,超过了32位整型的表达范围,故需要用64位整型或double。
#include <iostream>
#include <ctime>
#include <cstdlib>
using namespace std;
clock_t start__;
#define tic start__=clock()
#define toc cout<<(clock()-start__)*1000/CLOCKS_PER_SEC<<"ms\n"
int fun(int k)
{// k是当前要统计的位数范围
unsigned __int64 i;
int sum=0, l, r;
for(l=i=1;i<k;i++) l*=10;
r=l*10;
for(i=l;i<r;i++) if( (i*i)%r==i )
sum++,printf("%I64u\n",i);
return sum;
}
int main()
{
int k, sum=0;
cout << "守形数列表:" << endl;
tic;
for(k=1;k<10;k++) sum += fun(k);
cout << "用时"; toc;
cout << "1亿以内共有" << sum << "个守形数" << endl;
return 0;
}
网站的答案15是错误的!我用Mathematica软件验算过上述17个数,都是准确无误的。另一方面,无符号64位整型的计算范围达到18e18,完全超过题目计算最大中间值1e16,所以C++跑出的答案肯定也是精确的。