n = 73542616560647877565544036788738025202939381425158737721544398356851787401183516163221837013929559568993844046804187977705376289108065126883603562904941748653607836358267359664041064708762154474786168204628181667371305788303624396903323216279110685399145476916585122917284319282272004045859138239853037072761
p_given = 108294440701045353595867242719660522374526250640690193563048263854806748525172379331 * (10^(41+30)) + 341078269246532299656864881223

inv_delta = pow(10^30, -1, n)
print(inv_delta)

PR.<x> = PolynomialRing(Zmod(n))
f = p_given * inv_delta + x

ans = f.small_roots(X=10^41, beta=0.4)
print(ans)