ll m=100003; ll t[64][m]; rep(i,1,m){ rep(j,i,m,i){ t[0][j]+=i; } } rep(i,m){ t[0][i]%=m; } rep(l,1,64){ rep(i,m){ t[l][i]=t[l-1][t[l-1][i]]; } } ll@n,@k--; if(k){ n=DivisorSum(n)%m; --k; rrep(l,64){ if(k&1l<