a = vector(10 000)

{
	my (n=0);
	for prime (p=1, oo,
		n++;

		print (n, " ", a[n] = p - sum(k=1, n-1, if (gcd(k, n+1)==1, a[k], 0)));
		if (n==#a,
			break;
		);
	);
}

quit