"周期的に繰り返される繰り返し期間中の曜日と週との組合せ毎に、放射線治療の各プロトコールに対して基点となる治療工程の日を示すように当該各プロトコールを割当てたスケジュールパターンを生成する繰り返し型スケジュール生成装置であって、前記各治療工程における治療部位と照射治療工程の回数とを示す前記各プロトコールと、前記各プロトコールの患者数とを規定した割当作業情報を記憶する第1の記憶手段と、前記各プロトコールに対する各治療工程及び当該各治療工程のうちの基点となる治療工程を規定した治療工程情報を記憶する第2の記憶手段と、前記各治療工程間に空ける日数及び前記照射治療工程の回数に関する制約を規定した制約情報とを記憶する第3の記憶手段と、前記各プロトコールを割当て可能な曜日と週との組合せを個別に示す割当時間情報を記憶する第4の記憶手段と、前記割当時間情報内の曜日と週との組合せに個別に割り当てられた曜日種類を規定した割当時間種類情報を記憶する第5の記憶手段と、前記割当時間情報及び前記割当作業情報に基づいて、当該割当時間情報に示される曜日と週との組合せを個別に示す複数の割当時間ノードと、当該割当作業情報に規定された各プロトコールを個別に示す複数のプロトコールノードと、前記各割当時間ノードと前記各プロトコールノードとの組合せを個別に示す複数の割当アークと、前記各プロトコールの患者数とを含む割当ネットワークを作成する割当ネットワーク作成手段と、前記割当作業情報及び前記割当時間情報に基づいて、前記各割当アーク毎に、割当て可能な患者数としての流量を制限する流量制限値を算出する流量制限値算出手段と、前記割当時間情報及び前記割当時間種類情報に基づいて、前記曜日種類を個別に示す複数の割当時間種類ノードと、前記各割当時間種類ノードのうちの有効な割当時間種類ノードを個別に示す複数の割当時間種類アークとを生成する種類要素生成手段と、前記割当作業情報、前記割当時間情報及び前記割当時間種類情報に基づいて、前記各プロトコール毎に、前記各割当時間種類アークの使用枝数を制限する使用枝数制限値を算出する使用枝数制限値算出手段と、前記割当ネットワーク、前記流量制限値、前記各割当時間種類ノード、前記各割当時間種類アーク及び前記使用枝数制限値を統合することにより、前記各割当時間ノードと、前記各プロトコールノードと、前記各割当アークと、前記患者数と、前記流量制限値と、前記各割当時間種類ノードと、前記各割当時間種類アークと、前記使用枝数制限値とを含む制限ネットワークを作成する制限ネットワーク作成手段と、前記制限ネットワーク、前記治療工程情報及び前記制約情報に基づいて、前記スケジュールパターンを生成し、当該生成したスケジュールパターンをスケジュールパターン記憶手段に書込むスケジュールパターン生成手段と、前記書込まれたスケジュールパターンを記憶する前記スケジュールパターン記憶手段と、前記スケジュールパターン記憶手段からスケジュールパターンを読出し、当該読出したスケジュールパターンを出力するスケジュールパターン出力手段とを備えたことを特徴とする繰り返し型スケジュール生成装置。請求項1に記載の繰り返し型スケジュール生成装置において、前記流量制限値算出手段は、前記割当作業情報に規定された各プロトコールをjと表し、前記繰り返し期間中の日の各プロトコールjの患者数をdjと表し、前記割当時間情報に示される各日の集合をIと表し、前記集合Iに含まれる各日の個数を|I|と表し、前記流量制限値のうちの流量上限値をFLOWと表すとき、下記式に示すように、前記流量上限値FLOWを前記流量制限値として算出する流量上限値算出部を備え、前記使用枝数制限値算出手段は、前記割当時間情報に示される週の種類の数をhと表し、前記使用枝数制限値のうちの使用枝数上限値をUSEと表すとき、下記式に示すように、前記使用枝数上限値USEを前記使用枝数制限値として算出する使用枝数上限値算出部を備えたことを特徴とする繰り返し型スケジュール生成装置。請求項2に記載の繰り返し型スケジュール生成装置において、前記スケジュールパターンの出力後、ユーザの操作に応じて、前記流量制限値及び前記使用枝数制限値を個別に設定する制限値設定手段と、前記設定された流量制限値及び使用枝数制限値に基づいて、前記制限ネットワーク作成手段、前記スケジュールパターン生成手段及び前記スケジュールパターン出力手段を再実行するように制御する再実行制御手段とを備えたことを特徴とする繰り返し型スケジュール生成装置。請求項3に記載の繰り返し型スケジュール生成装置において、前記制限値設定手段は、前記流量制限値のうちの流量下限値をFLOWLと表すとき、下記式に示すように、前記流量下限値FLOWLを前記流量制限値として算出し、当該算出した流量下限値FLOWLを設定する流量下限値設定部と、前記使用枝数制限値のうちの使用枝数下限値をUSELと表すとき、下記式に示すように、前記使用枝数下限値USELを前記使用枝数制限値として算出し、当該算出した使用枝数下限値USELを設定する使用枝数下限値設定部とを備えたことを特徴とする繰り返し型スケジュール生成装置。請求項2に記載の繰り返し型スケジュール生成装置において、立案期間Tを示す立案期間情報を記憶する第6の記憶手段と、前記各治療工程の間に空ける日数種類lを規定した日数種類情報を記憶する第7の記憶手段と、前記各治療工程のうちの照射治療工程を休む日のパターンkを規定した休治療日パターン情報を記憶する第8の記憶手段とを更に備え、前記スケジュールパターン生成手段は、前記立案期間T中の日を{0,...,T}と表し、前記各治療工程をsと表し、当該各治療工程sの集合をSと表し、前記集合Iに含まれる各日をiと表し、前記パターンkの集合をKと表し、前記各プロトコールjの集合をJと表し、前記日数種類lの集合をLと表し、一週間中の曜日をwと表し、当該曜日wの集合をWと表し、前記繰り返し期間中の曜日wの日の集合をDAYと表し、前記各日iと前記パターンkと前記各プロトコールjと前記日数種類lとの組合せikjlにおける患者数をyikjlと表し、前記患者数yikjlのとりうる値の上界値をγと表し但し、r=max{FLOW|j∈J}、のとりうる値の上界値をυと表し但し、υ=max{USE|j∈J}、前記組合せikjlにおいて患者を一人割り当てた場合の日tにおける治療工程sの作業時間を示すコストをc’tsikjlと表し、前記組合せikjlにおいて患者を一人でも割り当てた場合を1とする0-1変数をxikjlと表し、前記曜日wにプロトコールjを割り当てた場合を1とする0-1変数をuwjと表し、整数全体をZと表し、下記式に示す混合整数線形計画問題の評価値をzと表すとき、前記評価値zを最小化するように前記混合整数線形計画問題を解くことにより、前記スケジュールパターンを生成するソルバ手段を備えたことを特徴とする繰り返し型スケジュール生成装置。請求項1に記載の繰り返し型スケジュール生成装置において、前記放射線治療に代えて、定期作業を用い、前記各プロトコールに代えて、各作業種類を用い、前記各治療工程に代えて、各作業工程を用い、前記患者数に代えて、作業対象数を用い、前記治療部位と照射治療工程の回数に代えて、点検作業又は工具の交換作業を用い、前記各治療工程間に空ける日数及び前記照射治療工程の回数に関する制約に代えて、前記作業工程間に空ける日数に関する制約を用い、前記曜日と週との組合せに代えて、作業点を用い、前記割当時間情報に代えて、作業点情報を用い、前記曜日種類に代えて、作業員を用い、前記割当時間種類情報に代えて、作業点種類情報を用い、前記割当時間ノードに代えて、作業点ノードを用い、前記プロトコールノードに代えて、作業種類ノードを用い、前記割当時間種類ノードに代えて、前記作業点種類ノードを用い、前記割当時間種類アークに代えて、作業点種類アークを用いたことを特徴とする繰り返し型スケジュール生成装置。情報入力手段、第1の記憶手段、第2の記憶手段、第3の記憶手段、第4の記憶手段、第5の記憶手段、割当ネットワーク作成手段、流量制限値算出手段、種類要素生成手段、使用枝数制限値算出手段、制限ネットワーク作成手段、スケジュールパターン生成手段、スケジュールパターン記憶手段、及びスケジュールパターン出力手段を備え、周期的に繰り返される繰り返し期間中の曜日と週との組合せ毎に、放射線治療の各プロトコールに対して基点となる治療工程の日を示すように当該各プロトコールを割当てたスケジュールパターンを生成する繰り返し型スケジュール生成装置が実行する繰り返し型スケジュール生成方法であって、前記情報入力手段が、前記各治療工程における治療部位と照射治療工程の回数とを示す前記各プロトコールと、前記各プロトコールの患者数とを規定した割当作業情報を前記第1の記憶手段に書込む工程と、前記情報入力手段が、前記各プロトコールに対する各治療工程及び当該各治療工程のうちの基点となる治療工程を規定した治療工程情報を前記第2の記憶手段に書込む工程と、前記情報入力手段が、前記各治療工程間に空ける日数及び前記照射治療工程の回数に関する制約を規定した制約情報を前記第3の記憶手段に書込む工程と、前記情報入力手段が、前記各プロトコールを割当て可能な曜日と週との組合せを個別に示す割当時間情報を前記第4の記憶手段に書込む工程と、前記情報入力手段が、前記割当時間情報内の曜日と週との組合せに個別に割り当てられた曜日種類を規定した割当時間種類情報を前記第5の記憶手段に書込む工程と、前記割当ネットワーク作成手段が、前記割当時間情報及び前記割当作業情報に基づいて、当該割当時間情報に示される曜日と週との組合せを個別に示す複数の割当時間ノードと、当該割当作業情報に規定された各プロトコールを個別に示す複数のプロトコールノードと、前記各割当時間ノードと前記各プロトコールノードとの組合せを個別に示す複数の割当アークと、前記各プロトコールの患者数とを含む割当ネットワークを作成する割当ネットワーク作成工程と、前記流量制限値算出手段が、前記割当作業情報及び前記割当時間情報に基づいて、前記各割当アーク毎に、割当て可能な患者数としての流量を制限する流量制限値を算出する流量制限値算出工程と、前記種類要素生成手段が、前記割当時間情報及び前記割当時間種類情報に基づいて、前記曜日種類を個別に示す複数の割当時間種類ノードと、前記各割当時間種類ノードのうちの有効な割当時間種類ノードを個別に示す複数の割当時間種類アークとを生成する種類要素生成工程と、前記使用枝数制限値算出手段が、前記割当作業情報、前記割当時間情報及び前記割当時間種類情報に基づいて、前記各プロトコール毎に、前記各割当時間種類アークの使用枝数を制限する使用枝数制限値を算出する使用枝数制限値算出工程と、前記制限ネットワーク作成手段が、前記割当ネットワーク、前記流量制限値、前記各割当時間種類ノード、前記各割当時間種類アーク及び前記使用枝数制限値を統合することにより、前記各割当時間ノードと、前記各プロトコールノードと、前記各割当アークと、前記患者数と、前記流量制限値と、前記各割当時間種類ノードと、前記各割当時間種類アークと、前記使用枝数制限値とを含む制限ネットワークを作成する制限ネットワーク作成工程と、前記スケジュールパターン生成手段が、前記制限ネットワーク、前記治療工程情報及び前記制約情報に基づいて、前記スケジュールパターンを生成するスケジュールパターン生成工程と、前記スケジュールパターン生成手段が、前記生成したスケジュールパターンを前記スケジュールパターン記憶手段に書込むスケジュールパターン書込工程と、前記スケジュールパターン出力手段が、前記スケジュールパターン記憶手段からスケジュールパターンを読出し、当該読出したスケジュールパターンを出力するスケジュールパターン出力工程とを備えたことを特徴とする繰り返し型スケジュール生成方法。請求項7に記載の繰り返し型スケジュール生成方法において、前記流量制限値算出工程は、前記割当作業情報に規定された各プロトコールをjと表し、前記繰り返し期間中の日の各プロトコールjの患者数をdjと表し、前記割当時間情報に示される各日の集合をIと表し、前記集合Iに含まれる各日の個数を|I|と表し、前記流量制限値のうちの流量上限値をFLOWと表すとき、下記式に示すように、前記流量上限値FLOWを前記流量制限値として算出する流量上限値算出工程を備え、前記使用枝数制限値算出工程は、前記割当時間情報に示される週の種類の数をhと表し、前記使用枝数制限値のうちの使用枝数上限値をUSEと表すとき、下記式に示すように、前記使用枝数上限値USEを前記使用枝数制限値として算出する使用枝数上限値算出工程を備えたことを特徴とする繰り返し型スケジュール生成装置。請求項8に記載の繰り返し型スケジュール生成方法において、前記繰り返し型スケジュール生成装置は、制限値設定手段及び再実行制御手段を更に備え、前記制限値設定手段が、前記スケジュールパターンの出力後、ユーザの操作に応じて、前記流量制限値及び前記使用枝数制限値を個別に設定する制限値設定工程と、前記再実行制御手段が、前記設定された流量制限値及び使用枝数制限値に基づいて、前記制限ネットワーク作成工程、前記スケジュールパターン生成工程、前記スケジュールパターン書込工程及び前記スケジュールパターン出力工程を再実行するように、前記制限ネットワーク作成手段、前記スケジュールパターン生成手段及び前記スケジュールパターン出力手段を制御する再実行制御工程とを備えたことを特徴とする繰り返し型スケジュール生成方法。請求項9に記載の繰り返し型スケジュール生成方法において、前記制限値設定工程は、前記流量制限値のうちの流量下限値をFLOWLと表すとき、下記式に示すように、前記流量下限値FLOWLを前記流量制限値として算出し、当該算出した流量下限値FLOWLを設定する流量下限値設定工程と、前記使用枝数制限値のうちの使用枝数下限値をUSELと表すとき、下記式に示すように、前記使用枝数下限値USELを前記使用枝数制限値として算出し、当該算出した使用枝数下限値USELを設定する使用枝数下限値設定工程とを備えたことを特徴とする繰り返し型スケジュール生成方法。請求項8に記載の繰り返し型スケジュール生成方法において、前記繰り返し型スケジュール生成装置は、第6の記憶手段、第7の記憶手段及び第8の記憶手段を更に備え、前記情報入力手段が、立案期間Tを示す立案期間情報を前記第6の記憶手段に書込む工程と、前記情報入力手段が、前記各治療工程の間に空ける日数種類lを規定した日数種類情報を前記第7の記憶手段に書込む工程と、前記情報入力手段が、前記各治療工程のうちの照射治療工程を休む日のパターンkを規定した休治療日パターン情報を前記第8の記憶手段に書込む工程とを更に備え、前記スケジュールパターン生成工程は、前記立案期間T中の日を{0,...,T}と表し、前記各治療工程をsと表し、当該各治療工程sの集合をSと表し、前記集合Iに含まれる各日をiと表し、前記パターンkの集合をKと表し、前記各プロトコールjの集合をJと表し、前記日数種類lの集合をLと表し、一週間中の曜日をwと表し、当該曜日wの集合をWと表し、前記繰り返し期間中の曜日wの日の集合をDAYと表し、前記各日iと前記パターンkと前記各プロトコールjと前記日数種類lとの組合せikjlにおける患者数をyikjlと表し、前記患者数yikjlのとりうる値の上界値をγと表し但し、r=max{FLOW|j∈J}、のとりうる値の上界値をυと表し但し、υ=max{USE|j∈J}、前記組合せikjlにおいて患者を一人割り当てた場合の日tにおける治療工程sの作業時間を示すコストをc’tsikjlと表し、前記組合せikjlにおいて患者を一人でも割り当てた場合を1とする0-1変数をxikjlと表し、前記曜日wにプロトコールjを割り当てた場合を1とする0-1変数をuwjと表し、整数全体をZと表し、下記式に示す混合整数線形計画問題の評価値をzと表すとき、前記評価値zを最小化するように前記混合整数線形計画問題を解くことにより、前記スケジュールパターンを生成するソルバ工程を備えたことを特徴とする繰り返し型スケジュール生成方法。請求項7に記載の繰り返し型スケジュール生成方法において、前記放射線治療に代えて、定期作業を用い、前記各プロトコールに代えて、各作業種類を用い、前記各治療工程に代えて、各作業工程を用い、前記患者数に代えて、作業対象数を用い、前記治療部位と照射治療工程の回数に代えて、点検作業又は工具の交換作業を用い、前記各治療工程間に空ける日数及び前記照射治療工程の回数に関する制約に代えて、前記作業工程間に空ける日数に関する制約を用い、前記曜日と週との組合せに代えて、作業点を用い、前記割当時間情報に代えて、作業点情報を用い、前記曜日種類に代えて、作業員を用い、前記割当時間種類情報に代えて、作業点種類情報を用い、前記割当時間ノードに代えて、作業点ノードを用い、前記プロトコールノードに代えて、作業種類ノードを用い、前記割当時間種類ノードに代えて、前記作業点種類ノードを用い、前記割当時間種類アークに代えて、作業点種類アークを用いたことを特徴とする繰り返し型スケジュール生成方法。以上説明したように、第1の先行技術では、作業が一週間単位で入るとは限らない場合に、作業員が認識しやすい規則性を持つスケジュールを作成することが困難である。また、第2の先行技術では、仮に、同じ作業員が担当する作業点を自由に与えたとしても、作業点における作業量が偏らないように平準化しつつ、規則性を持って特定の作業が集まりやすいスケジュールを、容易に作成することが困難である。本発明は上記実情を考慮してなされたもので、規則性のある偏らないスケジュールパターンを生成し得る繰り返し型スケジュール生成装置及び方法を提供することを目的とする。本発明の一つの局面は、周期的に繰り返される繰り返し期間中の曜日と週との組合せ毎に、放射線治療の各プロトコールに対して基点となる治療工程の日を示すように当該各プロトコールを割当てたスケジュールパターンを生成する繰り返し型スケジュール生成装置であって、前記各治療工程における治療部位と照射治療工程の回数とを示す前記各プロトコールと、前記各プロトコールの患者数とを規定した割当作業情報を記憶する第1の記憶手段と、前記各プロトコールに対する各治療工程及び当該各治療工程のうちの基点となる治療工程を規定した治療工程情報を記憶する第2の記憶手段と、前記各治療工程間に空ける日数及び前記照射治療工程の回数に関する制約を規定した制約情報とを記憶する第3の記憶手段と、前記各プロトコールを割当て可能な曜日と週との組合せを個別に示す割当時間情報を記憶する第4の記憶手段と、前記割当時間情報内の曜日と週との組合せに個別に割り当てられた曜日種類を規定した割当時間種類情報を記憶する第5の記憶手段と、前記割当時間情報及び前記割当作業情報に基づいて、当該割当時間情報に示される曜日と週との組合せを個別に示す複数の割当時間ノードと、当該割当作業情報に規定された各プロトコールを個別に示す複数のプロトコールノードと、前記各割当時間ノードと前記各プロトコールノードとの組合せを個別に示す複数の割当アークと、前記各プロトコールの患者数とを含む割当ネットワークを作成する割当ネットワーク作成手段と、前記割当作業情報及び前記割当時間情報に基づいて、前記各割当アーク毎に、割当て可能な患者数としての流量を制限する流量制限値を算出する流量制限値算出手段と、前記割当時間情報及び前記割当時間種類情報に基づいて、前記曜日種類を個別に示す複数の割当時間種類ノードと、前記各割当時間種類ノードのうちの有効な割当時間種類ノードを個別に示す複数の割当時間種類アークとを生成する種類要素生成手段と、前記割当作業情報、前記割当時間情報及び前記割当時間種類情報に基づいて、前記各プロトコール毎に、前記各割当時間種類アークの使用枝数を制限する使用枝数制限値を算出する使用枝数制限値算出手段と、前記割当ネットワーク、前記流量制限値、前記各割当時間種類ノード、前記各割当時間種類アーク及び前記使用枝数制限値を統合することにより、前記各割当時間ノードと、前記各プロトコールノードと、前記各割当アークと、前記患者数と、前記流量制限値と、前記各割当時間種類ノードと、前記各割当時間種類アークと、前記使用枝数制限値とを含む制限ネットワークを作成する制限ネットワーク作成手段と、前記制限ネットワーク、前記治療工程情報及び前記制約情報に基づいて、前記スケジュールパターンを生成し、当該生成したスケジュールパターンをスケジュールパターン記憶手段に書込むスケジュールパターン生成手段と、前記書込まれたスケジュールパターンを記憶する前記スケジュールパターン記憶手段と、前記スケジュールパターン記憶手段からスケジュールパターンを読出し、当該読出したスケジュールパターンを出力するスケジュールパターン出力手段とを備えた繰り返し型スケジュール生成装置である。以上説明したように本発明によれば、流量制限値により、特定の曜日と週との組合せに対して割り当てられる患者数を制限すると共に、使用枝数制限値により、曜日種類を制限する構成により、規則性のある偏らないスケジュールパターンを生成できる。以下、各実施形態について図面を用いて説明するが、その前に各実施形態の概要を述べる。第1及び第2の実施形態では、放射線治療のスケジューリング問題を、ネットワークフロー問題として捉え、例えば、その割当ネットワークに割当時間の種類という特別なノードとアークを追加し、各作業種類でこのアークを使用する本数に制限をかける構成により、生成されるスケジュールパターンに規則性を与えている。本数制限は、上限だけをゆるくかけることもできるし、上下限を設けることで厳しく制約することも可能であって、容易に調節できる。使用枝数の制約は、混合整数線形計画問題の制約式として定式化できるため、最適化が可能である。また、割当作業を割当時間に割り当てる割当アークを流れる流量に制限をかける流量制限値についても、制約として記述することで、目的関数例えば、平準化を最良化する最適化が可能となる。制約に用いる上下限値の算出の仕方は、第1及び第2の実施形態で一例を示している。また、第1及び第2の実施形態においては、放射線治療スケジューリングにおける最適化問題のMILPmixedintegerlinierprogramming混合整数線形計画法定式化の例を述べている。第3の実施形態では、第1の実施形態の変形例として、割当時間種類のグルーピングを作業員ごとに分けると、各作業員が、特定の作業種類の作業を繰返し行うスケジュールパターンを生成できることについて述べている。第3の実施形態に述べた例では、割当時間種類ノードに代えて作業点種類ノードを用い、割当時間ノードに代えて作業点ノードを用いている等の変更点がある。第3の実施形態は、情報の名称等が第1の実施形態から変更されているが、情報を加工又は演算する処理の内容については、第1の実施形態と同様となっている。以上が各実施形態の概要である。続いて、各実施形態について具体的に説明する。<第1の実施形態>図1は第1の実施形態に係る繰り返し型スケジュール生成装置の構成を示すブロック図であり、図2はこの繰り返し型スケジュール生成装置により生成されるスケジュールパターンの一例を示す模式図である。繰り返し型スケジュール生成装置は、各部1~12の処理により、周期的に繰り返される繰り返し期間のスケジュールパターンを生成する。繰り返し型スケジューリングにおけるスケジュールパターンとは、割当作業と割当時間を有し、所定数の各割当作業が割当時間のいずれかに割当てられたものを意味する。この種のスケジュールパターンにおいては、例えば、図2に示すように、繰り返し期間中の曜日と週との組合せ割当時間毎に、放射線治療の各プロトコール病状種別に対して基点となる治療工程割当作業の日を示すように当該各プロトコールが割当てられている。図2中、月1は第1週の月曜日を表す略称である。肺1は、治療部位が肺であり、照射治療工程が1回である旨のプロトコールを表す略称である。同様に、肝8及び前16は、それぞれ治療部位が肝臓及び前立腺であり、照射治療工程が8回及び16回である旨の各プロトコールを表す略称である。以下の説明では、これらの略称を用いる。また、スケジュールパターンは、各プロトコールを当該各プロトコールの患者数だけ含んでいる。スケジュールパターンに10個の肝8が割り当てられていることは、プロトコール肝8の患者数が10人であることに対応している。一方、繰り返し型スケジュール生成装置は、図1に示すように、情報入力部1、割当作業情報記憶部2、割当時間情報記憶部3、割当時間種類情報記憶部4、スケジュール情報記憶部5、割当ネットワーク作成部6、流量制限値算出部7、種類要素生成部8、使用枝数制限値算出部9、制限ネットワーク作成部10、スケジュールパターン生成部11及びスケジュールパターン記憶及び出力部12を備えている。情報入力部1は、ユーザと各部2~12との間の入力インタフェースであり、例えば、ユーザの操作に応じて、種々の情報を各記憶部2~5に書込む機能をもっている。各記憶部2~5は、各部1,6~12から読出/書込可能な記憶装置又は記憶装置の記憶領域であり、例えば、情報入力部1から書込まれた情報を各部6~12から読出可能に記憶する。また例えば、各記憶部2~5は、各部6~11の処理途中のデータを記憶してもよい。具体的には、割当作業情報記憶部第1の記憶部2は、図3に示すように、各治療工程における治療部位と照射治療工程の回数とを示す各プロトコールjと、各プロトコールjの患者数djとを規定した割当作業情報を記憶する。割当時間情報記憶部第4の記憶部3は、図4に示すように、各プロトコールjを割当て可能な曜日と週との組合せを個別に示す割当時間情報iを記憶する。割当時間種類情報記憶部第5の記憶手段4は、図5に示すように、割当時間情報i内の曜日と週との組合せに個別に割り当てられた曜日種類wを規定した割当時間種類情報を記憶する。スケジュール情報記憶部第2、第6~第8の記憶部5は、図6~図11に示すように、情報入力部1から書き込まれた立案期間情報、プロトコール情報、治療工程情報、制約情報、休治療日パターン情報及び日数種類情報を記憶する。なお、スケジュール情報記憶部5内の各情報のうち、例えば、立案期間情報、プロトコール情報、休治療日パターン情報及び日数種類情報は、スケジュールパターンを作成する観点からは、任意の付加的事項であって省略可能である。但し、立案期間情報、休治療日パターン情報及び日数種類情報は、後述する混合整数線形計画問題を定式化する観点から省略しない方が好ましい。また、各情報は、説明の便宜上、スケジュール情報記憶部5にまとめて記憶しているが、適宜、他の記憶部に記憶させるように変形してもよい。立案期間情報は、図6に示すように、立案期間Tを示す情報である。プロトコール情報は、図7に示すように、各プロトコールjと照射治療工程の回数とを規定した情報である。治療工程情報は、図8に示すように、各治療工程sを規定した情報であり、例えば、各プロトコールjに対する各治療工程s及び当該各治療工程sのうちの基点となる治療工程以下、基準治療工程ともいうを規定した情報である。なお、基準治療工程では、各日の基準治療工程を実施する作業時間の上限が指定可能となっている。他の治療工程は、基準治療工程からの相対位置間に空ける日数でスケジュールが生成されることが規則で定められている。基準治療工程としては、最初の治療工程である固定具作成に限らず、CTシミュレーションや第1回照射を指定することも可能である。制約情報は、図9に示すように、治療工程の制約を規定した情報であり、具体的には、各治療工程間に空ける日数及び照射治療工程の回数に関する制約を規定した情報である。休治療日パターン情報は、図10に示すように、各治療工程のうちの照射治療工程を休む日図中、照射の休治療日のパターンkを規定した情報である。日数種類情報は、図11に示すように、各治療工程の間に空ける日数種類lを規定した情報である。割当ネットワーク作成部6は、各記憶部2,3内の割当時間情報及び割当作業情報に基づいて、割当ネットワークを作成する機能をもっている。割当ネットワークは、図12に示すように、当該割当時間情報iに示される曜日と週との組合せを個別に示す複数の割当時間ノード月1,...,金2と、当該割当作業情報に規定された各プロトコールjを個別に示す複数のプロトコールノード肺1,...,前16と、各割当時間ノードと各プロトコールノードとの組合せを個別に示す複数の割当アークαと、各プロトコールjの患者数djとを含んでいる。なお、割当ネットワークは、図13に示すように、休治療日パターンkを個別に示す休治療日パターンノードと、間に空ける日数種類lを個別に示す日数種類ノードとを更に含んでもよい。本実施形態では、主に、図13に示す割当ネットワークの場合を例に挙げて述べる。流量制限値算出部7は、各記憶部3,5内の割当作業情報及び割当時間情報iに基づいて、各割当アークα毎に、割当て可能な患者数djとしての流量を制限する流量制限値を算出する流量制限値算出機能をもっている。なお、流量制限値は、図2中、同一の日同一の列における同一プロトコールの個数例、木1の日の肝8の個数は2の上限値又は上下限値と解釈することも可能である。ここで、流量制限値算出機能は、割当作業情報に規定された各プロトコールをjと表し、繰り返し期間中の日の各プロトコールjの患者数をdjと表し、割当時間情報iに示される各日の集合をIと表し、この集合Iに含まれる各日の個数を|I|と表し、流量制限値のうちの流量上限値をFLOWと表すとき、下記式1に示すように、流量上限値FLOWを流量制限値として算出する流量上限値算出機能を備えてもよい。種類要素生成部8は、各記憶部3,4内の割当時間情報i及び割当時間種類情報に基づいて、図14の一部に示す如き、曜日種類wを個別に示す複数の割当時間種類ノード月,...,金と、各割当時間種類ノード月,...,金のうちの有効な割当時間種類ノードを個別に示す複数の割当時間種類アークβとを生成する機能をもっている。なお、種類要素の語は割当時間種類ノード及び割当時間種類アークと読み替えてもよい。使用枝数制限値算出部9は、割当作業情報、割当時間情報i及び割当時間種類情報に基づいて、各プロトコールj毎に、各割当時間種類アークの使用枝数を制限する使用枝数制限値を算出する使用枝数制限値算出機能をもっている。なお、使用枝数制限値は、図2中、曜日及び週の組合せ月1,...,金2に関し、有効な曜日の種類月、火、木、金の個数図2では4の上限値又は上下限値としても解釈可能である。ここで、使用枝数制限値算出機能は、割当時間情報iに示される週の種類の数をhと表し、使用枝数制限値のうちの使用枝数上限値をUSEと表すとき、下記式2に示すように、使用枝数上限値USEを使用枝数制限値として算出する使用枝数上限値算出機能を備えてもよい。制限ネットワーク作成部10は、割当ネットワーク、流量制限値、各割当時間種類ノード、各割当時間種類アーク及び使用枝数制限値を統合することにより、制限ネットワークを作成する機能をもっている。制限ネットワークは、図14に示すように、各割当時間ノード月1,...,金2と、各プロトコールノード肺1,...,前16と、各割当アークαと、患者数djと、流量制限値図示せずと、各割当時間種類ノード月,...,金と、各割当時間種類アークβと、使用枝数制限値図示せずとを含んでいる。なお、制限ネットワークは、図15に示すように、休治療日パターンkを個別に示す休治療日パターンノードと、間に空ける日数種類lを個別に示す日数種類ノードとを更に含んでもよい。本実施形態では、主に、図15に示す制限ネットワークの場合を例に挙げて述べる。また、制限ネットワークの語は、割当時間種類をもつネットワークと読み替えてもよい。スケジュールパターン生成部11は、制限ネットワーク、治療工程情報及び制約情報に基づいて、スケジュールパターンを生成するスケジュールパターン生成機能と、当該生成したスケジュールパターンをスケジュールパターン記憶及び出力部12に書込むスケジュールパターン書込機能とをもっている。ここで、スケジュールパターン生成機能は、立案期間T中の日を{0,...,T}と表し、各治療工程をsと表し、当該各治療工程sの集合をSと表し、集合Iに含まれる各日をiと表し、パターンkの集合をKと表し、各プロトコールjの集合をJと表し、日数種類lの集合をLと表し、一週間中の曜日をwと表し、当該曜日wの集合をWと表し、繰り返し期間中の曜日wの日の集合をDAYと表し、各日iとパターンkと各プロトコールjと日数種類lとの組合せikjlにおける患者数をyikjlと表し、患者数yikjlのとりうる値の上界値をγと表し但し、r=max{FLOW|j∈J}、のとりうる値の上界値をυと表し但し、υ=max{USE|j∈J}、組合せikjlにおいて患者を一人割り当てた場合の日tにおける治療工程sの作業時間を示すコストをc’tsikjlと表し、組合せikjlにおいて患者を一人でも割り当てた場合を1とする0-1変数をxikjlと表し、曜日wにプロトコールjを割り当てた場合を1とする0-1変数をuwjと表し、整数全体をZと表し、下記式3に示す混合整数線形計画問題の評価値をzと表すとき、評価値zを最小化するように混合整数線形計画問題を解くことにより、スケジュールパターンを生成するソルバ機能を備えてもよい。スケジュールパターン記憶及び出力部12は、スケジュールパターン生成部11から書込まれたスケジュールパターンを記憶するスケジュールパターン記憶部図示せずを有し、当該スケジュールパターン記憶部からスケジュールパターンを読出し、当該読出したスケジュールパターンを出力する機能をもっている。次に、以上のように構成された繰り返し型スケジュール生成装置の動作について図16を参照しながら説明する。以下の説明は、重粒子線治療における2週間繰返しのスケジュールを作成する場合を例に挙げて述べる。始めに、重粒子線におけるスケジューリングの概要について述べる。ここでのスケジューリングとは、来院患者の頻度と患者のプロトコール病状種別の割合から、立案期間分の予定枠を作成するものである。病院では、日々の負荷を超えないように、例えば半期分の予定枠治療枠を予め作成しておき、新患の予約時に枠を埋めていくという方法を取る場合が多い。治療枠は、他の療法との併用で数ヶ月先に予約を入れる場合や、来院患者数を見積もる場合にも有用である。予定枠を作成して埋めていく方法は、スタッフにとって理解しやすい規則性と類似のスケジュールが、週単位や2週間単位に繰り返される性質と、負荷平準化による効率性とを同時に満たすことで、より多くの患者を安全に治療することを図っている。重粒子線治療では図17に示す治療工程があり、このうちスケジューリングを行う治療工程は、例えば固定具作成、CTシミュレーション、照射の3つとする。他に、リハーサルなどの工程をスケジューリングに含む場合もある。照射治療工程は複数回行われ、プロトコールの種類によって回数が異なる。スケジューリングを行う治療工程について、患者1人分の実施日スケジュールがきまったものを、1つのスケジュール枠以下、枠と呼ぶ。半期分のすべての患者の枠の集合を、ここでは治療枠と呼ぶ。各枠はプロトコールが決まっているが、どの患者に割当てるかについては決まっていない。重粒子線治療における枠作成問題以下、枠作成問題とは、部屋の数などの施設の規模あるいは病院スタッフの人数、各種スケジュールなどによって各プロトコールにおける各治療工程の実行可能な条件と最適化目的が与えられているとき、実行可能な条件を満たし最適化目的を良くする、各患者における枠の集合治療予定枠を得ることである。安全性と利便性の観点から、治療枠には週単位の規則性を持たせることが好ましい。実際には週1回未満のプロトコールもあるため、頻度が最小のプロトコールに合わせて繰返し期間を設定し、立案期間中はこの繰返し期間のスケジュールが周期的に繰り返されることとしている。さらに、繰返し期間中のスケジュールについても、週単位の規則性及び、人数の多いプロトコールは各日でなるべく均等に割り当てる規則性が重視される。一方、スタッフの負荷を軽減するためには、各日の作業時間がなるべく平準化されるのが望ましい。祝日や機器のメンテナンス等で不規則なケースが発生することから、作業時間は立案期間の全体で考慮する必要がある。立案期間中で照射の作業にかかる作業時間が最も多い日の作業時間を最小にする。または、繰り返し期間中の各日ごとに割り当てられた各治療工程の作業時間を合計した合計作業時間のうち、最大の合計作業時間を持つ日の当該合計作業時間を最小にするようにしてもよい。工程間に空ける日数が与えられるため、繰返し期間中の各日に各プロトコールの治療開始固定具作成を1人割り当てたときの立案期間中の各日の照射にかかる作業時間を予め計算できる。治療予定枠は、繰り返し期間中の各プロトコールの開始日の指定を含むスケジュールパターンを決定することで生成できる。このようなスケジューリングを行うにあたって、情報入力部1は、ユーザの操作に応じて割当作業情報、割当時間情報i及び割当時間種類情報をそれぞれ割当作業情報記憶部2、割当時間情報記憶部3及び割当時間種類情報記憶部4に書込む。同様に、情報入力部1は、ユーザの操作に応じて、立案期間情報、プロトコール情報、治療工程情報、制約情報、休治療日パターン情報及び日数種類情報をスケジュール情報記憶部5に書込む。続いて、割当ネットワーク作成部6は、各記憶部2,3内の割当時間情報i及び割当作業情報に基づいて、当該割当時間情報iに示される曜日と週との組合せを個別に示す複数の割当時間ノード月1,...,金2と、当該割当作業情報に規定された各プロトコールjを個別に示す複数のプロトコールノード肺1,...,前16と、各割当時間ノードと各プロトコールノードとの組合せを個別に示す複数の割当アークαと、各プロトコールjの患者数djとを含む割当ネットワークを作成するS11。このステップS11の詳細については後述する。割当ネットワークは、例えば図12に示したように、どの割当時間にどの割当作業をいくつずつ割り当てるかを二部グラフのマッチングとして表現するもので、割当時間と割当作業間のアークである割当アークαに流れる流量の組が、ひとつのスケジュールパターン要素を表す。割当ネットワークは、図13に示したように、各割当時間と各割当作業に、それぞれ休治療日パターンkと間に空ける日数種類lの選択肢を拡張したネットワークでも良い。間に空ける日数種類lの選択肢とは、基準治療工程以外の治療工程の作業日が一意に定まらなくなる問題を解決するため、治療工程間に空ける日数種類の中で最適化により良い日数を選択可能とするものである。休治療日パターンkの選択肢は、1週間内の照射可能日が患者一人あたりの1週間の照射回数上限よりも多い週がある場合に、複数通りの照射日の入れ方のパターンの中から最適化により良いパターンを選択可能とするものである。また、日数種類lと休治療日パターンkとの2点を選択可能としたことに加え、同じ開始日iの同じプロトコールjの各患者に対し、図18に示すように、治療工程間に空ける日数種類lと、休治療日パターンkとの組合せを1種類にしたノードkとノードl間の割当アークαを1本だけ使用する。これにより、基準治療工程以外の治療工程の作業日が一意に定まらなくなる問題を解決している。流量制限値算出部7は、割当ネットワークの各割当アークαと各割当作業pの各組αpについて、流量制限値を算出するS12、S13。各割当作業pは、各プロトコールjに相当する。流量制限は、割当てが特定の割当アークαに偏らないように制限をかけるものであり、例えば、上限値が下記のように算出可能となっている。すなわち、流量上限値FLOWは、プロトコールjを繰り返し期間中の日1日あたりに割り当ててよい上限を示すため、集合Iの要素数を|I|とすると、前述した式1に基づいて、プロトコールjの患者数djを繰り返し期間中の日数|I|で除算し、この除算結果以上の最小の整数として算出できる。なお、除算結果以上の最小の整数は、除算結果を下回らない最小の整数又は除算結果の小数点以下を切り上げた整数と読み替えてもよい。ここで、図3及び図4に示す情報に基づく流量上限値FLOWの算出例を図19に示す。例えば、プロトコールjが肝8の場合、dj=10、|I|=8、として式1に基づき、FLOWは、dj/|I|=10/8以上の最小の整数2として算出される。次に、種類要素生成部8は、各記憶部3,4内の割当時間情報及び割当時間種類情報に基づいて、図15に示した如き、割当時間種類ノード月,...,金を生成すると共にS14、各割当時間種類ノードとソースノード、各割当時間種類ノードと割当時間種類に対応する割当時間ノード間にアークを生成するS15,S16。使用枝数制限値算出部9は、生成した各割当時間種類アークβ図20と各割当作業pの各組β,pについて、使用枝数制限値を算出するS17,S18。使用枝数制限値は、割当作業pが割り当てられる割当時間の種類数割当時間種類アークβの本数を制限するもので、例えば上限値が下記のように算出可能となっている。すなわち、使用枝数上限値USEは、繰返し期間中の週の数をhとし、前述した式2に基づいて、プロトコールjの患者数djを繰り返し期間中の週の数hで除算し、この除算結果以上の最小の整数として算出できる。図3~図5に示す情報に基づく使用枝数上限値USEの算出例を図21に示す。例えば、プロトコールjが肝8の場合、dj=10、h=2、として式2に基づき、USEは、dj/h=10/2以上の最小の整数5として算出されるが、割当時間種類アークβの使用枝数が最大で4本のため、算出結果を4としている。制限ネットワーク作成部10は、割当ネットワーク、流量上限値FLOW、各割当時間種類ノード月,...,金、各割当時間種類アークβ及び使用枝数上限値USEを統合することにより、図15に示した如き、制限ネットワークを作成するS19。この制限ネットワークは、混合整数線形計画問題として上記式3に示したように定式化でき、CPLEXなどの混合整数線形計画問題を解く汎用ソルバや、厳密解法、ヒューリスティクス解法、メタヒューリスティクス解法、によって最適化できる。定式化においては、入力データに関する記号を以下に定義する。{0,...,T}:立案期間中の日。S:治療工程の集合。I:繰り返し期間中の日の集合。K:照射パターンの集合。J:プロトコールの集合。L:間に空ける日数の集合。W:一週間中の曜日の集合。dj:繰り返し期間中の日のプロトコールjの患者数。次に、入力データから算出する値に関する記号を定義する。DAY:曜日がwである繰り返し期間中の日の集合を返す関数。USE:プロトコールjが取ってよい曜日の種類数の上限値を返す関数。FLOW:プロトコールjを繰り返し期間中の日の1日あたりに割り当ててよい上限量を返す関数。γ:yikjlのとりうる値の上界値。c’tsikjl:ikjlに患者を一人割り当てたときの、日tにおける治療工程sのコスト以下に、定式化中の変数に関する記号を定義する。yikjl:ikjlに患者を割り当てた人数を示す、整数変数。xikjl:ikjlに患者が1人でも割り当てられた場合に1を取る、0-1変数。uwj:曜日wにプロトコールjを割り当てたときに1を取る、0-1変数。ただし、Zは整数全体を示す記号とする。以上のような定義を用い、制限ネットワークは、式3に示した如き、混合整数線形計画問題として定式化される。スケジュールパターン生成部11は、このように定式化可能な制限ネットワークと、スケジュール記憶部5内の治療工程情報及び制約情報に基づいて、前述した式3の評価値zを最小化するように混合整数線形計画問題を解くことにより、図2に示した如き、スケジュールパターンを生成するS20。しかる後、スケジュールパターン生成部11は、生成したスケジュールパターンをスケジュールパターン記憶及び出力部12に書込む。スケジュールパターン記憶及び出力部12は、書込まれたスケジュールパターンをスケジュールパターン記憶部図示せずが記憶し、当該スケジュールパターン記憶部からスケジュールパターンを読出し、当該読出したスケジュールパターンを出力する。次に、前述したステップS11の割当ネットワークの作成処理について詳細に述べる。放射線治療の例では、割当時間は例えば固定具作成などの基準となる治療工程の開始日を示し、割当作業は、プロトコールの種類を示す。図13に示した割当ネットワークは、割当作業情報記憶部2内の割当作業情報と、割当時間情報記憶部3内の割当時間情報iと、スケジュール情報記憶部5内の各情報とに基づいて、割当ネットワーク作成部6により作成される。このような割当ネットワークの作成処理について図22を用いて詳細に説明する。始めに、割当ネットワーク作成部6は、シンクsinkノードとソースsourceノードという、ダミーのノードを1つずつ生成するS11-1。割当ネットワーク作成部6は、割当作業情報記憶部2内の割当作業情報に基づいて、各プロトコールjの割当作業ノード肺1,...,前16を生成するS11-2。続いて、割当ネットワーク作成部6は、以下のステップS11-4~S11-6の処理を全ての割当作業ノード肺1,...,前16について繰り返すS11-3。すなわち、割当ネットワーク作成部6は、各割当作業ノード肺1,...,前16について、割当作業ノードとシンクノード間のアークを生成しS11-4、スケジュール記憶部5内の日数種類情報に基づいて、間に空ける日数種類ノードを種類数分生成するS11-5。また、割当ネットワーク作成部6は、日数種類ノードと割当作業ノード間のアークを生成するS11-6。ステップS11-4~S11-6を繰り返すステップS11-3の終了後、割当ネットワーク作成部6は、割当時間情報記憶部3内の割当時間情報iに基づいて、割当時間ノード月1,...,金2を生成するS11-7。しかる後、前述同様に、以下のステップS11-9~S11-11の処理を全ての割当時間ノード月1,...,金2について繰り返すS11-8。すなわち、割当ネットワーク作成部6は、各割当時間ノード月1,...,金2について、割当時間ノードとソースノード間のアークを生成しS11-9、スケジュール記憶部5内の休治療日パターン情報に基づいて、休治療日パターンノードを生成するS11-10。また、割当ネットワーク作成部6は、休治療日パターンノードと割当時間ノード間のアークを生成するS11-11。ステップS11-9~S11-11を繰り返すステップS11-8の終了後、割当ネットワーク作成部6は、割当アークαを生成するS11-12。割当アークαの生成処理としては、図13に示す割当ネットワークの場合、各割当作業に関する間に空ける日数種類ノードと、各割当時間に関する休治療日パターンノードで完全2部グラフを構成すればよい。図12に示す割当ネットワークのように休治療日パターンノードと、間に空ける日数種類ノードとがない場合には、各割当作業ノード肺1,...,前16と各割当時間ノード月1,...,金2で完全2部グラフを構成すればよい。割当作業に関するノードと割当時間に関するノードで完全2部グラフを構成するように生成されたアークを割当アークαと呼ぶ。割当ネットワーク作成部6は、各割当アークαについて、スケジュール記憶部5内の各情報に基づいて、図23に示すように、多次元コストを算出しS11-13,S11-14、割当ネットワークを作成する。ここで、ステップS11-13の多次元コストの算出方法について詳細に説明する。繰り返し期間中の日i、プロトコールj、休治療日パターンk、間に空ける日数種類lに患者を1人割り当てた場合の、立案期間中の日tにおけるコストc’tsikjlは、次のように計算できる。始めに、繰り返し期間中の日iに対して、治療の開始日の集合Qiを以下で定義する。Qi={t|TIME=it∈1T}次に、各開始日tで開始される1人分の枠を下記の手順[1]~[3]で計算する。[1]開始日tに基準治療工程を設定する。[2]残りの治療工程の開始日は基準治療工程からの間に空ける日数より算出する。[3]照射治療工程の開始日から、照射回数プロトコールjによって決まっている分の照射治療工程を各日に設定する。照射治療工程の設定は各施設で定める照射のスケジュールの割当て方法に従って行う。例えば下記の割当て方法が使用可能である。もし、その週に照射が可能な日数が、一人当たりの週の照射回数上限値より多く、かつ、残りの照射回数が一人当たりの週の照射回数上限値以上の場合、照射治療工程は、休治療日パターンkで定められた割当て方法で設定する。それ以外の場合、照射治療工程の開始日から、翌日以降の照射可能な日に順次照射治療工程を割当てていく。各t∈Qiの開始日について、それぞれ1人分の枠を計算し、各立案期間中の日で照射治療工程が設定された人数を算出する。c’tsikjlは、コストを算出するときに固定具作成とCTシミュレーションについても各日の作業時間を算出した値を用い、プロトコールjの照射の作業時間をかけて計算できる図23。照射等の作業が複数日にわたる治療工程の作業時間は、回数によって作業時間を変えて設定しても良く、上記の計算方法と同様にして一人当たりの多次元コストを算出できる。上述したように本実施形態によれば、割当ネットワークの各割当アークα毎に、割当て可能な患者数としての流量を制限する流量制限値を算出し、曜日種類wを個別に示す各割当時間種類ノード月,...,金のうちの有効な割当時間種類ノードを個別に示す各割当時間種類アークβの使用枝数を各プロトコールj毎に算出し、割当ネットワーク、流量制限値、各割当時間種類ノード月,...,金、各割当時間種類アークβ及び使用枝数制限値を統合することにより、制限ネットワークを作成し、制限ネットワーク、治療工程情報及び制約情報に基づいて、スケジュールパターンを生成する構成により、規則性のある偏らないスケジュールパターンを生成することができる。補足すると、本実施形態では、割当時間種類ノード月,...,金および割当時間種類アークβを生成して、各割当時間種類アークβに使用枝数制限をつけることにより、例えば数週間単位のスケジュールパターンを生成するに当たって、週単位の規則性を容易に作成できるという効果を持つ。元のスケジュールの制約から前の週と次の週は必ずしも同じにはならないが、その各週が曜日によって行う作業がなるべくそろうスケジュールが生成できれば、作業員がスケジュールパターンを覚えやすく、作業員の負担が減って安全性が増す。また、複数か月単位の大まかな計画は上記先行技術で作成し、1日、あるいは単位時間単位の詳細な計画を本実施形態の繰り返し型スケジュール生成装置により作成しても良い。また、作業が一週間単位で必ず入るわけではない場合にも、作業員が認識しやすい規則性を持つスケジュールを作成できることに加え、ユーザによる目的関数のウェイト調節などが不要であるため、規則性の持たせ方を容易に調節できる図24、図25。例えば、流量上限値FLOW及び使用枝数上限値USEが無い場合、スケジュールパターンは、図24aに示すように、同一のプロトコール0の基準治療工程が、前の週と次の週とでは異なる曜日に割当てられており、曜日の規則性が無い。また、他の同一のプロトコール2又は7についても同様に規則性が無い。一方、流量上限値FLOW及び使用枝数上限値USEがある場合、スケジュールパターンは、図24bに示すように、同一のプロトコール0の基準治療工程が、前の週と次の週とで同一の曜日に割当てられており、曜日の規則性がある。他の同一のプロトコール2又は7についても同様に規則性がある。また例えば、流量上限値FLOW及び使用枝数上限値USEが無い場合で日によって割当て人数が異なるとき、スケジュールパターンは、図25aに示すように、同一のプロトコール0,4火2-木2又は7の基準治療工程について、曜日の規則性が無い。但し、他の同一のプロトコール4月1-月2、金1-金2については曜日の規則性がある。また他方、流量上限値FLOW及び使用枝数上限値USEがある場合で日によって割当て人数が異なるとき、スケジュールパターンは、図25bに示すように、同一のプロトコール0の基準治療工程について、曜日の規則性がある。他の同一のプロトコール4又は7についても同様に規則性がある。なお、本実施形態においては、図13に示したklの選択肢をもつ割当ネットワークに基づいて図15に示した制限ネットワークを混合整数線形計画問題として定式化した場合を例に挙げて説明したが、これに限らず、図12に示したklの選択肢をもたない割当ネットワークに基づいて図14に示す制限ネットワークを混合整数線形計画問題として定式化した場合としても実施できる。この場合、休治療日パターンk及び日数種類lを予め定めておけばよい。このことは以下の各実施形態でも同様である。<第2の実施形態>図26は第2の実施形態に係る繰り返し型スケジュール生成装置の構成を示すブロック図であり、図1と同一部分には同一符号を付してその詳しい説明を省略し、ここでは異なる部分について主に述べる。なお、以下の各実施形態も同様にして重複した部分の説明を省略する。第2の実施形態は、第1の実施形態の変形例であり、流量制限値及び使用枝数制限値を調節可能としたものであって、図1に示した構成に比べ、制限値設定部13、スケジュールパターン表示部14及び制限値表示部15を付加した構成となっている。ここで、制限値設定部13は、スケジュールパターン記憶及び出力部12によるスケジュールパターンの出力後、ユーザの操作に応じて、流量制限値及び使用枝数制限値を個別に設定する制限値設定機能と、当該設定された流量制限値及び使用枝数制限値に基づいて、制限ネットワーク作成部10、スケジュールパターン生成部11、及びスケジュールパターン記憶及び出力部12を再実行するように制御する再実行制御機能とをもっている。また、制限値設定部13の制限値設定機能は、流量制限値のうちの流量下限値をFLOWLと表すとき、下記式4に示すように、流量下限値FLOWLを流量制限値として算出し、当該算出した流量下限値FLOWLを設定する流量下限値設定機能を有していてもよい。なお、下記式4は、流量下限値FLOWLを算出する際に、流量上限値FLOWに基づき、FLOW-1>0のときにFLOW-1とし、FLOW-1≦0のときに0として求める手順をも包含している。同様に、制限値設定部13の制限値設定機能は、使用枝数制限値のうちの使用枝数下限値をUSELと表すとき、下記式5に示すように、使用枝数下限値USELを使用枝数制限値として算出し、当該算出した使用枝数下限値USELを設定する使用枝数下限値設定機能を有していてもよい。なお、下記式5は、使用枝数下限値USELを算出する際に、使用枝数上限値USEに基づき、USE-1>0のときにUSE-1とし、USE-1≦0のときに0として求める手順をも包含している。スケジュールパターン表示部14は、現在のスケジュールパターンを表示するものであり、具体的には、スケジュールパターン記憶及び出力部12により出力されたスケジュールパターンを表示する機能をもっている。制限値表示部15は、現在の流量制限値及び使用枚数制限値を表示するものであり、具体的には、流量制限値算出部7に算出又は設定された流量制限値と、使用枝数制限値算出部9に算出又は設定された使用枝数制限値とを表示する機能をもっている。次に、以上のように構成された繰り返し型スケジュール生成装置の動作を説明する。いま、第1の実施形態に述べた動作が終了し、スケジュールパターン記憶及び出力部12からスケジュールパターンが出力されたとする。このとき、制限値表示部15は、図27に一例を示すように、現在の制限値として、流量上限値FLOW、流量下限値FLOWL、使用枝数上限値USE、使用枝数下限値USELをプロトコールj毎に表示する。また、スケジュールパターン表示部14は、図28に一例を示すように、現在のスケジュールパターンとして、当該出力されたスケジュールパターンを表示する。ここで、ユーザが規則性の制約をもっと緩くして平準化を進めたいと考えた場合、制限値設定部13は、図29に一例を示すように、ユーザの操作に応じて、一部の制約を、例えばプロトコールj前16の使用枝数上限値USEを3から4に緩めて設定する。しかる後、制限値設定部13は、この使用枝数上限値USE4を使用枝数制限値算出部9を介して制限ネットワーク作成部10に送出することにより、制限ネットワーク作成部10、スケジュールパターン生成部11、及びスケジュールパターン記憶及び出力部12を再実行するように制御する。再実行の結果、繰り返し型スケジュール生成装置は、使用枝数上限値USEを緩めた設定によってスケジュールパターンのバリエーションが増えるため、もし図28に示したスケジュールパターンよりも平準化が良い解、即ち、最大作業時間が短い解が存在すれば、図30に一例を示すように最大作業時間が4時間の解をもつスケジュールパターンを生成する。逆に、制約を厳しくして規則性をより増したい場合、制限値設定部13は、ユーザの操作に応じて、流量上限値FLOWに加えて流量下限値FLOWLを流量制限値算出部7に設定すると共に、使用枝数上限値USEに加えて使用枝数下限値USELを使用枝数制限値算出部9に設定する。ユーザは、これら下限値FLOWL,USELの制約を規則性を増したい特定のプロトコールjにのみ適用する。例えば図31に示すように、プロトコールj肝8の各下限値FLOWL,USELの制約を適用すると、図32に示すように、規則性を増したスケジュールパターンが得られる。なお、制限値設定部13による設定は、ユーザが手動で数字を入力して行っても良いし、図29及び図31に示すように、値の増加ボタン及び減少ボタン等の入力画面を設けて、値を増加又は減少させたときに当該値を自動的に設定しても良い。また、制限値設定部13は、ユーザの操作に応じて、増減された流量上限値FLOWに対応して式4等から算出した流量下限値FLOWLを設定してもよく、同様に、ユーザの操作に応じて、増減された使用枝数上限値USEに対応して式5等から算出した使用枝数下限値USELを設定してもよい。また同様に、制限値設定部13は、これら設定された流量下限値FLOWL及び使用枝数下限値USELをそれぞれ流量制限値算出部7及び使用枝数制限値算出部9を介して制限ネットワーク作成部10に送出することにより、制限ネットワーク作成部10、スケジュールパターン生成部11、及びスケジュールパターン記憶及び出力部12を再実行するように制御する。これにより、流量上限値FLOW及び使用枝数上限値USEの増減に対応して流量下限値FLOWL及び使用枝数下限値USELの制約を含めたスケジュールパターンを出力することができる。上述したように本実施形態によれば、スケジュールパターンの出力後、ユーザの操作に応じて、流量制限値及び使用枝数制限値を個別に設定し、当該設定された流量制限値及び使用枝数制限値に基づいて、制限ネットワーク作成部10、スケジュールパターン生成部11、及びスケジュールパターン記憶及び出力部12を再実行するように制御する構成により、第1の実施形態の効果に加え、流量制限値及び使用枝数制限値による制約を調整して、スケジュールパターンにおける規則性と偏り具合を調整することができる。補足すると、本実施形態は、割当時間種類ノード月,...,金及び割当時間種類アークβを生成して、各割当時間種類アークβに使用枝数制限をつけてスケジュールパターンが取り得るパターン種類の上限値、または上限値と下限値を算出して制約をかけることで、どの程度週単位の規則性が保たれるかが保証でき、各日の患者数をなるべく一定にする又は偏りをへらすという目的関数については、制約を満たす範囲で最適化できる。また、ユーザが複数ある目的関数のウェイトを調節するよりも容易に操作ができるという効果がある。さらに、流量制限値及び使用枝数制限値の設定は、式1、式2、式4及び式5に示したデフォルト設定の他にユーザが個別に調節することができるため、少し制約をゆるくしたい場合などにも対応できる。<第3の実施形態>次に、第3の実施形態に係る繰り返し型スケジュール生成装置について説明する。第3の実施形態は、第1又は第2の実施形態の変形例であり、前述した繰り返し型スケジュール生成装置を、放射線治療に代えて、定期作業のスケジュール立案で、作業員に分かりやすいスケジュールを立案する場合に適用した形態となっている。これに伴い、前述した繰り返し型スケジュール生成装置は、以下の点1~14が変更されている。1放射線治療に代えて、定期作業を用いる。2各プロトコールに代えて、各作業種類を用いる。3各治療工程に代えて、各作業工程を用いる。4患者数に代えて、作業対象数を用いる。5治療部位と照射治療工程の回数に代えて、点検作業又は工具の交換作業を用いる。6各治療工程間に空ける日数及び前記照射治療工程の回数に関する制約に代えて、作業工程間に空ける日数に関する制約を用いる。7曜日と週との組合せに代えて、作業点を用いる。8割当時間情報に代えて、作業点情報を用いる。9曜日種類に代えて、作業員を用いる。10割当時間種類情報に代えて、作業点種類情報を用いる。11割当時間ノードに代えて、作業点ノードを用いる。12プロトコールノードに代えて、作業種類ノードを用いる。13割当時間種類ノードに代えて、作業点種類ノードを用いる。14割当時間種類アークに代えて、作業点種類アークを用いる。これらの点1~14の変更は、以下のような構成を表している。スケジュールパターンは、図33に示すように、繰り返し期間中の作業点割当時間毎に、定期作業の各作業種類に対して基点となる作業工程割当作業の日を示すように当該各作業種類が割当てられている。スケジュールパターンは、各作業種類を当該各作業種類の作業対象数だけ含んでいる。なお、図33では、分かりやすさのため、スケジュールパターンの表示の上部に作業員の種類を記載したが、この記載は省略してもよい。割当作業情報記憶部第1の記憶部2は、各作業工程における点検作業又は工具の交換作業を示す各作業種類と、各作業種類の作業対象数とを規定した割当作業情報を記憶する。作業点情報記憶部第4の記憶部3は、各作業種類を割当て可能な作業点を個別に示す作業点情報を記憶する。作業点種類情報記憶部第5の記憶手段4は、作業点情報内の作業点に個別に割り当てられた作業員を規定した作業点種類情報を記憶する。作業工程情報は、各作業工程を規定した情報であり、例えば、各作業種類に対する各作業工程及び当該各作業工程のうちの基点となる作業工程以下、基準作業工程ともいうを規定した情報である。なお、基準作業工程では、各日の基準作業工程を実施する作業時間の上限が指定可能となっている。他の作業工程は、基準作業工程からの相対位置間に空ける日数でスケジュールが生成される。制約情報は、作業工程の制約を規定した情報であり、具体的には、各作業工程間に空ける日数に関する制約を規定した情報である。割当ネットワーク作成部6は、各記憶部2,3内の作業点情報及び割当作業情報に基づいて、割当ネットワークを作成する機能をもっている。割当ネットワークは、図34に示すように、当該作業点情報に示される作業点を個別に示す複数の作業点ノード1,...,12と、当該割当作業情報に規定された各作業種類を個別に示す複数の作業種類ノードA,...,Dと、各作業点ノードと各作業種類ノードとの組合せを個別に示す複数の割当アークαと、各作業種類の作業対象数6,...,1とを含んでいる。流量制限値算出部7は、各記憶部3,5内の割当作業情報及び作業点情報に基づいて、各割当アークα毎に、割当て可能な作業対象数としての流量を制限する流量制限値を算出する流量制限値算出機能をもっている。ここで、流量制限値算出機能は、割当作業情報に規定された各作業種類をjと表し、繰り返し期間中の日の各作業種類jの作業対象数をdjと表し、作業点情報に示される各作業点の集合をIと表し、この集合Iに含まれる各作業点の個数を|I|と表し、流量制限値のうちの流量上限値をFLOWと表すとき、前述した式1に示すように、流量上限値FLOWを流量制限値として算出する流量上限値算出機能を備えてもよい。種類要素生成部8は、各記憶部3,4内の作業点情報及び作業点種類情報に基づいて、作業員を個別に示す複数の作業点種類ノードと、各作業点種類ノードのうちの有効な作業点種類ノードを個別に示す複数の作業点種類アークとを生成する機能をもっている。なお、種類要素の語は作業点種類ノード及び作業点種類アークと読み替えてもよい。使用枝数制限値算出部9は、割当作業情報、作業点情報及び作業点種類情報に基づいて、各作業種類毎に、各作業点種類アークの使用枝数を制限する使用枝数制限値を算出する使用枝数制限値算出機能をもっている。ここで、使用枝数制限値算出機能は、作業点情報に示される作業点の個数|I|を作業点種類情報に示される作業員の個数|W|で除した値をh=|I|/|W|と表し、使用枝数制限値のうちの使用枝数上限値をUSEと表すとき、前述した式2に示すように、使用枝数上限値USEを使用枝数制限値として算出する使用枝数上限値算出機能を備えてもよい。制限ネットワーク作成部10は、割当ネットワーク、流量制限値、各作業点種類ノード、各作業点種類アーク及び使用枝数制限値を統合することにより、制限ネットワークを作成する機能をもっている。制限ネットワークは、図35に示すように、各作業点ノード1,...,12と、各作業種類ノードA,...,Dと、各割当アークαと、作業対象数6,...,1と、流量制限値図示せずと、各作業点種類ノード図中、作業員ノードa,...,cと、各作業点種類アークβと、使用枝数制限値図示せずとを含んでいる。なお、制限ネットワークの語は、作業点種類をもつネットワークと読み替えてもよい。スケジュールパターン生成部11は、制限ネットワーク、作業工程情報及び制約情報に基づいて、図33に示した如き、スケジュールパターンを生成するスケジュールパターン生成機能と、当該生成したスケジュールパターンをスケジュールパターン記憶及び出力部12に書込むスケジュールパターン書込機能とをもっている。スケジュールパターン記憶及び出力部12は、前述した機能をもっている。次に、以上のように構成された繰り返し型スケジュール生成装置の動作を説明する。例えば、定期作業の回数によって、定期的に点検や工具の取り替えが必要な作業が複数種類あり、各定期作業は作業種類に応じて取り替えのタイミングが異なる。作業点とは、定期作業毎の作業間隔の公約数である単位作業間隔で設定されている、定期作業が入る可能性のあるタイミングのことである。作業種類を作業点に割り当てることでスケジュールパターンが生成可能となっている。第2の先行技術では、定期作業が必ず一定間隔で割当てられるスケジューリングを行っている。一方、本実施形態では、定期作業を作業点の繰り返し期間ごとで繰返すスケジュールとして、繰返し期間内は、必ずしも定期作業が一定間隔でなくても良いとし、作業員が作業を行う作業点が作業員ごとに得られている場合、作業員に分かりやすいスケジュールパターンを生成できる。このようなスケジューリングを行うにあたって、情報入力部1は、ユーザの操作に応じて割当作業情報、作業点情報及び作業点種類情報をそれぞれ割当作業情報記憶部2、作業点情報記憶部3及び作業点種類情報記憶部4に書込む。同様に、情報入力部1は、ユーザの操作に応じて、作業工程情報及び制約情報をスケジュール情報記憶部5に書込む。続いて、割当ネットワーク作成部6は、各記憶部2,3内の作業点情報及び割当作業情報に基づいて、例えば、繰り返し期間における作業点が12個、作業種類が4種類の場合、図34に示す如き、当該作業点情報に示される作業点を個別に示す複数の作業点ノード1,...,12と、当該割当作業情報に規定された各作業種類を個別に示す複数の作業種類ノードA,...,Dと、各作業点ノードと各作業種類ノードとの組合せを個別に示す複数の割当アークαと、各作業種類の作業対象数6,...,1とを含む割当ネットワークを作成する。各作業種類を各作業点に割り当てたときの多次元コストは、立案期間中に各作業が繰り返し期間における作業点数毎に繰り返すとすれば、図23と同様に算出できるため、各日の作業時間の平準化を可能としている。ここで、前述同様に、流量制限値算出部7が流量制限値を算出する。しかる後、種類要素生成部8は、各記憶部3,4内の作業点情報及び作業点種類情報に基づいて、図35に示す如き制限ネットワークに統合される、3つの作業点種類ノード図中、作業員a,...,cを生成すると共に、各作業点種類ノードとソースノード、各作業点種類ノードと作業種類ノード間にアークを生成する。このように、3つの作業点種類ノードa,...,cが生成された場合、使用枝数制限値算出部9は、図36に示すように、各作業点種類アークβの使用枝数上限値USEを算出することができる。例えば、作業種類Aの場合、dj=6、h=|I|/|W|=12/3=4、として式2に基づき、USEは、dj/h=6/4以上の最小の整数2として算出される。なお、図36の場合、作業員aは作業点3つ、作業員bは5つ、作業員cは4つと受け持つ作業点の数が不均等になっているので、繰り返し回数としてのh=4が必ずしも好ましいわけではない。ここで述べたhの計算の仕方は一例であり、これに限定されない。続いて、制限ネットワーク作成部10は、割当ネットワーク、流量制限値、各作業点種類ノードa,...,c、各作業点種類アークβ及び使用枝数制限値を統合することにより、図35に示した如き、制限ネットワークを作成する。スケジュールパターン生成部11は、前述同様に定式化可能な制限ネットワークと、スケジュール記憶部5内の作業工程情報及び制約情報に基づいて、前述した式3の評価値zを最小化するように混合整数線形計画問題を解くことにより、図33に示したように、スケジュールパターンを生成する。また、作業員の数及び担当作業点が図37に示すように変更された場合にも、各作業点種類アークβの使用枝数上限値USEを図38に示すように算出でき、図39に示す如き、スケジュールパターンを生成することができる。上述したように本実施形態によれば、第1の実施形態に比べ、曜日と週との組合せに代えて作業点を用い、割当時間情報に代えて作業点情報を用い、曜日種類に代えて作業員を用い、割当時間種類情報に代えて作業点種類情報を用いる構成により、第1の実施形態の効果に加え、作業点で作業を行う作業員を作業点種類として用いることで、その作業点に特定の作業種類が集まりやすいスケジュールパターンを作成できる。また、3つの作業点種類ノードa,...,cに対応する3つの作業点種類アークβを生成し、これら作業点種類アークβの使用枝数制限をつける構成により、作業点における作業量を平準化してかつ、作業点種類に応じて特定の作業種類が集まりやすいスケジュールパターンを容易に作成できるため、より作業者にわかりやすい計画を作成することができる。また、同じ作業員が担当する作業点を自由に与えた場合にも同じ作業員になるべく同じ作業が当たるスケジュールを作成することができる。なお、上記の各実施形態に記載した手法は、コンピュータによって実行させることができ、また、コンピュータによって実行させることのできるプログラムとして、磁気ディスクフレキシブルディスク、ハードディスクなど、光ディスクCD-ROM、DVDなど、光磁気ディスクMO、半導体メモリなどの記憶媒体に格納して頒布することもできる。また、各実施形態に記載した装置は、ハードウェア構成、又はハードウェア資源とソフトウェアとの組合せ構成のいずれでも実施可能となっている。組合せ構成のソフトウェアとしては、予めネットワーク又は記憶媒体からコンピュータにインストールされ、各実施形態に記載した装置の機能を当該コンピュータに実現させるためのプログラムが用いられる。すなわち、各実施形態に記載した装置は、コンピュータのメモリ内の当該プログラムを当該コンピュータのプロセッサが実行することにより、実施してもよい。また、本発明のいくつかの実施形態を説明したが、これらの実施形態は、例として提示したものであり、発明の範囲を限定することは意図していない。これら新規な実施形態は、その他の様々な形態で実施されることが可能であり、発明の要旨を逸脱しない範囲で、種々の省略、置き換え、変更を行うことができる。これら実施形態やその変形は、発明の範囲や要旨に含まれるとともに、特許請求の範囲に記載された発明とその均等の範囲に含まれる。"