// x1=16598013x+12820163 mod 2^24
#include <stdio.h>
#include <math.h>
void main()
{
long y, x, xi, m ;
x=16598013 ;
m=16777216 ; // 2^24
for (xi=1; xi<m; xi++)
{
y=(x*xi)%m ;
if (y==1 || y==-1)
{
printf ("Found y=%d : xi=%d, x*xi=%ld\n", y, xi, x*xi) ;
}
}
}