Advertisement
PyTimur

Untitled

Apr 16th, 2021
53
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 1.14 KB | None | 0 0
  1. def main():
  2. def f(n):
  3. d = 1; k = r = 0
  4. while d*d < n:
  5. if n % d == 0:
  6. k += 2
  7. if isprime(d) == 1: r += 1
  8. if isprime(n//d) == 1: r += 1
  9. d += 1
  10. if d*d == n:
  11. k += 1
  12. if isprime(d) == 1: r += 1
  13. if (k == 4) and (r == 2):
  14. return 1
  15. return 0
  16.  
  17. def isprime(n):
  18. d = 2;k = 0
  19. while d*d < n:
  20. if n % d == 0:
  21. k += 2
  22. return 0
  23. d += 1
  24. if d*d == n: return 0
  25. return 1
  26.  
  27. t = r = 0
  28. for i in range(1,1000000):
  29. if int(round(pow(i,1/3)))**3 != i:
  30. continue
  31. if f(i) == 1:
  32. t += 1
  33. if t % 3 == 1: r += i
  34. a = [];d = 1
  35. while d*d < r:
  36. if r % d == 0:
  37. a.append(d)
  38. a.append((r//d))
  39. d += 1
  40. if d*d == r: a.append(d)
  41. for i in range(len(a)-1):
  42. for j in range(i, len(a)-1):
  43. if a[j] > a[j + 1]:
  44. a[j], a[j + 1] = a[j + 1], a[j]
  45. print(a)
  46. if __name__ == '__main__':
  47. main()
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement