twitterで見つけた誰かの解法です.誰だったか忘れました...><
lcm(x,y,z)=k となる (x,y,z) の個数を C(k) で表すと,
Sn=k=1∑nC(k)×⌊kn⌋
である.
ここで,k=p1e1p2e2⋯pmem とすると
C(k)=((e1+1)3−e13)×((e2+1)3−e23)×⋯×((em+1)3−em3)
なので,常に奇数となる.従って
Sn≡k=1∑n⌊kn⌋(mod2)
であり,あとは他の方のユーザー解説と同様.
いろんな解法があって面白いなぁと思いました.