「状压DP」「暴力搜索」排列perm

Posted 614685877--aakennes

tags:

篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了「状压DP」「暴力搜索」排列perm相关的知识,希望对你有一定的参考价值。

「状压DP」「暴力搜索」排列

题目描述:

题目描述

给一个数字串 s 和正整数 d, 统计 sss 有多少种不同的排列能被 d 整除(可以有前导 0)。例如 123434 有 90 种排列能被 2 整除,其中末位为 2 的有 30 种,末位为 4 的有 60 种。

输入格式

输入第一行是一个整数 TTT,表示测试数据的个数,以下每行一组 s 和 d,中间用空格隔开。s 保证只包含数字 0,1,2,3,4,5,6,7,8,9

输出格式

每个数据仅一行,表示能被 d 整除的排列的个数。

输入输出样例

输入 #1

7
000 1
001 1
1234567890 1
123434 2
1234 7
12345 17
12345678 29

输出 #1

1
3
3628800
90
3
6
1398

说明/提示

100% 的数据满足:s 的长度不超过 10,1≤d≤1000,1≤T≤15。

在前三个例子中,排列分别有 1,3,36288001 种,它们都是 1 的倍数。

解法1:状压DP

思路:

s的长度很短,不是暴搜就是状压,然鹅这道题都可以用

大多数状压中都是当前某个状态对之后于此相关状态产生影响

因此可以对整个序列的长度进行状压,详见代码

代码:

#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
#include<cmath>
#include<stack>
using namespace std;
const int maxn=(1<<10)+5,INF=0x3f3f3f3f;
int n,m,f[maxn][1000],sum[11],mol,a[maxn],vis[maxn];
char str[1000];
inline int read(){
	int s=0,w=1;
	char ch=getchar();
	while(ch<‘0‘||ch>‘9‘){if(ch==‘-‘)w=-1;ch=getchar();}
	while(ch>=‘0‘&&ch<=‘9‘)s=s*10+ch-‘0‘,ch=getchar();
	return s*w;
}
int main(){
	freopen("a.in","r",stdin);
	int t=read();
	while(t--){
		memset(f,0,sizeof(f));
		memset(a,0,sizeof(a));
		memset(sum,0,sizeof(sum));
		scanf("%s",str);
		int n=strlen(str);
		mol=read();
		for(int i=1;i<=n;i++){
			a[i]=str[i-1]-‘0‘;
		}
		int maxs=(1<<n)-1;
		f[0][0]=1;//0序列为000....的方案肯定只有一种
		for(int s=0;s<=maxs;s++){//s代表所选的数的状态,00101代表选了第一个数和第三个数的状态
			memset(vis,0,sizeof(vis));//vis记录这个数是否被访问过,因为相同的数会被重复计入方案,
			                          //如,对11230排列,11230和11230会被重复计算,用vis就可以去重,当然也可以不用vis用数学方法去重,即记录所有数出现的个数,最后ans/=cnt[i]!(i=0~9)
                        for(int i=1;i<=n;i++){   //i表示下一个要选的数
				if(s&(1<<(i-1))||vis[a[i]])continue;//s&(i<<(i-1))表示当前枚举的数被包含在了当前状态s里面,不用再次计算
				vis[a[i]]=1;
				for(int k=0;k<mol;k++){//枚举所有余数k,f[s][k]代表所选的数的状态为s,且余数为k时的总排列数
					f[s|(1<<(i-1))][(k*10+a[i])%mol]+=f[s][k];//解释:因为我们的状态是从小到大枚举的,所以f[s][k]会对
                                         //它的所有下一个状态产生影响,而下一个状态的余数恰恰是(k*10+a[i])%mol
                                         //举个例子,序列1234,模数是4,当前状态0011代表12的排列,12会转移到124和123即 原数*10+a[i],相应的余数就会变成 (余数*10+a[i])%mol,即0和3
					//	cout<<f[s|(1<<(i-1))][(k*10+a[i])%mol]<<endl;
				}
			}
		}
		cout<<f[maxs][0]<<endl;
	}
}

解法2:暴搜

相对于状压而言,暴搜更好想一些,就是从低位依次枚举至高位,但是时间消耗更大,洛谷上会T两个点,可能剪剪枝会过

#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
#include<cmath>
#include<stack>
using namespace std;
const int maxn=(1<<10)+5,INF=0x3f3f3f3f;
int n,m,f[maxn][1000],sum[11],mol,a[maxn],vis[maxn],ans;
char str[1000];
inline int read(){
	int s=0,w=1;
	char ch=getchar();
	while(ch<‘0‘||ch>‘9‘){if(ch==‘-‘)w=-1;ch=getchar();}
	while(ch>=‘0‘&&ch<=‘9‘)s=s*10+ch-‘0‘,ch=getchar();
	return s*w;
}
void DFS(int now,long long x){//看数据范围,不开long long会炸,x代表当前的数,now代表当前位数
	if(now>n){//now>n说明x是末状态
		if(x%mol==0)ans++;
		return;
	}
	for(int i=0;i<=9;i++){
		if(sum[i]){
			sum[i]--;//为保证当前这一位选取这个元素对这一位选取其他元素没影响,所以自减后要自加回来
			DFS(now+1,x*10+i);
			sum[i]++;
		}
	}
}
int main(){
	freopen("a.in","r",stdin);
	int t=read();
	while(t--){
		memset(a,0,sizeof(a));
		memset(sum,0,sizeof(sum));
		string str;
		cin>>str;
		mol=read();
		n=str.size();
		ans=0;
		for(int i=0;i<str.size();i++){
			a[i]=str[i]-‘0‘;
			sum[a[i]]++;
		}
		DFS(1,0);
		cout<<ans<<endl;
	}
}

解法3:next_permutation生成全排列

应该所有人都想过用全排列来写,但是时间开销很大,吾日观洛谷,发现STL中的几个比较有趣的函数:

next_permutation:从原递增序列中求出所有全排列

prev_permutation:从原递减序列中求出所有全排列

atol:将字符串转换为数列

详见代码

代码1:不用atol正常写

#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
#include<cmath>
#include<stack>
using namespace std;
const int maxn=(1<<10)+5,INF=0x3f3f3f3f;
int n,m,f[maxn][1000],sum[11],mol,a[maxn],vis[maxn],ans;
char str[1000];
inline int read(){
	int s=0,w=1;
	char ch=getchar();
	while(ch<‘0‘||ch>‘9‘){if(ch==‘-‘)w=-1;ch=getchar();}
	while(ch>=‘0‘&&ch<=‘9‘)s=s*10+ch-‘0‘,ch=getchar();
	return s*w;
}
int main(){
	freopen("a.in","r",stdin);
	int t=read();
	while(t--){
		memset(a,0,sizeof(a));
		memset(sum,0,sizeof(sum));
		string str;
		cin>>str;
		mol=read();
		n=str.size();
		ans=0;
		for(int i=0;i<n;i++){
			a[i]=str[i]-‘0‘;
		}
		sort(a,a+n);//严格要求必须递增,否则全排列不对
		do{
			long long temp=0;
			for(int i=0;i<n;i++)temp=temp*10+a[i];//全排列是保存在数组a里的,通过这种方法取出来
			if(temp%mol==0)ans++;
		}while(next_permutation(a,a+n));
		cout<<ans<<endl;
	}
}

代码2:用atol

#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
#include<cmath>
#include<stack>
using namespace std;
const int maxn=(1<<10)+5,INF=0x3f3f3f3f;
int n,m,f[maxn][1000],sum[11],mol,a[maxn],vis[maxn],ans;
char str[1000];
inline int read(){
	int s=0,w=1;
	char ch=getchar();
	while(ch<‘0‘||ch>‘9‘){if(ch==‘-‘)w=-1;ch=getchar();}
	while(ch>=‘0‘&&ch<=‘9‘)s=s*10+ch-‘0‘,ch=getchar();
	return s*w;
}
int main(){
	freopen("a.in","r",stdin);
	int t=read();
	while(t--){
		memset(a,0,sizeof(a));
		memset(sum,0,sizeof(sum));
		char str[maxn];//string排序排的是字符串,char排的是字符
		cin>>str;
		mol=read();
		n=strlen(str);
		ans=0;
		sort(str,str+n);
		do{
			long long temp=atoll(str);//直接转
			
			if(temp%mol==0)ans++;
		}while(next_permutation(str,str+n));
		cout<<ans<<endl;
	}
}













以上是关于「状压DP」「暴力搜索」排列perm的主要内容,如果未能解决你的问题,请参考以下文章

[Bzoj1072][SCOI2007]排列perm(状压dp)

BZOJ1072: [SCOI2007]排列perm 状压DP

[BZOJ1072][SCOI2007]排列perm 状压dp

B1072 [SCOI2007]排列perm 状压dp

bzoj 1072: [SCOI2007]排列perm 状压dp

BZOJ 1072: [SCOI2007]排列perm [DP 状压 排列组合]