n = int(input())
mod = 10**9+7

ans = [1,12,77,237]
if n <= 3:
    print(ans[n])
    exit()

le = 1+6*n
sm = le-6
edge = n*(n+1)//2
edge %= mod
edge2 = edge-n
edge2 %= mod

ans = sm**2-edge2*4
ans %= mod

base = (le**2-edge*4)-(sm**2-edge2*4)
base %= mod

ans += 36*n-12
ans %= mod
print(ans)