Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- import sys
- from decimal import *
- setcontext(Context(prec=5000,rounding=ROUND_DOWN))
- factorial=[1]
- fac=Decimal('1')
- for i in range(2,225):
- fac*=Decimal(i)
- factorial.append(fac)
- while True:
- s=sys.stdin.readline()
- if not s:
- break
- s=Decimal(s)
- sin=s
- for i in range(5,225,4):
- sin+=(Decimal((s**Decimal(i)))/Decimal(factorial[i-1]))
- for i in range(3,225,4):
- sin-=(Decimal((s**Decimal(i)))/Decimal(factorial[i-1]))
- print("{:.50f}".format(sin))
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement