求n的数字重排且整除m的数（无前导0）
这道题我的写法挺鬼畜（用了一直想写却从未成功的11维动态数组）
dp[x][st][……]表示前x位，原数位st，每种数字的个数为……的数。
这里注意一下，0-9一共10个数字，但是只要9维状态就行，因为剩下的可以用x和其他数表示，所以也是确定的，不用写进状态了。
（常数级别的优化，但是这里不优化空间会被卡常）。
下附鬼畜代码TT:
1 #include<bits/stdc++.h> 2 #define ll long long 3 using namespace std; 4 int c[10]; 5 ll n,m; 6 int len; 7 ll ***********dp; 8 ll dfs(int x,int st,int fz,int c0,int c1,int c2,int c3,int c4,int c5,int c6, int c7,int c8){ 9 if (x==-1){ 10 return (st==0); 11 } 12 if (!fz && dp[x][st][c0][c1][c2][c3][c4][c5][c6][c7][c8]!=-1) return dp[x][st][c0][c1][c2][c3][c4][c5][c6][c7][c8]; 13 int maxn=9; 14 ll ret=0; 15 for (int i=0; i<=maxn; i++){ 16 if (fz){ 17 switch (i){ 18 case 1:{if (c1>0) ret+=dfs(x-1,(st*10+i)%m,0,c0,c1-1,c2,c3,c4,c5,c6,c7,c8); break;} 19 case 2:{if (c2>0) ret+=dfs(x-1,(st*10+i)%m,0,c0,c1,c2-1,c3,c4,c5,c6,c7,c8); break;} 20 case 3:{if (c3>0) ret+=dfs(x-1,(st*10+i)%m,0,c0,c1,c2,c3-1,c4,c5,c6,c7,c8); break;} 21 case 4:{if (c4>0) ret+=dfs(x-1,(st*10+i)%m,0,c0,c1,c2,c3,c4-1,c5,c6,c7,c8); break;} 22 case 5:{if (c5>0) ret+=dfs(x-1,(st*10+i)%m,0,c0,c1,c2,c3,c4,c5-1,c6,c7,c8); break;} 23 case 6:{if (c6>0) ret+=dfs(x-1,(st*10+i)%m,0,c0,c1,c2,c3,c4,c5,c6-1,c7,c8); break;} 24 case 7:{if (c7>0) ret+=dfs(x-1,(st*10+i)%m,0,c0,c1,c2,c3,c4,c5,c6,c7-1,c8); break;} 25 case 8:{if (c8>0) ret+=dfs(x-1,(st*10+i)%m,0,c0,c1,c2,c3,c4,c5,c6,c7,c8-1); break;} 26 case 9:{if ((x+1-c0-c1-c2-c3-c4-c5-c6-c7-c8)>0) ret+=dfs(x-1,(st*10+i)%m,0,c0,c1,c2,c3,c4,c5,c6,c7,c8); break;} 27 } 28 } 29 else { 30 switch (i){ 31 case 0:{if (c0>0) ret+=dfs(x-1,(st*10+i)%m,0,c0-1,c1,c2,c3,c4,c5,c6,c7,c8); break;} 32 case 1:{if (c1>0) ret+=dfs(x-1,(st*10+i)%m,0,c0,c1-1,c2,c3,c4,c5,c6,c7,c8); break;} 33 case 2:{if (c2>0) ret+=dfs(x-1,(st*10+i)%m,0,c0,c1,c2-1,c3,c4,c5,c6,c7,c8); break;} 34 case 3:{if (c3>0) ret+=dfs(x-1,(st*10+i)%m,0,c0,c1,c2,c3-1,c4,c5,c6,c7,c8); break;} 35 case 4:{if (c4>0) ret+=dfs(x-1,(st*10+i)%m,0,c0,c1,c2,c3,c4-1,c5,c6,c7,c8); break;} 36 case 5:{if (c5>0) ret+=dfs(x-1,(st*10+i)%m,0,c0,c1,c2,c3,c4,c5-1,c6,c7,c8); break;} 37 case 6:{if (c6>0) ret+=dfs(x-1,(st*10+i)%m,0,c0,c1,c2,c3,c4,c5,c6-1,c7,c8); break;} 38 case 7:{if (c7>0) ret+=dfs(x-1,(st*10+i)%m,0,c0,c1,c2,c3,c4,c5,c6,c7-1,c8); break;} 39 case 8:{if (c8>0) ret+=dfs(x-1,(st*10+i)%m,0,c0,c1,c2,c3,c4,c5,c6,c7,c8-1); break;} 40 case 9:{if ((x+1-c0-c1-c2-c3-c4-c5-c6-c7-c8)>0) ret+=dfs(x-1,(st*10+i)%m,0,c0,c1,c2,c3,c4,c5,c6,c7,c8); break;} 41 } 42 } 43 } 44 return dp[x][st][c0][c1][c2][c3][c4][c5][c6][c7][c8]=ret; 45 } 46 int main(){ 47 scanf("%lld%lld",&n,&m); 48 len=0; 49 while (n){ 50 c[n%10]++; 51 n/=10; 52 len++; 53 } 54 dp=new ll **********[len]; 55 for (int i=0; i<len; i++){ 56 dp[i]=new ll*********[m]; 57 for (int i0=0; i0<m; i0++){ 58 dp[i][i0]=new ll********[c[0]+1]; 59 for (int i1=0; i1<c[0]+1; i1++){ 60 dp[i][i0][i1]=new ll*******[c[1]+1]; 61 for (int i2=0; i2<c[1]+1; i2++){ 62 dp[i][i0][i1][i2]=new ll******[c[2]+1]; 63 for (int i3=0; i3<c[2]+1; i3++){ 64 dp[i][i0][i1][i2][i3]=new ll*****[c[3]+1]; 65 for (int i4=0; i4<c[3]+1; i4++){ 66 dp[i][i0][i1][i2][i3][i4]=new ll ****[c[4]+1]; 67 for (int i5=0; i5<c[4]+1; i5++){ 68 dp[i][i0][i1][i2][i3][i4][i5]=new ll ***[c[5]+1]; 69 for (int i6=0; i6<c[5]+1; i6++){ 70 dp[i][i0][i1][i2][i3][i4][i5][i6]=new ll **[c[6]+1]; 71 for (int i7=0; i7<c[6]+1; i7++){ 72 dp[i][i0][i1][i2][i3][i4][i5][i6][i7]=new ll *[c[7]+1]; 73 for (int i8=0; i8<c[7]+1; i8++){ 74 dp[i][i0][i1][i2][i3][i4][i5][i6][i7][i8]=new ll [c[8]+1]; 75 for (int i9=0; i9<c[8]+1; i9++){ 76 dp[i][i0][i1][i2][i3][i4][i5][i6][i7][i8][i9]=-1; 77 } 78 } 79 } 80 } 81 } 82 } 83 } 84 } 85 } 86 } 87 } 88 ll res=dfs(len-1,0,1,c[0],c[1],c[2],c[3],c[4],c[5],c[6],c[7],c[8]); 89 printf("%lld",res); 90 }