关于数位 dp 的 dfs 和 for 写法

ll dfs(ll pos, ll pre1, ll pre2, bool finish, bool not_zero, bool up) 
{
    if (pos > n) return finish;
    if (!up && ~f[pos][pre1][pre2][finish]) return f[pos][pre1][pre2][finish];
    ll end = up ? num[pos] : 9, ans = 0;
    for (int i = 0; i <= end; ++i) 
        ans = (ans + dfs(pos + 1, i, not_zero ? pre1 : -1, 
            finish || (i == pre1 || i == pre2) && not_zero, not_zero || i, up && end == i)) % mod;
    if (!up && not_zero && ~pre1 && ~pre2) f[pos][pre1][pre2][finish] = ans;
    return ans;
}

void _(int&x,int y){x+=x+y<P?y:y-P;}
int _(char*s){
	int n=strlen(s);
	if(n<2)return 0;
	for(i=n;i--;)for(j=10;j--;)for(k=10;k--;f[i][j][k]=g[i][j][k]=p[i][j][k]=q[i][j][k]=0);
	for(x=*s^48,y=s[1]^48,j=0;++j^x;)for(k=10;k--;)(j^k?p:f)[1][j][k]=1;
	for(k=y;k--;(x^k?p:f)[1][x][k]=1);
	for(i=1,(x^y?q:g)[1][x][y]=1;++i^n;){
		for(j=10;j--;)for(k=10;k--;)for(l=10;l--;)_(f[i][j][k],f[i-1][l][j]),_((l^k&&j^k?p:f)[i][j][k],p[i-1][l][j]);
		for(x=s[i-2]^48,y=s[i-1]^48,j=z=s[i]^48;j--;)_(f[i][y][j],g[i-1][x][y]),_((x^j&&y^j?p:f)[i][y][j],q[i-1][x][y]);
		for(j=10,_(g[i][y][z],g[i-1][x][y]),_((x^z&&y^z?q:g)[i][y][z],q[i-1][x][y]);--j;)for(k=10;k--;)_((j^k?p:f)[i][j][k],1);
	}
	for(S=0,j=10;j--;)for(k=10;k--;_(S,g[n-1][j][k]))_(S,f[n-1][j][k]);return S;
}