Projecteuler 169

最新推荐文章于 2025-06-13 14:58:29 发布

JavaSxy

最新推荐文章于 2025-06-13 14:58:29 发布

阅读量244

点赞数

分类专栏：数学&编程文章标签： F# Scala PHP .net

数学&编程专栏收录该内容

11 篇文章

订阅专栏

本文探讨了Project Euler 169问题的解决方法，该问题要求计算特定数值下表达方式的数量。通过递归关系式和优化算法，最终得出10^25的情况下的解答。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

[url]http://projecteuler.net/index.php?section=problems&id=169[/url]
[quote]Define f(0)=1 and f(n) to be the number of different ways n can be expressed as a sum of integer powers of 2 using each power no more than twice.

For example, f(10)=5 since there are five different ways to express 10:

1 + 1 + 8
1 + 1 + 4 + 4
1 + 1 + 2 + 2 + 4
2 + 4 + 4
2 + 8

What is f(10^(25))?
[/quote]
注意到
f(2k+1) = f(k)
f(2k) = f(k) + f(k-1) k >= 1
容易得到，对 n = p*2^m(其中p%2==1,m>=0)，有：
f(n) = f(p) + m*f(p-1)
因此
f(n+1) = f(p*2^m+1)
= f(p*2^(m-1))
= f(p) + (m-1)*f(p-1)
从而：
x*f(n)+y*f(n+1) = (mx+(m-1)y)f(p-1)+(x+y)f(p)

因此得到下面的代码：

/**
  compute x*f(n)+y*f(n+1) where n%2 == 0
*/
def count(x: BigInt, y: BigInt, n: BigInt): BigInt = {
   if(n == 0) return x + y
   val m = n.lowestSetBit    // n = p*2^m where p%2 == 1
   val p = n >> m
   count(x*m + y*(m-1),x+y,p-1)
}

[quote]scala> count(1,0,(10: BigInt).pow(25))
res1: BigInt = 178653872807
[/quote]