from math import gcd
import numpy as np
import networkx as nx
import matplotlib.pyplot as plt

np.random.seed(123)

n = 1000  # number of nodes
p = (np.log(n) + 1) / n  # probability of edge creation
p = 1 / n
# p = 1 / n**2

G = nx.erdos_renyi_graph(n, p)

nx.draw(G, node_size=10)
plt.savefig("erdos_renyi.png", dpi=300)