1divisorslenstra(245784, 1, 65)
2divisorslenstra(245784, 19, 65)
3divisorslenstra(288288, 1, 71)
4divisorslenstra(288288, 28, 71)
5divisorslenstra(2841696, 2, 175)
6divisorslenstra(2841696, 23, 175)
7divisorslenstra(1680,1,13)
8divisorslenstra(4,1,2) \\ #1930
9divisorslenstra(2,1,2)
10
11divisors(1)
12divisors(12)
13divisors(1, 1)
14divisors(12, 1)
15{
16for (i=1,100,
17  d = divisors(i,1);
18  for(j=1,#d, if (factorback(d[j][2])!=d[j][1],error(i))));
19}
20
21fordivfactored(100, d, print(d))
22