projecteuler No.100 Arranged probability

题目链接：

Problem 100 Arranged probability

题目100：要使取到两个蓝碟子的几率是50%需要有多少个蓝碟子？

通过人数：7708

题目分析：

这道题记盘子数为n，蓝盘子为p，则相当于要解方程：

n(n-1)=2p(p-1)的第一组使得n>10^12的正整数解的p值。

对于这个方程我最开始还是尝试蛮力解，不过几分钟的时间只从10^12之后搜了1亿个数...于是显然蛮力解法不行...

于是恒等变形成，等价于如下方程：

(2n-1)^2-2(2p-1^2)=-1

这是第二型佩尔方程。

而关于第二型佩尔方程的正整数解的通用解法，好像还是一个未解决问题。

毕竟在谷歌学术中搜索第II型佩尔方程，会查到不少关于这个问题的子问题或者交集问题的论文。

于是就考虑曲线的求解方法。

一方面，通过一篇博文（Project Euler 100: Finding the number of blue discs for which there is 50% chance of taking two blue.）,得知了这个网站：http://www.alpertron.com.ar/QUAD.HTM，给出了二元二次不定方程正整数解的在线求解工具。

另一方面，考虑到(1,1)为这个方程的最小解，记为(x1,y1)则第二型佩尔方程有通解公式：

xn=(1/2)((x1+y1sqrt(d))^(2n+1)+(x1-y1sqrt(d))^(2n+1))

yn=(1/(2sqrt(d)))((x1+y1sqrt(d))^(2n+1)-(x1-y1sqrt(d))^(2n+1))

其中，d=2.

于是，构成线性地推：

x_(n+1) = 3x_n+4y_n

y_(n+1) = 2x_n+3y_n

这样就可以求解了。

解题过程：

利用上面的分析，可以写出如下代码：

long long int x[100];
long long int y[100];

x[0] = 1;
y[0] = 1;

for (int i = 1;i <= 16; i++)
{
	x[i] = 3 * x[i-1] + 4 * y[i-1];
	y[i] = 2 * x[i-1] + 3 * y[i-1];

	if (x[i] > (long long int)1000000000000)
	{
		printf("%lld\n",(y[i] + 1) / 2);
		break;
	}
}
return 0;

最后打印出结果： 756872327473

以上只是我做题时的解法。

如果有更好的解法、更好的思路，欢迎评论讨论~O(∩_∩)O~