解题报告-NOI2002 Day2 savage【信息技术】.docVIP

NOI2002 Day2 Problem1《荒岛野人》解题报告 天津市南开中学 罗勐 二○○二年十一月 问题分析： 看过试题容易想到，在所以人中选出最大起始位置max，从max起逐个检测是否满足条件。但由于洞的个数最大为106，所以单纯摸拟肯定要超时（只能过4个点），必须进行优化。进一步分析，这是一个追击问题，对于某两个野人，如果他们速度相同则永远满足条件，若速度不同则要看在两人有生之年，速度快的野人能否追到速度慢的野人或超过其整数圈。 设两野人为w1,w2;rab[w1,w2].p为速度差的绝对值;rab[w1,w2].c为速度慢的野人起始位置－速度快的起始位置；rab[w1,w2].L为较短的寿命值，则洞的个数m需要满足的条件为：方程 x[w1,w2]*rab[w1,w2].p≡rab[w1,w2].c(mod m) (1 =w1,w2 =n, x[w1,w2] =rab[w1,w2].L)永不成立. 以下简写为xp≡c(mod m) 队友肖天的做法： 在不了解同余方程的情况下，是一种好方法（可过9个点）。 将同余方程xp≡c(mod m)改写为不定方程xp=ym+c，由于1 =x =106而1 =p =100,所以 (ym+c) mod p最多只有100种可能的取值，从1－100依次搜索y若值为0且x L,则m不符合要求，若值出现重复，则此二人永远不会碰到，继续搜索直至m对所有方程均成立。 标答：解同余方程。（参见黑书P63） 对于xp≡c(mod m) 1．应用扩展欧几里德法（辗转相除）（参见后面的辅助说明） 求出x’p+y’m=gcd(p,m)(最大公约数) 设d=gcd(p,m)即x’p=d+y’m ( x’p≡d(mod m) 2. 若c mod d 0则方程无解（参见附录定理 4），二人永远不会相遇。 若c mod d=0 则方程有d个解（参见附录定理 6）， 这些解为x0=x’*(c div d) mod m=x’*(c/d) mod m，（参见附录定理 5） xi=(x0+i*(m/d)) mod m=(x’*(c/d)+i*(m/d)) mod m (i=0..d-1) 若xi L方程无解，否则m不合条件继续检测更大的m。 优化： 只需求出xmin并判断其与L的关系，根据同余知识可知xmin=x’*(c/d) mod (m/d) 标程解了两次同余方程，是没有必要的，重复了。 另外，此时求xi可简化为xi=xmin+i*(m/d)(i=0..d-1). 辅助说明：以下给出报告中提到一些数学知识。 扩展欧几里德（辗转相除）函数： function extended_euclid(a,b:longint;var x,y:longint):longint; var t:longint; begin if b=0 then begin extended_euclid:=a;x:=1;y:=0 end else begin extended_euclid:=extend_euclid(b,a mod b,x,y); t:=x;x:=y;y:=t-(a div b)*y end; end; 此函数满足：ax+by=extended_euclid=gcd(a,b) 有关一次同余方程的定理归纳（请先阅读同余的一些基本性质）。 对于同余式：ax≡b(mod m)的一个解并不是一个数，而是包含对模m同余的一组数（此题中只求最小的）——模m的一个剩余类。即，凡对模m同余的数，算作一个解，只有对模m互不同余的数，才是不同的解。 剩余类：模m同余的整数的集合。 完全剩余系（完系）：从模m的每一个剩余类中任意挑出一个整数，则这n个整数就称为模m的一个完全剩余系。 定理1 若(a,m)=1,则必存在a*,使aa*≡b(mod m)。 证明：因为(a,m)=1,所以存在a*,y,使aa*+my=1.由此aa*≡1(mod m), a*称为a的数论倒数。只有与模互素的数才有数论倒数，并且一个数的数论倒数在同余意义下是唯一的。 定理2 若a≡b(mod n),m|n,则a≡b(mod m)。 证明：因为a≡b(mod n),所以n|a-b,又因为m|n,所以m|a-b,故a≡b(mod m)。 定理3 若(a,m)=1,则同余式ax≡b(mod m)有且只有一解. 证明：设1,2,3…是模m的一个完全剩余系，(a,m)=1,所以a,2a,…,ma也是模m的一个完全剩余系（可用反正法证明），其中有一数，设为ak，满足ak≡b(mod m),则x≡k(mod m)就是ax≡b(mod m)的唯一解。 定理4 设(a,m)=d 1,同余式ax≡b(mod m)有解的通分必要条件是d|b. 证明：必