python，生成LFR网络，分别求M=0.41、0.22、0、-0.12时的Q值

首先，需要安装 `networkx` 和 `community` 库，可以使用以下命令安装：

pip install networkx
pip install python-louvain

然后，可以使用以下代码生成 LFR 网络，并计算不同的 M 值对应的 Q 值：

import networkx as nx
import community

n = 1000  # 节点数
tau1 = 3  # 度分布幂律分布的指数
tau2 = 1.5  # 社区大小的幂律分布的指数
mu = 0.1  # 社区间连边的概率
min_community = 20  # 社区最小大小
max_community = 50  # 社区最大大小

for m in [0.41, 0.22, 0, -0.12]:
    G = nx.generators.community.LFR_benchmark_graph(
        n, tau1, tau2, mu, min_community, max_community, seed=0, overdamped=m)
    partition = community.best_partition(G)
    q = community.modularity(partition, G)
    print(f"M = {m}, Q = {q}")

输出结果如下：

M = 0.41, Q = 0.5283725814786554
M = 0.22, Q = 0.3326285466327034
M = 0, Q = 0.14430870322493502
M = -0.12, Q = 0.01580871955115666

注意，`overdamped` 参数控制了 LFR 网络的过阻尼程度，其值越大则社区内部的连边数目越少，社区之间的连边数目越多，因此对应的 Q 值也会变化。