first(m) = {res = List(); 
a333538 = List(); a332558 = List(); 
my(n = 1, nk1 = 2, p = 1, r = 0); 
	for(i = 1, m, 
		if(i % 10^6 == 0, 
			print1(i \ 10^6", ");
		);
		while(p % nk1 != 0, 
			p*=nk1; 
			nk1++; 
		); 
		c = hpf(nk1); 
		if(c > r, 
			r = c; 
			listput(a333538, i); 
			\\print(); print(); 
			\\ print(a333538); \\ A333538
			\\print(); print(); 
			
		); 
		\\ print1(hpf(nk1)", "); \\ A333537 
		p \= i; 
		listput(a332558, nk1 - i - 1); 
		print1(nk1 - i - 1", "); \\ A332558
		\\print1(nk1", "); \\ A332559
		
	); a332558
}

hpf(n)= my(f = factor(n)); f[#f~, 1]