max.l=1NAl=l=1Ni{dlid0}NG(dli,d0)Ri=l=1Ni{dlid0}NG(dli,d0)Wi+Oil{dild0}G(dil,d0)Plmax. \quad \sum_{l=1}^N A_l=\sum_{l=1}^N\sum_{i\in \{d_{li}\le d_0\}}^NG(d_{li},d_0)R_i \\ =\sum_{l=1}^N\sum_{i\in \{d_{li}\le d_0\}}^NG(d_{li},d_0)\frac{W_i+O_i}{\sum_{l\in \{d_{il}\le d_0\}}G(d_{il},d_0)P_l}

其中NN表示网格总数,OiO_i表示格网ii 原始的服务设施能力总和,WiW_i为格网ii 更新的服务设施的服务能力总和,计算公式为:

Wi=TxiW_i= \sum{T}\cdot x_{i}

其中,xix_i 表示数量矩阵XX在第ii 个格网的分量,记作:xi=[xi1,,xij]x_i=[x_{i1},\cdots,x_{ij}],TT表示修正后的设施价值参数,记作:

T=[T1,,Tj]T=[T_1,\cdots,T_j]

(2)(2)中,xix_i的原始矩阵XX记作:

X=x11x1jxi1xijX=\left| \begin{matrix} x_{11} & \cdots & x_{1j} \\ \vdots &\ddots & \vdots \\ x_{i1} & \cdots & x_{ij} \\ \end{matrix} \right|

矩阵XX 记录设施的数量分配信息,向量TT 记录设施的服务能力,他们蕴含了问题中自变量的最基本信息。