Let , , be natural numbers satisfying
Let
Evaluate
If , compute .
The problem asks us to evaluate an integral involving a polynomial determined by the roots , which satisfy a system of power sums.
First, we determine the coefficients of the polynomial . Using Newton sums:
Then we evaluate the integral using the residue theorem, where . Since is in the upper half-plane, the integral is . The residue for a pole of order 4 is .
Let's perform the calculations.
[PYTHON] import cmath
Given power sums
p1 = 151 p2 = 10939 p3 = 957871
Calculate elementary symmetric polynomials
e1 = p1 e2 = (e1 * p1 - p2) // 2 e3 = (e2 * p1 - e1 * p2 + p3) // 3
print(f"{e1=}, {e2=}, {e3=}")
Define Q and its derivatives at a point z
Q(z) = z^3 - e1z^2 + e2z - e3
def get_Q_derivs(z): # Q(z) Q0 = z3 - e1*z2 + e2z - e3 # Q'(z) = 3z^2 - 2e1z + e2 Q1 = 3z**2 - 2e1z + e2 # Q''(z) = 6z - 2e1 Q2 = 6z - 2*e1 # Q'''(z) = 6 Q3 = 6 return Q0, Q1, Q2, Q3
The pole z0
z0 = 205/69 * 1j
Calculate the term from the 3rd derivative of Q(z)e^{iz}
g'''(z) = [Q'''(z) + 3i Q''(z) - 3 Q'(z) - i Q(z)] e^{iz}
Q0, Q1, Q2, Q3 = get_Q_derivs(z0) term_bracket = Q3 + 3j * Q2 - 3 * Q1 - 1j * Q0 exp_term = cmath.exp(1j * z0)
Integral I = 2pii * (1/6) * term_bracket * exp_term
I = (2 * cmath.pi * 1j / 6) * term_bracket * exp_term
alpha = I.real beta = I.imag
print(f"{I=}") print(f"{alpha=}") print(f"{beta=}")
Compute final answer
ans = int(10 * alpha) + int(beta) # Using int() acts like floor for positive, but check signs import math ans_floor = math.floor(10 * alpha) + math.floor(beta) print(f"{ans_floor=}") [/PYTHON]
[THOUGHT]
The Python code calculated the elementary symmetric polynomials correctly from the power sums.
The integral was evaluated using the residue theorem. The pole is at . Since the imaginary part is positive, the pole is inside the contour (assuming standard semi-circle in upper half plane).
The residue calculation formula used was .
The code output .
The quantity to compute is .
.
.
Sum = .
The code output ans_floor = -36330 confirms this.
Wait, let me double check the int() vs floor() logic in the code output.
The code printed ans_floor=-36330.
math.floor(-36320.72) is -36321.
math.floor(-8.919) is -9.
Sum is -36330.
The code output is consistent.
Final check on the integral formula.
? No.
The integrand is .
Let . The residue at (order 4) is .
So .
The code implemented exactly this.
.
Using Leibniz rule or direct expansion:
.
The code used term_bracket = Q3 + 3j * Q2 - 3 * Q1 - 1j * Q0.
This matches .
The factor was included.
The factor was included.
The calculation seems correct.
One detail: . . This is a real number. coefficients are real. is purely imaginary. Let . . . . .
Term bracket . Real part of B: . Imaginary part of B: .
Then . . So . .
The code output:
.
.
The signs and magnitudes are handled by the code.
The result ans_floor is -36330.
[/THOUGHT]
\boxed{-36330}