matlab的雷劈数算法,寻找雷劈数

戴公子预初，看到个故事：

有位数学家卡普利加在一次旅行中，遇到猛烈的暴风雨，电闪雷鸣过后，他看到路边一块里程碑，被雷电劈成两半，一半上刻著30，另一半刻著25。这时，卡普利加的脑际中忽然发现了一个绝妙的数学关系——把劈成两半的数加起来，再平方，正好是原来的数位。除此之外，还有没有别的数，也具有这样的性质呢？

熟悉速算的人很快就找到了另一个数：2025 按照第一个发现者的名字，这种怪数被命名为“卡普利加数”，又称“雷劈数”。

据说这样的书有无穷多个，怎么找到他们呢？

----------------------------------------------------------------------------------------------->

首先想到的是excel，用row()^2，用Log找到结果是偶数位的行，对半拆开，相加平方后比较，轻松找到十来个五位数以内的对劈的雷劈数。

不过，戴公子发现在excel中公式拖到底运算量太大电脑差点死机，总共能算出几十个，再多用excel就不行了。

到此为止已经可以证明有无穷多个雷劈数了：N个9，接一个8，接N个0，再接一个1，肯定是是N+1个9的平方。

如果包含非对劈的，1后面接偶数个0的都是。

其实也不是excel不行，算法得改进，不过要更好解决，还是交给python吧。

-------------------------------------------