C# 语言Pagerank两种实现

  最近正在学习Hadoop相关的知识，参考了很多材料，如下是C#的两种实现。

  目前对这个算法只是知道了运算的过程和思路，理解的还不是太深刻，我看到论坛上很多朋友算的pr值超过1，总感觉是有问题，具体哪有问题反倒倒说不好了，请知晓的朋友帮忙分析下。
        附件是两种实现方式，感兴趣的朋友可以自行下来看看，有不对的地方，欢迎批评指正。
        这是第二种实现的代码,我是通过判断两个矩阵的相差的范围进而判断是否退出迭代的。
       C# 语言的实现，参考http://f.dataguru.cn/forum.php?mod=viewthread&tid=58606             输出结果中的S是构建的原始矩阵

1. using System;
2. using System.Collections.Generic;
3. using System.Linq;
4. using System.Text;
6. namespace Pagerank
7. {
8.     public class  PagerankMapReduce
9.     {
10.         const double EPS = 0.000001;//精度值设置越小，迭代的次数就越多
11.         const double A = 0.85;
12.         const int NODE_COUNT = 5;
14.         public void Rank()
15.         {
16.             double[,] temp = new double[NODE_COUNT, NODE_COUNT] { 
17.             { 0,1 / 2.0, 1 / 2.0, 0, 1 / 2.0 }, 
18.             { 1 / 4.0, 0, 0, 0, 0 },
19.             { 1 / 4.0, 0, 0, 1, 1 / 2.0 }, 
20.             {1/4.0,1/2.0,1/2.0,0,0},
21.             { 1 / 4.0, 0, 0, 0, 0 } };
22.             double[] last_q_vec = new double[] { 0, 0, 0, 0, 0 };
23.             bool canBreak = false;
24.             int interCount = 1;
25.             while (!canBreak)
26.             {
27.                 double[] next_q_vec = new double[] { 0, 0, 0, 0, 0 };
28.                 for (int row = 0; row < NODE_COUNT; row++)
29.                 {
30.                     for (int col = 0; col < NODE_COUNT; col++)
31.                     {
32.                         next_q_vec[row] += temp[row, col] * last_q_vec[col];
33.                     }
34.                     next_q_vec[row] = A * next_q_vec[row] + (1 - A) / NODE_COUNT;
35.                 }
36.                 canBreak = true;
37.                 for (int i = 0; i < last_q_vec.Length; i++)
38.                 {
39.                     if (Math.Abs(last_q_vec[i] - next_q_vec[i]) > EPS)
40.                     {
41.                         canBreak = false;
42.                         break;
43.                     }
44.                 }
45.                 if (!canBreak)
46.                 {
47.                     for (int i = 0; i < next_q_vec.Length; i++)
48.                     {
49.                         last_q_vec[i] = next_q_vec[i];
50.                     }
51.                 }
52.                 interCount++;
53.             }
54.             Console.WriteLine("InterCount:{0}", interCount);
55.             PrintArray(last_q_vec);
56.         }
57.         public void PrintArray(double[] row)
58.         {
59.             int dimension = row.Length;
60.             Console.Write("[");
61.             for (int col = 0; col < dimension; col++)
62.             {
63.                 Console.Write(string.Format("{0,8}", row[col]));
64.                 if (col == dimension - 1)
65.                 {
66.                     Console.Write("]");
67.                 }
68.                 else
69.                 {
70.                     Console.Write(" ");
71.                 }
72.             }
73.             Console.WriteLine();
74.         }
75.     }
76. }

复制代码

运算结果: 

 

代码下载