Luogu 1092(DFS)(虫食算)(NOIP 2004)

传送门

题意:补全一个成立的加法竖式,每一个大写字母代表一个n进制数且不同大写字母代表的数字一定不同,(字母与数值为双射)

题解:每一位相当于一个模方程,第一反应是高斯消元但是忘了咋写于是考虑搜索。由于要尽可能降低复杂度:①选择从低位往高位搜②对于每一位,尽量由已知数算出未知数,而不是靠枚举(直到枚举到(a+b+nxt)%n=c为止)

100多行的dfs函数一遍写下来只有一个小笔误得了80分也算是对得起这把老骨头了,暴力地讨论每一位可能的八种情况(三个数分别已知/未知),虽然代码冗长但是结构清晰?(狗头保命)

P.S.不知道第九个点是什么毒瘤构造的数据,对于未知数位从0到n-1枚举会TLE,反向枚举直接秒过......

#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;
const int N=30;
int n;
char a[N],b[N],c[N];
int mp[100];
int nxt[N];
bool vis[N];
inline void dfs(int p) {
	if (!p) {
		if (!nxt[p]) {
			for (int i=0;i<n;++i)
				printf("%d ",mp[i+'A']);
			exit(0);
		}
		else return ;
	}
	bool f1=(mp[a[p]]==-1),f2=(mp[b[p]]==-1),f3=(mp[c[p]]==-1);
	if (!f3) {//the sum of the two digit is known
		if (!f1&&!f2) {
			if ((mp[a[p]]+mp[b[p]]+nxt[p])%n==mp[c[p]]) {
				nxt[p-1]=(mp[a[p]]+mp[b[p]]+nxt[p])/n;
				dfs(p-1);
				nxt[p-1]=0;
			}
		} else if (!f1&&f2) {//fill b
			int t=(mp[c[p]]-mp[a[p]]-nxt[p]+n)%n;
			if (!vis[t]) {
				mp[b[p]]=t;
				vis[t]=true;
				nxt[p-1]=(mp[a[p]]+mp[b[p]]+nxt[p])/n;
				dfs(p-1);
				mp[b[p]]=-1;
				vis[t]=false;
				nxt[p-1]=0;
			}
		} else if (f1&&!f2) {//fill a
			int t=(mp[c[p]]-mp[b[p]]-nxt[p]+n)%n;
			if (!vis[t]) {
				mp[a[p]]=t;
				vis[t]=true;
				nxt[p-1]=(mp[a[p]]+mp[b[p]]+nxt[p])/n;
				dfs(p-1);
				mp[a[p]]=-1;
				vis[t]=false;
				nxt[p-1]=0;
			}
		} else {//fill a,b
			for (int i=n-1;~i;--i) {
				int t=(mp[c[p]]-i-nxt[p]+n)%n;
				if (!vis[i]&&!vis[t]&&(a[p]==b[p])==(i==t)) {
					mp[a[p]]=i;
					mp[b[p]]=t;
					vis[i]=vis[t]=true;
					nxt[p-1]=(mp[a[p]]+mp[b[p]]+nxt[p])/n;
					dfs(p-1);
					mp[a[p]]=-1;
					mp[b[p]]=-1;
					vis[i]=vis[t]=false;
					nxt[p-1]=0;
				}
			}
		}
	} else {//the sum of the two digit is unknown
		if (!f1&&!f2) {//fill c
			int t=(mp[a[p]]+mp[b[p]]+nxt[p])%n;
			if (!vis[t]) {
				mp[c[p]]=t;
				vis[t]=true;
				nxt[p-1]=(mp[a[p]]+mp[b[p]]+nxt[p])/n;
				dfs(p-1);
				mp[c[p]]=-1;
				vis[t]=false;
				nxt[p-1]=0;
			}
		} else if (f1&&!f2) {//fill a,c
			for (int i=n-1;~i;--i) {
				int t=(mp[b[p]]+i+nxt[p])%n;
				if (!vis[i]&&!vis[t]&&(a[p]==c[p])==(i==t)) {
					mp[a[p]]=i;
					mp[c[p]]=t;
					vis[i]=vis[t]=true;
					nxt[p-1]=(mp[a[p]]+mp[b[p]]+nxt[p])/n;
					dfs(p-1);
					mp[a[p]]=-1;
					mp[c[p]]=-1;
					vis[i]=vis[t]=false;
					nxt[p-1]=0;
				}
			}
		} else if (!f1&&f2) {//fill b,c
			for (int i=n-1;~i;--i) {
				int t=(mp[a[p]]+i+nxt[p])%n;
				if (!vis[i]&&!vis[t]&&(b[p]==c[p])==(i==t)) {
					mp[b[p]]=i;
					mp[c[p]]=t;
					vis[i]=vis[t]=true;
					nxt[p-1]=(mp[a[p]]+mp[b[p]]+nxt[p])/n;
					dfs(p-1);
					mp[b[p]]=-1;
					mp[c[p]]=-1;
					vis[i]=vis[t]=false;
					nxt[p-1]=0;
				}
			}
		} else {//fill a,b,c
			for (int i=n-1;~i;--i)
				for (int j=n-1;~j;--j) {
					if (vis[i]||vis[j]||(a[p]==b[p])^(i==j)) continue;
					int t=(i+j+nxt[p])%n;
					if (vis[t]||(a[p]==c[p])^(i==t)||(b[p]==c[p])^(j==t)) continue;
					mp[a[p]]=i;
					mp[b[p]]=j;
					mp[c[p]]=t;
					vis[i]=vis[j]=vis[t]=true;
					nxt[p-1]=(mp[a[p]]+mp[b[p]]+nxt[p])/n;
					dfs(p-1);
					mp[a[p]]=-1;
					mp[b[p]]=-1;
					mp[c[p]]=-1;
					vis[i]=vis[j]=vis[t]=false;
					nxt[p-1]=0;
				}
		}
	}
}
int main() {
//	freopen("in.txt","r",stdin);
	memset(mp,-1,sizeof(mp));
	scanf("%d",&n);
	scanf("%s%s%s",a+1,b+1,c+1);
	dfs(n);
	return 0;
}

 

 

 

 

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值