ネットワークを構成するノードとして機能するコンピュータであって、自ノードに隣接する隣接ノードの仮圧力値および近似圧力値を前記隣接ノードから取得する隣接ノード通信部と、前記自ノードに接続する各リンクの管径値および長さと、前記隣接ノードの前記仮圧力値と、予め与えられた前記自ノードの湧き出しまたは吸い込みの流量とを用いて、前記自ノードの近似圧力値を計算する近似圧力値計算部と、前記管径値と、前記長さと、前記隣接ノードの前記近似圧力値と、前記自ノードの前記近似圧力値とを用いて、前記各リンクを流れる流量を計算する流量値計算部と、前記各リンクを流れる前記流量と、前記管径値と、粘菌の特徴を表す関数とを用いて、前記管径値を更新する管径値更新部とを備え、少なくとも前記管径値が収束するまで、前記近似圧力値の計算、前記流量の計算、および、前記管径値の更新を繰り返すコンピュータ。前記近似圧力値計算部は、前記自ノードの仮圧力値の初期値を設定する仮圧力初期値設定部と、前記自ノードの前記仮圧力値に所定の反復式を適用し、得られた仮圧力値を前記自ノードの前記近似圧力値として出力する仮圧力値反復式適用部とを有する請求項1に記載のコンピュータ。前記仮圧力値反復式適用部は、前記自ノードをiとし、前記管径値をDijとし、前記長さをLijとし、前記隣接ノードの前記仮圧力値をπjk-1とし、前記自ノードの湧き出しまたは吸い込みの前記流量をIiとし、k-1回目の反復における前記自ノードの前記仮圧力値をπik-1とし、前記隣接ノードの集合をnbriとし、予め与えられた0<sp<1を満たす値をspとしたときに、次式に従って、k回目の反復における前記自ノードの前記仮圧力値πikを求める請求項2に記載のコンピュータ。前記仮圧力初期値設定部は、前記仮圧力値の前記初期値が設定される時点の前記自ノードの前記近似圧力値を前記仮圧力値の前記初期値として設定する請求項2又は3に記載のコンピュータ。複数のコンピュータを備え、前記複数のコンピュータのそれぞれは請求項1~請求項4のいずれか1項に記載のコンピュータであって、前記自ノードの湧き出しまたは吸い込みの前記流量として、前記複数のコンピュータのそれぞれが、湧き出しノードか、吸い込みノードか、または、前記湧き出しノード及び前記吸い込みノード以外のノードかに応じて決められた値を前記複数のコンピュータのそれぞれに与え、前記複数のコンピュータのそれぞれに保持されている前記管径値に基づいて、前記管径値が所定の値に収束したリンクを抽出し、抽出された前記リンクに基づいて、前記湧き出しノードと前記吸い込みノードとの間の経路を決定する経路探索システム。ネットワークを構成するノードとして機能するコンピュータを制御する制御方法であって、自ノードに隣接する隣接ノードの仮圧力値および近似圧力値を前記隣接ノードから取得し、前記自ノードに接続する各リンクの管径値および長さと、前記隣接ノードの前記仮圧力値と、予め与えられた前記自ノードの湧き出しまたは吸い込みの流量とを用いて、前記自ノードの近似圧力値を計算し、前記管径値と、前記長さと、前記隣接ノードの前記近似圧力値と、前記自ノードの前記近似圧力値とを用いて、前記各リンクを流れる流量を計算し、前記各リンクを流れる前記流量と、前記管径値と、粘菌の特徴を表す関数とを用いて、前記管径値を更新し、少なくとも前記管径値が収束するまで、前記近似圧力値の計算、前記流量の計算、および、前記管径値の更新を繰り返す制御方法。ネットワークを構成するノードとして機能するコンピュータ用のプログラムであって、前記コンピュータを、自ノードに隣接する隣接ノードの仮圧力値および近似圧力値を前記隣接ノードから取得する隣接ノード通信部と、前記自ノードに接続する各リンクの管径値および長さと、前記隣接ノードの前記仮圧力値と、予め与えられた前記自ノードの湧き出しまたは吸い込みの流量とを用いて、前記自ノードの近似圧力値を計算する近似圧力値計算部と、前記管径値と、前記長さと、前記隣接ノードの前記近似圧力値と、前記自ノードの前記近似圧力値とを用いて、前記各リンクを流れる流量を計算する流量値計算部と、前記各リンクを流れる前記流量と、前記管径値と、粘菌の特徴を表す関数とを用いて、前記管径値を更新する管径値更新部として機能させ、少なくとも前記管径値が収束するまで、前記近似圧力値の計算、前記流量の計算、および、前記管径値の更新を繰り返させるプログラム。ところが、非特許文献1及び非特許文献2の方式では、管径発展のシミュレーションを各ノードで分散させて行うことができない。これは、式2を厳密に解いて各ノードの圧力を求めようとすると、ネットワーク内における全てのリンクの長さと管径を知る必要があるためである。本発明の目的は、上述した課題を解決する経路探索システム、コンピュータ、制御方法、及びプログラムを提供することにある。上記課題を解決するために、本発明は、ネットワークを構成するノードとして機能するコンピュータであって、自ノードに隣接する隣接ノードの仮圧力値および近似圧力値を前記隣接ノードから取得する隣接ノード通信部と、前記自ノードに接続する各リンクの管径値および長さと、前記隣接ノードの前記仮圧力値と、予め与えられた前記自ノードの湧き出しまたは吸い込みの流量とを用いて、前記自ノードの近似圧力値を計算する近似圧力値計算部と、前記管径値と、前記長さと、前記隣接ノードの前記近似圧力値と、前記自ノードの前記近似圧力値とを用いて、前記各リンクを流れる流量を計算する流量値計算部と、前記各リンクを流れる前記流量と、前記管径値と、粘菌の特徴を表す関数とを用いて、前記管径値を更新する管径値更新部とを備え、少なくとも前記管径値が収束するまで、前記近似圧力値の計算、前記流量の計算、および、前記管径値の更新を繰り返す。また、本発明は、ネットワークを構成するノードとして機能するコンピュータを制御する制御方法であって、自ノードに隣接する隣接ノードの仮圧力値および近似圧力値を前記隣接ノードから取得し、前記自ノードに接続する各リンクの管径値および長さと、前記隣接ノードの前記仮圧力値と、予め与えられた前記自ノードの湧き出しまたは吸い込みの流量とを用いて、前記自ノードの近似圧力値を計算し、前記管径値と、前記長さと、前記隣接ノードの前記近似圧力値と、前記自ノードの前記近似圧力値とを用いて、前記各リンクを流れる流量を計算し、前記各リンクを流れる前記流量と、前記管径値と、粘菌の特徴を表す関数とを用いて、前記管径値を更新し、少なくとも前記管径値が収束するまで、前記近似圧力値の計算、前記流量の計算、および、前記管径値の更新を繰り返す。また、本発明は、ネットワークを構成するノードとして機能するコンピュータ用のプログラムであって、前記コンピュータを、自ノードに隣接する隣接ノードの仮圧力値および近似圧力値を前記隣接ノードから取得する隣接ノード通信部と、前記自ノードに接続する各リンクの管径値および長さと、前記隣接ノードの前記仮圧力値と、予め与えられた前記自ノードの湧き出しまたは吸い込みの流量とを用いて、前記自ノードの近似圧力値を計算する近似圧力値計算部と、前記管径値と、前記長さと、前記隣接ノードの前記近似圧力値と、前記自ノードの前記近似圧力値とを用いて、前記各リンクを流れる流量を計算する流量値計算部と、前記各リンクを流れる前記流量と、前記管径値と、粘菌の特徴を表す関数とを用いて、前記管径値を更新する管径値更新部として機能させ、少なくとも前記管径値が収束するまで、前記近似圧力値の計算、前記流量の計算、および、前記管径値の更新を繰り返させる。なおまた、上記の発明の概要は、本発明の必要な特徴の全てを列挙したものではなく、これらの特徴群のサブコンビネーションもまた、発明となり得る。以上の説明から明らかなように、この発明は、ネットワーク上の経路探索問題を、各ノードがネットワーク全体の情報を知ることなく、自律分散的に解くシステムを提供できる。以下、実施形態および実施例を通じて本発明を説明する。ただし、以下の実施形態および実施例は、請求の範囲にかかる発明を限定するものではない。また、実施形態および実施例の中で説明されている特徴の組み合わせの全てが本発明の解決手段として必須であるとは限らない。本発明の実施形態は、ネットワーク内の各ノードが、図1に示す機能を持つことを特徴とする。図1を参照すると、各ネットワークノード100は、自ノードの近似圧力値を求める近似圧力値計算部110と、自ノードに接続するリンク隣接リンクを流れる流量を求める流量値計算部120と、自ノードに接続するリンクの管径を更新する管径値更新部130と、自ノードに隣接するノード隣接ノードと通信する隣接ノード通信部140とから構成されている。近似圧力計算部110は、仮圧力値の初期値を設定する仮圧力値初期値設定部111と、仮圧力値を反復更新する仮圧力値反復式適用部112と、を含む。これらの機能ブロックはそれぞれ概略次のように動作する。なお、以下の説明ではノードiにおける動作について述べるが、どのノードについても動作は同様である。仮圧力初期値設定部111は、自ノードの近似圧力値を求めるための、仮圧力値の初期値を設定する。初期値の設定値は任意である。例えば、初期値を設定する時点での近似圧力値を用いた場合、仮圧力値の初期値は式4で表される。ここで、piはノードiの近似圧力値、πikはkkは0以上の整数回の反復式を適用したときのノードiの仮圧力値を示す。仮圧力値反復式適用部112は、式5で表される仮圧力値に関する反復式を、予め与えられた所定の回数Nだけ適用し、得られた仮圧力値πiNを近似圧力値piとして出力する式6を参照。ただし、spは予め与えられた、0<sp<1を満たす所定の値である。式5では、式中の2箇所において、ネットワーク内のノードiを除く全てのノードに関して和が取られている。ただ、ノードiとの間にリンクが存在しないノードjに関しては、Dij=0であるため、式5は式7と等価である。ここで、nbriはノードiと隣接するノードのインデックス集合を表す。流量値計算部120は、式1を用いてノードiと隣接するノードとの間の各リンクの流量Qijを求める。管径値更新部130は、式8を用いてノードiと隣接するノードとの間の各リンクの管径Dijを更新する。なお、sdは予め与えられた0<sd<1を満たす所定の値である。隣接ノード通信部140は、ノードiに隣接する各ノードと通信するブロックである。隣接ノード通信部140は、仮圧力値反復式適用部112又は流量値計算部120からの要求に応じて、仮圧力値又は近似圧力値の送受信を行う。なお、ネットワークノード100とは別に設けられた装置が、探索された経路を構築する場合、隣接ノード通信部140は、管径値更新部130から管径値を取得し、取得した管径値をネットワークノード100の外部へ送信する。通信のプロトコルは任意のものを用いてよい。次に、図1及び図2のフローチャートを参照して本実施形態の全体の動作について詳細に説明する。まず、初期化として、自ノードに接続するリンクのリンク長Lij及び管径の初期値Dij、自ノードからの湧き出し又は吸い込み量を表すIi及び近似圧力値の初期値piが与えられるステップS201。なお、Iiに関しては、湧き出しノードに相当するコンピュータに対して-1が与えられ、吸い込みノードに相当するコンピュータに対して1が与えられ、これら以外のノードに相当するコンピュータに対して0が与えられる。次に、仮圧力初期値設定部111が所定の方法例えば、式4に従った方法でノードiの仮圧力値πi0を設定するステップS202。次に、仮圧力値反復式適用部112が反復回数kの値を1に初期化するステップS203。仮圧力値反復式適用部112は、隣接ノード通信部140に各隣接ノードの仮圧力値πjk-1を取得させるステップS204。更に、仮圧力値反復式適用部112は、自ノードと隣接ノードとの間の各リンク隣接リンクの管径Dijおよびリンク長Lijと、各隣接ノードの仮圧力値πjk-1と、自ノードの仮圧力値πik-1と、自ノードから湧き出し又は自ノードへ吸い込まれる流量Iiと、値spとを用いて、式7に従って、自ノードの仮圧力値πikを計算するステップS205。仮圧力値反復式適用部112は、ステップS204とステップS205を予め与えられた所定の回数Nに達するまで繰り返すステップS206、ステップS207。N回の繰り返しが終了後、仮圧力値反復式適用部112は、得られた仮圧力値πiNを自ノードの近似圧力値に代入するステップS208。次に、流量値計算部120が隣接ノード通信部140に各隣接ノードの近似圧力値pjを取得させるステップS209。更に、流量値計算部120は、各隣接リンクの管径Dijおよびリンク長Lijと、各隣接ノードの近似圧力値pjと、自ノードの近似圧力値piとを用いて、式1に従って、各隣接リンクに対応する流量Qijの値を更新するステップS210。次に、管径値更新部130が、各隣接リンクの現在の管径Dijと、各隣接リンクの流量Qijと、値sdと、関数fとを用いて、式8に従って、各隣接リンクに対応する管径値Dijを更新するステップS211。以上のステップS202からステップS211までの処理が所定の終了条件が満たされるまで繰り返されるステップS212。終了条件としては、例えば、予め与えられた所定の回数だけステップS202からステップS211までの処理が繰り返されたとき等の条件を用いる。ただし、少なくとも各隣接リンクの管径が0または1に収束するまでは、ステップS202からステップS211までの処理を繰り返す。上述したように、管径が1に収束したリンクのみを抽出することで経路が探索されたことになるため、抽出されたリンクから探索された経路を構築する。例えば、原形質の湧き出しノードと吸い込みノードとの間の経路上の各ネットワークノード100が、管径が1に収束した2本の隣接リンクを抽出し、このうち一方のリンクから受け取った通信を、他方のリンクへと転送する。その結果、原形質の湧き出しノードと吸い込みノード間の通信経路が構築される。あるいは、図5に示したように、ネットワークノード100とは別にネットワーク内にホストコンピュータ200を設ける。なお、図5では、煩雑さを避けるため、一部のネットワークノードについてのみ参照符号100を付している。ホストコンピュータ200は、各ネットワークノード100から各隣接リンクに対応する管径値を収集し、収集した管径値をもとに管径が1のリンクを抽出し、抽出されたリンクを接続して湧き出しノードと吸い込みノードとの間の経路を構築する。あるいは、各ネットワークノード100が各隣接リンクに対応する管径値をもとに管径が1のリンクを抽出し、ホストコンピュータ200が当該リンクに関する情報を収集し、収集された情報で示されるリンクを接続して経路を構築するようにしても良い。あるいは、ネットワークノードとは別にホストコンピュータ200を設けるのではなく、特定のネットワークノード例えば、予め決められたノード、原形質の湧き出しノード、あるいは、原形質の吸い込みノードが、ホストコンピュータ200と同様の動作を行うようにしても良い。なお、ホストコンピュータ200または特定のネットワークノードが各ネットワークノード100または他のネットワークノードから管径値または管径が1のリンクに関する情報を取得するためには、ネットワークの分野で知られている種々の手法を用いることができる。次に、本実施形態の効果について説明する。本実施形態では、自ノードが保持している値と、隣接ノードの仮圧力値からのみ構成される式7を用いて近似圧力値を求めている。また、近似圧力値を求めるためのステップ以外の他のステップも全て自ノードが保持している値及び隣接ノードから取得できる値のみを用いて実行可能である。その結果、本実施形態では、真性粘菌のシミュレーションを各ノードで自律分散的に実現させている。なお本発明は、経路探索システムを構築するノードとしてコンピュータを機能させる経路探索用プログラムとしても実現することができる。コンピュータは、経路探索用プログラムを読み込み、読み込まれた経路探索用プログラムを実行する中央処理装置CentralProcessingUnit;CPU、自ノードの仮圧力値および近似圧力値、各隣接リンクのリンク長と管径などを保存する記憶装置ハードディスク等、及び隣接するコンピュータと通信するための通信インタフェースを備える。次に、本発明の具体的な実施例を示す。ここでは、sp=sd=0.5、N=1、全てのノードiについて近似圧力値の初期値piを0、全てのリンクijについて管径の初期値Dij=0とし、式8の関数fはfx=absxとした。対象のネットワークとしてはノード数100のBarabashi-Albertモデルを用いており、ランダムに選んだ2点間の最短経路を求めた。図3に、最短経路shortest1またはshortest2上および非最短経路non-shortest1またはnon-shortest2上の時間発展の様子を示す。また、図4に得られた経路を示す。最短経路上の管径が1に収束し、非最短経路上の管径が0に収束しており、正しく最短経路が得られている。以上、本発明を実施形態および実施例を用いて説明したが、本発明の技術的範囲は、上記実施形態および実施例に記載の範囲には限定されない。上記実施形態および実施例に、多様な変更又は改良を加えることが可能であることは当業者に明らかである。そのような変更又は改良を加えた形態も本発明の技術的範囲に含まれ得ることが、請求の範囲の記載から明らかである。この出願は、2010年9月16日に出願された日本出願特願2010-208083号を基礎とする優先権を主張し、その開示の全てをここに取り込む。