不同的点

【题目描述】
陶陶为了给一道平面几何题出数据，需要产生 N 个点(x[i],y[i])。已知 x,y 是由伪随机函数顺序产生，即：
X[i+1] = (X[i]*Ax+Bx+i) mod Cx (X[1], Ax,Bx,Cx 是事先给定的）
Y[i+1] = (Y[i]*Ay+By+i) mod Cy (Y[1], Ay,By,Cy 是事先给定的）
这样，就可以快速连续产生很多点坐标(X[i], Y[i])。
不幸的是，这样产生的点有可能有相同的，虽然这种几率很少，但严谨的陶陶不允许这种事发生。陶陶要求你帮助他解决最少要产生前多少项时，正好有 N 个不相同的点。
【输入格式】
第一行:一个整数 N
第二行：4 个整数 X[1]、 Ax、Bx、Cx
第三行：4 个整数 Y[1]、 Ay、By、Cy
【输出格式】
一个整数 M表示最少要连续产生 M 个点，正好有 N 个不相同的点。数据保证有答案。
【样例输入】
21
2 4 3 6
5 2 3 13
【样例输出】
24
【数据范围】
1<=N<=1,000,000, 其它所有数据都在[1…1,000,000,000]范围内。
【分析】
如果直接模拟会超时，故使用哈希优化。

const
  inf=1123357;
var
  i,n,xi,ax,bx,cx,yi,ay,by,cy,m:qword;
    xx,yy,h:array[0..inf]of qword;
procedure insert(x,y:qword);
var
  t,i:qword;
begin
  t:=(x*319+y*131) mod inf;
    while h[t]<>0 do begin
      i:=h[t];
        if (xx[i]=x)and(yy[i]=y) then exit;
        t:=(t+1) mod inf;
    end;
    inc(m);
    xx[m]:=x;yy[m]:=y;
    h[t]:=m;
end;
begin
  readln(n);
    readln(xi,ax,bx,cx);
    readln(yi,ay,by,cy);
    fillchar(h,sizeof(h),0);
    m:=0;
    i:=1;
    insert(xi,yi);
    while m<n do begin
      xi:=(xi*ax+bx+i) mod cx;
      yi:=(yi*ay+by+i) mod cy;
        insert(xi,yi);
      inc(i);
    end;
    write(i);
end.