fork download
  1. p=[1,2,3,8]
  2. q=[0,1,1,3]
  3. n=int(input());
  4. for i in range(0,4100) :
  5. a=1;
  6. if i%3==2 :
  7. a=2*(i//3 + 2)
  8. p.append(a*p[-1]+p[-2])
  9. q.append(a*q[-1]+q[-2])
  10. def sa(b,i,m) :
  11. res=m*b*p[i]
  12. if i&1 :
  13. res+=(b*p[i]*q[i]-q[i]+p[i]-1)*b // 2
  14. else :
  15. res+=(b*p[i]*q[i]-q[i]+p[i]+1)*b // 2
  16. return res
  17. ans=0
  18. for i in range(4000,0,-1) :
  19. if q[i]<=n and q[i+1]>n :
  20. if q[i]==n :
  21. ans+=sa(1,i,0)
  22. n=0;
  23. else :
  24. b=n//q[i]
  25. m=n%q[i]
  26. ans+=sa(b,i,m)
  27. n=m
  28. print(ans)
Success #stdin #stdout 0.02s 27704KB
stdin
3
stdout
14