import random
import scipy.stats as stats
for x in range(20):
nums = []
for _ in range(4):
nums.append(random.randint(1,201))
_, two_sided_pvalue = stats.fisher_exact([[nums[0], nums[1]], [nums[2], nums[3]]], 'two-sided')
_, less_sided_pvalue = stats.fisher_exact([[nums[0], nums[1]], [nums[2], nums[3]]], 'less')
_, greater_sided_pvalue = stats.fisher_exact([[nums[0], nums[1]], [nums[2], nums[3]]], 'greater')
print("([{}, {}, {}, {}], {}, {}, {}),".format(
nums[0], nums[1], nums[2], nums[3],
less_sided_pvalue,
greater_sided_pvalue,
two_sided_pvalue))