Compute $I=x_sH_p(G)$.
Randomly pick $\alpha,c_1,…,c_{s-1},c_{s+1},…,c_n$ and compute:
$L=\alpha G+\sum_{i\neq s}^{n}{c_i}P_i$.
$R=\alpha H_p(G)+\sum_{i\neq s}^{n}c_iI$.
Compute $c_s=H_s(L,R,M)-\sum_{i\neq s}^nc_i$.
Compute $r=\alpha -c_sx_s$.
Output the signature $\sigma=(I,r,c_1,…,c_n)$.
Compute $L=rG+\sum_i^nc_iP_i$ and $R=rH_p(G)+\sum_i^nc_iI$.
Check whether $\sum_i^nc_i=H_s(L,R,M)$.