a=Mod[1504170715041707*1,4503599627370517]

b=Mod[1504170715041707*3,4503599627370517]

a+b

d=4503599627370517

Grid[Table[{i,BaseForm[d,i]},{i,2,32}]]

m=1504170715041707

Grid[Table[{i,BaseForm[m,i]},{i,2,32}]]

FindInstance[Mod[1504170715041707*n,4503599627370517]==10487287,n]

BaseForm[ModularInverse[m,d],16]

inv=ModularInverse[m,d]

36198519674524953353568

Mod[1504170715041707/3,4503599627370517]

Mod[1504170715041707*inv,d]

Mod[10487287*inv,d]

Mod[1504170715041707*3732049906,4503599627370517]

Mod[10487281*inv,d]

Mod[1504170715041707*1808058673188507,4503599627370517]

Mod[107159*inv,d]

Total[Import["~/tmp/pe700.txt","List"]]

t=Table[Mod[i*inv,d],{i,1,668000}]

Sort[t]

Mod[m*10543898410,d]

Mod[m*10827725431,d]

Mod[m*21655450862,d]

first=1504170715041707

Log[10,first]//N