As in the previous solution, label the irreducible lattice points (x1,y1),…,(xn,yn) and assume without loss of generality that no two of the points are collinear with the origin. We induct on n to construct a homogenous polynomial f(x,y) such that f(xi,yi)=1 for all 1⩽i⩽n.
If n=1 : Since x1 and y1 are relatively prime, there exist some integers c,d such that cx1+dy1=1. Then f(x,y)=cx+dy is suitable.
If n⩾2 : By the induction hypothesis we already have a homogeneous polynomial g(x,y) with g(x1,y1)=…=g(xn−1,yn−1)=1. Let j=degg,
gn(x,y)=k=1∏n−1(ykx−xky)
and an=gn(xn,yn). By assumption, an=0. Take some integers c,d such that cxn+dyn=1. We will construct f(x,y) in the form
f(x,y)=g(x,y)K−C⋅gn(x,y)⋅(cx+dy)L
where K and L are some positive integers and C is some integer. We assume that L=Kj−n+1 so that f is homogenous.
Due to g(x1,y1)=…=g(xn−1,yn−1)=1 and gn(x1,y1)=…=gn(xn−1,yn−1)=0, the property f(x1,y1)=…=f(xn−1,yn−1)=1 is automatically satisfied with any choice of K,L, and C.
Furthermore,
f(xn,yn)=g(xn,yn)K−C⋅gn(xn,yn)⋅(cxn+dyn)L=g(xn,yn)K−Can.
If we have an exponent K such that g(xn,yn)K≡1(modan), then we may choose C such that f(xn,yn)=1. We now choose such a K.
Consider an arbitrary prime divisor p of an. By
p∣an=gn(xn,yn)=k=1∏n−1(ykxn−xkyn)
there is some 1⩽k<n such that xkyn≡xnyk(modp). We first show that xkxn or ykyn is relatively prime with p. This is trivial in the case xkyn≡xnyk≡0(modp). In the other case, we have xkyn≡xnyk≡0(modp). If, say p∣xk, then p∤yk because (xk,yk) is irreducible, so p∣xn; then p∤yn because (xn,yn) is irreducible. In summary, p∣xk implies p∤ykyn. Similarly, p∣yn implies p∤xkxn.
By the homogeneity of g we have the congruences
xkd⋅g(xn,yn)=g(xkxn,xkyn)≡g(xkxn,ykxn)=xnd⋅g(xk,yk)=xnd(modp)
and
ykd⋅g(xn,yn)=g(ykxn,ykyn)≡g(xkyn,ykyn)=ynd⋅g(xk,yk)=ynd(modp).
If p∤xkxn, then take the (p−1)st power of the first congruence; otherwise take the (p−1)st power of the second; by Fermat's theorem, in both cases we get
g(xn,yn)p−1≡1(modp)
If pα∣m, then we have
g(xn,yn)pα−1(p−1)≡1(modpα)
which implies that the exponent K=n⋅φ(an), which is a multiple of all pα−1(p−1), is a suitable choice. (The factor n is added only so that K⩾n and so L>0.)