program-id.y.data division.working-storage section. 1 i pic 99. 1 k pic 9(8). 1 n pic 9(11). 1 s pic x(8). procedure division.a.accept s move s to k move 0 to n i.b.add 1 to i compute k=k/i compute n=n+function mod(k,1+i)*10^i if k>0 go b.display n(11- function log10(n):)if 13=s(6:)display 1000 else if s<>120 if s<>508 go a.
Note that non-ascii characters in the above source code will be escaped (such as \x9f).