BZOJ 1072 [SCOI2007]排列perm
1072: [SCOI2007]排列perm
Time Limit: 10 Sec Memory Limit: 162 MBSubmit: 1268 Solved: 782
[Submit][Status][Discuss]
Description
给一个数字串s和正整数d, 统计s有多少种不同的排列能被d整除(可以有前导0)。例如123434有90种排列能被2整除,其中末位为2的有30种,末位为4的有60种。
Input
输入第一行是一个整数T,表示测试数据的个数,以下每行一组s和d,中间用空格隔开。s保证只包含数字0, 1, 2, 3, 4, 5, 6, 7, 8, 9.
Output
每个数据仅一行,表示能被d整除的排列的个数。
Sample Input
7
000 1
001 1
1234567890 1
123434 2
1234 7
12345 17
12345678 29
000 1
001 1
1234567890 1
123434 2
1234 7
12345 17
12345678 29
Sample Output
1
3
3628800
90
3
6
1398
3
3628800
90
3
6
1398
HINT
在前三个例子中,排列分别有1, 3, 3628800种,它们都是1的倍数。
【限制】
100%的数据满足:s的长度不超过10, 1<=d<=1000, 1<=T<=15
Source
题解:隐藏了一下的"n小思状压",f[S][i]表示状态为S,模d为i的数量,在S后面补东西来转移。注意,如果原数中有重复的,那么算出来也就会有重复的,用最后的答案除以cnt[i]的阶乘就可以了。
1 #include<iostream> 2 #include<cstdio> 3 #include<cmath> 4 #include<algorithm> 5 #include<stack> 6 #include<queue> 7 #include<cstring> 8 #define PAU putchar(' ') 9 #define ENT putchar('\n') 10 #define MSE(a,b) memset(a,b,sizeof(a)) 11 #define REN(x) for(ted*e=fch[x];e;e=e->nxt) 12 #define TIL(x) for(int i=1;i<=x;i++) 13 #define ALL(x) for(int j=1;j<=x;j++) 14 using namespace std; 15 const int maxn=1000+10,maxd=(1<<10)+10; 16 int f[maxd][maxn],cnt[10]; 17 inline int read(){ 18 int x=0;bool sig=true;char ch=getchar(); 19 for(;!isdigit(ch);ch=getchar())if(ch=='-')sig=false; 20 for(;isdigit(ch);ch=getchar())x=10*x+ch-'0';return sig?x:-x; 21 } 22 inline void write(int x){ 23 if(x==0){putchar('0');return;}if(x<0)putchar('-'),x=-x; 24 int len=0;static int buf[20];while(x)buf[len++]=x%10,x/=10; 25 for(int i=len-1;i>=0;i--)putchar(buf[i]+'0');return; 26 } 27 int T,d,lim,len,A[10],ans;char s[10]; 28 int main(){ 29 T=read(); 30 while(T--){ 31 scanf("%s",s);d=read();MSE(f,0);MSE(cnt,0); 32 for(int i=0;s[i];i++)cnt[A[i]=s[i]-'0']++;lim=1<<(len=strlen(s)); 33 f[0][0]=1; 34 for(int S=0;S<lim;S++) 35 for(int i=0;i<d;i++) 36 for(int j=0;j<len;j++)if(!(S&(1<<j))) 37 f[S|(1<<j)][(10*i+A[j])%d]+=f[S][i]; 38 ans=f[lim-1][0]; 39 for(int i=0;i<10;i++) 40 for(int j=1;j<=cnt[i];j++) 41 ans/=j; 42 write(ans);ENT; 43 } 44 return 0; 45 }