金蝉平方数

       参考文献：郭继展,郭勇,苏辉《程序算法与技巧精选》,机械工业出版社,2008年5月第一版,ISBN 978-7-111-23816-4,第5.6节。

       定义：设s是一个平方数，由1,2,3,4,5,6,7,8,9,0中无重复的数字组成。s经一次脱壳（即同时去掉其最高位与最低位数字）后剩下的数还是平方数。如果剩下的数仍多于2位，再行脱壳，剩下的数仍是平方数。以此类推，直至最后剩下的数（一位或两位）仍是平方数，这样的数s称为金蝉平方数。

       例如：38416是196的平方，38416脱壳后的841是29的平方，841脱壳后的4是2的平方，所以38416是金蝉平方数。

       问题：求解出所有的金蝉平方数。

       算法分析

       （1）s的值必定小于9876543210，否则数字必定重复。

       （2）本章思想，平方数问题——算法多从数的平方下手，枚举0~9876543210内的平方数（3万多个），而不是枚举0~9876543210内的每一个数，这是算法优化的关键。

       （3）得到平方数后，先判断各个数字的出现有无重复。（用isudig函数实现）

       （4）再用一个递归函数判断当前或多次脱壳后是否为平方数。（用isjcpfs函数实现，jcpfs即"金蝉平方数"拼音首字母组合）

      实现代码

       “mymath.h”：http://blog.csdn.net/code4101/article/details/17993747

       “debug.h”：http://blog.csdn.net/code4101/article/details/17993293

#include <iostream>
#include "mymath.h" 
#include "debug.h"
using namespace std;

int isjcpfs(__int64 v)
{
	if(fmod(sqrt(v),1)) return 0;		// 当且仅当fmod(sqrt(v),1)值为0时,v为平方数
	else if(v<100) return 1;
	v/=10;					// 脱去最低位
	v%=(int)(pow(10, (int)log10(v)));	// 脱去最高位
	return isjcpfs(v);
}

int main()
{
	int		c=0;
	__int64 i;

	tic;
	cout << "所有的金蝉平方数:\n";
	for(i=0; i<sqrt(9876543210); i++) if(isudig(i*i)) if(isjcpfs(i*i))
		printf("%2d: %6I64d\n",++c,i*i);
	cout<<"计算用时"; toc;

	return 0;
}

        运行结果截图