数论复习-补充【信息技术】.docVIP

1、数据结构 设 const maxN=200;maxP=100;{房间数和转盘格子数的上限} LimitAns=1000000000;｛庆典最早日期的初始值｝ type Tarr = array[1..maxN] of integer;{置换的数据类型} var ans:longint;{ 第一次举行庆典的日期} data:array[1..maxP] of Tarr;{data[i,j]为寺院建立时房间j转盘的第i个格子（顺时针）的数字，由输入得到} NowP,NextC:Tarr;{ NowP为转盘当前格的置换ci；NextC为轮换性的置换Ｃ} n, p:integer;{ 房间数和转盘格子数} aa,nn:Tarr；{当前置换中祭司i到位的天数为aa[i]，周期为nn[i]} ２、设置初始时各祭司的房间号 procedure initArr(var a : Tarr); var i : integer; begin fillchar(a,sizeof(a),0); for i←1 to n do a[i]←i; end;{ initArr } 3、计算一次置换 我们输入当前祭司的房间序列a和各房间转盘的数字序列b，通过调用过程NextPosition计算置换c： procedure NextPosition(const a,b:Tarr;var c:Tarr); var i : integer; begin fillchar(c, sizeof(c), 0); for i←1 to n do c[i]←b[a[i]];{祭司i根据所在房间a[i]的转盘数字b[a[i]]走入c[i]房间} end;{NextPosition} ４、建立同余方程aa和nn（w=aa[i] mod bb[i](1≤i≤n)） 　　我们输入第FA年的置换Posi。如果n个祭司通过若干次轮换性置换NextC回到自己编号的房间，则返回true，并建立对应的同余方程aa和nn；若有祭司无法回到自己编号的房间，则返回false： function FillFC(const Posi :Tarr; FA:longint;var aa,nn:Tarr):boolean; var i,j,k:integer; ok: array[1..maxN] of boolean; begin FillFC←False; for i←1 to n do{搜索每一个祭司} begin fillchar(ok, sizeof(ok), true);{所有房间未访问} j←Posi[i];{取出祭司i所在的房间} k←0;{计算祭司i经过轮换性置换的次数k} while(ji)and ok[j]do begin ok[j]←false;inc(k);j←NextC[j];end;{while} 　if　ji　then exit；｛若祭司i不能回到自己编号的房间，则失败退出｝ Aa[i]←k*p+FA｛设置祭司i到位的年数Aa[i]｝ end;{for} for i←1 to n do{计算每个祭司的周期} begin fillchar(ok, sizeof(ok), true); j←i;k←0; {计算祭司i经过轮换性置换的次数k} repeat inc(k);j←NextC[j]; if not ok[j] then exit; ok[j]←false; until j = i; Nn[i]←k*p{计算祭司i的周期Nn[i]} end;{for} for i←1 to n-1 do{若出现祭司对的周期非互质的情况，则返回失败信息} for j←i to n do if gcd(nn[i],nn[j])1 then exit; FillFC←True;{返回成功信息} end;{ FillFC } 主算法 procedure Main; var i, j, k,s : longint; temp, temp2 : Tarr; f:boolean； begin initArr(NowP);{初始时祭司走入与自己编号相同的房间} for i←1 to p do{转盘逆时针转p格} begin NextPosition(NowP,data[i],temp); NowP←temp;{计算第i年的置换NowP} initArr(temp); k←i;{计算轮换性的置换NextC } fo