後ろから見ていく.

後ろから \(i\) 桁目まで見たときに \(13\) で割った余りが \(j\) である組み合わせの数を \(C(i, j)\) とすると,

\[C(i, (j+k)*10^{i-1}) \mod 13) += C(i-1, k)\]

として DP で更新していく. \(i\) 桁目が ? ならば \(j\) には \(0\) から \(9\) まで入れてループを回し, 特定の数値ならば \(j\) はその数値のみが有効となる.

\(10^{i-1}\) の部分は \(i\) が増えるたびに \(10\) 倍して更新していけばいい.