posts - 5,comments - 8,views - 6906

想到了一个新的解法,求20亿以内的双立方数只用1秒。

最近在园子里看到了几篇关于这个双立方数的文章,自己也很有兴趣,希望能找到一个好的算法。

题目要求:双倍超立方数是指一个正整数可以正好被拆分为两种不同的a^3+b^3的方式,其中a,b均为整数且0<a<=b。对于任何一个指定的 int n, 返回所有的小于等于n的双倍超立方数的个数。

最开始想到的算法是这样的:

 1。找到所有的小于给定的最大整数N的所有立方数,就是满足 y=x*x*x (x,y都为正整数)的 y,并把它们都存在一个Hash表里。

 2。 遍历所有从1到N的每一个整数n,作如下处理:

    2.1 遍历Hash表中的每一个小于等于 n/2 的立方数 y,如果 n-y 存在于 Hash表中,则说明 n可以表述为两个立方数的和。

    2.2 如果正好存在两个 y ,满足2.1的条件,则n满足题设条件。

 3。输出满足条件的n的个数。

 但是这样的思路,无论怎么优化,一旦N>10^8,都不会小于一分钟。

后来,突发奇想,我为什么要测试每一个n是不是满足条件呢,为什么不倒过来,通过y来计算n呢?

于是有了现在的这个算法:

1。 找到所有的小于给定的最大整数N的所有立方数,就是满足 y=x*x*x (x,y都为正整数)的 y,并把它们都存在一个链表(或者是数组)List里。

2。对List的每一个立方数y1,作如下处理:

  对List里的每一个大于等于y1的立方数y2,计算y1+y2,并把计算结果保存在一个字典 Dic 里,key为y1+y2,value为 1 ,如果字典中已存在key=y1+y2,

  则 Dic[y1+y2]的value加1。

 3。找到Dic中所有value为2的key,并输出其个数。

由于List的中的元素数量(小于N的开立方)远远小于N,所以算法效率有了质的提高,在N=2*10^9的情况下,小于1秒。

下面是实现:

 1        private static int GetDoubleCubeCount(int n)
 2        {
 3            int[] arCubes = new int[(int)Math.Floor(Math.Pow(n, 1.0 / 3))];
 4            Dictionary<intint> dicDblTubes = new Dictionary<intint>();
 5
 6            int t = 1;
 7            int index = 0;
 8            int doubleCube = 0;
 9
10            while (index < arCubes.Length)
11            {
12                arCubes[index++= t * t * t;
13                t++;
14            }

15
16            for (int i = 0; i < arCubes.Length; i++)
17            {
18                for (int j = i; j < arCubes.Length; j++)
19                {
20                    doubleCube = arCubes[i] + arCubes[j];
21                    if (doubleCube > n) break;
22
23                    if (dicDblTubes.ContainsKey(doubleCube))
24                    {
25                        dicDblTubes[doubleCube]++;
26                    }

27                    else
28                    {
29                        dicDblTubes.Add(doubleCube, 1);
30                    }

31                }

32            }

33
34            n = 0;
35
36            foreach (int value in dicDblTubes.Values)
37            {
38                if (value == 2)
39                {
40                    n++;
41                }

42            }

43            return n;
44        }

 

上面的实现仅仅是对思路的完整实现,应该还有很大的优化空间,哪位朋友要是有更好的想法,或者是优化的方法,还希望能多多交流。 由于俺的电脑是单核的,没有优化算法适应多核,多核的话,应该更快啦。

 

posted @ 2009-06-05 17:30 yujiasw 阅读(2637) 评论(6) 推荐(0) 编辑
摘要: 在上一篇文章 NET 并行扩展(ParallelFX) 试用(上) 中,我对任务并行库做了一个简单的介绍,对Parallel.For() 方法作了几个实验,对其线程的分配和异常的处理有了一定的了解。虽然只是针对这一个方法,但是它其实在某种程度上反映了任务并行库一些内部的处理方法,有一定的普遍性。 上一篇文章的示例中,我们的例子是并行循环,它特点是共享索引,各个任务之间没有任何的关... 阅读全文
posted @ 2007-12-21 17:38 yujiasw 阅读(2181) 评论(2) 推荐(0) 编辑
摘要: 单核处理器单靠提高频率来提升性能变得越来越困难,多核处理器已经是大势所趋。目前绝大部分新买的电脑都已经是双核的了,今后还会有更多的核心。 与硬件的多核形成对比的是,软件开发领域,还没有成熟的技术让我们方便的使用多核所带来的性能提升,大多数的执行过程仍然是单线程的,因为写多线程的应用很困难,除非涉及到大量的并发IO,否则一般不会使用多线程(至少我是这样,因为我..懒)。 值得... 阅读全文
posted @ 2007-12-20 13:09 yujiasw 阅读(1015) 评论(0) 推荐(0) 编辑
摘要: 通常我们见到的单件模式的实现只限于单个执行过程.然而如果要求某个公共对象对多个程序实例,甚至是不同程序都只能有一个实例,该怎么办呢? 首先,我们想到这个对象一定是跨进程,跨应用程序域的.应用程序域是系统内程序执行的安全和资源边界,不同的应用程序域其内存都是独立的,使用普通的方式,是不可能实现的. 要跨越应用程序域,我们似乎只有一个选择,使用Remoting技术,使不同的程... 阅读全文
posted @ 2007-12-17 14:10 yujiasw 阅读(647) 评论(0) 推荐(0) 编辑
摘要: 在CSDN 上 看了一篇Daniel Cazzulinode的文章:Linq的超越 (原文链接http://blog.csdn.net/programmer_editor/archive/2006/09/29/1305859.aspx) 挺有意思的,自己也研究了一下. 题目可能没有说清楚,那么咱么现在考虑一个场景: 为了能够显示不同数据对象的数据,提高可复用性,你写了一个... 阅读全文
posted @ 2007-12-01 18:25 yujiasw 阅读(426) 评论(0) 推荐(0) 编辑
< 2025年3月 >
23 24 25 26 27 28 1
2 3 4 5 6 7 8
9 10 11 12 13 14 15
16 17 18 19 20 21 22
23 24 25 26 27 28 29
30 31 1 2 3 4 5

点击右上角即可分享
微信分享提示