题目链接:
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~