Арсен Жуматай,304 ArbTangent 6429			Жангирхан Шаку, 404 ArbTangent 6418
f	1	from decimal import Decimal, getcontext, ROUND_HALF_UP	f	1	from decimal import Decimal, getcontext, ROUND_HALF_UP
	2	import sys		2	import sys
	3			3
n	4	def compute_pi(prec):	n	4	def compute_pi(precision):
	5	getcontext().prec = prec + 5		5	getcontext().prec = precision + 5
	6	C = 426880 * Decimal(10005).sqrt()		6	coeff = 426880 * Decimal(10005).sqrt()
	7	M = Decimal(1)		7	M = Decimal(1)
	8	L = Decimal(13591409)		8	L = Decimal(13591409)
	9	X = Decimal(1)		9	X = Decimal(1)
	10	K = Decimal(6)		10	K = Decimal(6)
n	11	S = L	n	11	total_sum = L
	12	for k in range(1, prec):		12	for k in range(1, precision):
	13	M = M * (K ** 3 - 16 * K) / k ** 3		13	M = M * (K ** 3 - 16 * K) / k ** 3
	14	L += 545140134		14	L += 545140134
	15	X *= -262537412640768000		15	X *= -262537412640768000
	16	term = M * L / X		16	term = M * L / X
n	17	S += term	n	17	total_sum += term
	18	if term == 0:		18	if term == 0:
	19	break		19	break
	20	K += 12		20	K += 12
n	21	pi = C / S	n	21	return coeff / total_sum
	22	return pi
	23			22
n	24	def sin_series(x, terms):	n	23	def series_sine(x, num_terms):
	25	getcontext().prec += 2		24	getcontext().prec += 2
	26	x_power = x		25	x_power = x
	27	factorial = Decimal(1)		26	factorial = Decimal(1)
n	28	sin_x = x	n	27	sine_result = x
	29	sign = -1		28	sign = -1
n	30	for i in range(3, terms * 2, 2):	n	29	for i in range(3, num_terms * 2, 2):
	31	x_power *= x * x		30	x_power *= x * x
	32	factorial *= i * (i - 1)		31	factorial *= i * (i - 1)
	33	term = x_power / factorial		32	term = x_power / factorial
n	34	sin_x += sign * term	n	33	sine_result += sign * term
	35	sign *= -1		34	sign *= -1
	36	if term == 0:		35	if term == 0:
	37	break		36	break
	38	getcontext().prec -= 2		37	getcontext().prec -= 2
n	39	return +sin_x	n	38	return +sine_result
	40			39
n	41	def cos_series(x, terms):	n	40	def series_cosine(x, num_terms):
	42	getcontext().prec += 2		41	getcontext().prec += 2
	43	x_power = Decimal(1)		42	x_power = Decimal(1)
	44	factorial = Decimal(1)		43	factorial = Decimal(1)
n	45	cos_x = Decimal(1)	n	44	cosine_result = Decimal(1)
	46	sign = -1		45	sign = -1
n	47	for i in range(2, terms * 2, 2):	n	46	for i in range(2, num_terms * 2, 2):
	48	x_power *= x * x		47	x_power *= x * x
	49	factorial *= i * (i - 1)		48	factorial *= i * (i - 1)
	50	term = x_power / factorial		49	term = x_power / factorial
n	51	cos_x += sign * term	n	50	cosine_result += sign * term
	52	sign *= -1		51	sign *= -1
	53	if term == 0:		52	if term == 0:
	54	break		53	break
	55	getcontext().prec -= 2		54	getcontext().prec -= 2
n	56	return +cos_x	n	55	return +cosine_result
	57			56
n	58	def round_to_significant_digits(x, digits):	n	57	def significant_digit_round(value, digits):
	59	if x.is_zero():		58	if value.is_zero():
	60	return Decimal(0)		59	return Decimal(0)
	61	else:		60	else:
n	62	sign, digits_tuple, exponent = x.as_tuple()	n	61	sign, digits_tuple, exponent = value.as_tuple()
	63	num_digits = len(digits_tuple)		62	num_digits = len(digits_tuple)
	64	adjusted_exponent = exponent + num_digits - digits		63	adjusted_exponent = exponent + num_digits - digits
	65	quantize_exp = Decimal('1e{}'.format(adjusted_exponent))		64	quantize_exp = Decimal('1e{}'.format(adjusted_exponent))
n	66	return x.quantize(quantize_exp, rounding=ROUND_HALF_UP)	n	65	return value.quantize(quantize_exp, rounding=ROUND_HALF_UP)
	67			66
n	68	def main():	n	67	def perform_calculations():
	69	A = Decimal(sys.stdin.readline().strip())		68	angle_degrees = Decimal(sys.stdin.readline().strip())
	70	E = int(sys.stdin.readline().strip())		69	precision = int(sys.stdin.readline().strip())
	71	if E < 4 or E > 1000:		70	if precision < 4 or precision > 1000:
	72	print('Precision E must be between 4 and 1000')		71	print('Precision must be between 4 and 1000')
	73	return		72	return
n	74	getcontext().prec = E + 10	n	73	getcontext().prec = precision + 10
	75	pi = compute_pi(getcontext().prec)		74	pi_result = compute_pi(getcontext().prec)
	76	angle_rad = A * pi / Decimal(200)		75	angle_radians = angle_degrees * pi_result / Decimal(200)
	77	terms = E + 5		76	terms_required = precision + 5
	78	sin_x = sin_series(angle_rad, terms)		77	sine_value = series_sine(angle_radians, terms_required)
	79	cos_x = cos_series(angle_rad, terms)		78	cosine_value = series_cosine(angle_radians, terms_required)
	80	tan_x = sin_x / cos_x		79	tangent_value = sine_value / cosine_value
	81	tan_x = round_to_significant_digits(tan_x, E)		80	tangent_value = significant_digit_round(tangent_value, precision)
	82	tan_str = '{0}'.format(tan_x.normalize())		81	tangent_string = '{0}'.format(tangent_value.normalize())
	83	if '.' in tan_str:		82	if '.' in tangent_string:
	84	integer_part, fractional_part = tan_str.split('.')		83	integer_part, fractional_part = tangent_string.split('.')
	85	significant_digits = len(integer_part.lstrip('-') + fractional_p		84	significant_digits = len(integer_part.lstrip('-') + fractional_p
	>	art)		>	art)
	86	else:		85	else:
n	87	integer_part = tan_str.lstrip('-')	n	86	integer_part = tangent_string.lstrip('-')
	88	significant_digits = len(integer_part)		87	significant_digits = len(integer_part)
n	89	if significant_digits < E:	n	88	if significant_digits < precision:
	90	zeros_needed = E - significant_digits		89	zeros_needed = precision - significant_digits
	91	if '.' in tan_str:		90	if '.' in tangent_string:
	92	tan_str += '0' * zeros_needed		91	tangent_string += '0' * zeros_needed
	93	else:		92	else:
n	94	tan_str += '.' + '0' * zeros_needed	n	93	tangent_string += '.' + '0' * zeros_needed
	95	print(tan_str)		94	print(tangent_string)
	96	if __name__ == '__main__':		95	if __name__ == '__main__':
t	97	main()	t	96	perform_calculations()

Legends
Colors  Added  Changed Deleted	Links (f)irst change (n)ext change (t)op