{
	/* produces A072871 in disorder */
	lim  = 10 000;
	rem  = lim; /* number of missing terms */
	miss = vector(lim); /* index of missing terms */
	for(i=1,lim,miss[i]=i);

	n=0;
	while (rem>0,
		n=n+1;
		p=numbpart(n);
		forstep(m=rem,1,-1,
			if(p % miss[m]==0,
				print(miss[m] " " p);
				miss[m]=miss[rem];
				rem=rem-1;
			)
		)
	)
}