Noip提高组 2004 T4 虫食算搜索优化-CSDN博客

本文链接：https://blog.csdn.net/ThinFatty/article/details/78221507

题目描述

所谓虫食算，就是原先的算式中有一部分被虫子啃掉了，需要我们根据剩下的数字来判定被啃掉的字母。来看一个简单的例子：

http://paste.ubuntu.com/25448822/

其中#号代表被虫子啃掉的数字。根据算式，我们很容易判断：第一行的两个数字分别是5和3，第二行的数字是5。

现在，我们对问题做两个限制：

首先，我们只考虑加法的虫食算。这里的加法是N进制加法，算式中三个数都有N位，允许有前导的0。

其次，虫子把所有的数都啃光了，我们只知道哪些数字是相同的，我们将相同的数字用相同的字母表示，不同的数字用不同的字母表示。如果这个算式是N进制的，我们就取英文字母表午的前N个大写字母来表示这个算式中的0到N-1这N个不同的数字：但是这N个字母并不一定顺序地代表0到N-1)。输入数据保证N个字母分别至少出现一次。

http://paste.ubuntu.com/25448824/

上面的算式是一个4进制的算式。很显然，我们只要让ABCD分别代表0123，便可以让这个式子成立了。你的任务是，对于给定的N进制加法算式，求出N个不同的字母分别代表的数字，使得该加法算式成立。输入数据保证有且仅有一组解

输入输出格式

输入格式：

包含四行。第一行有一个正整数N(N<=26)，后面的3行每行有一个由大写字母组成的字符串，分别代表两个加数以及和。这3个字符串左右两端都没有空格，从高位到低位，并且恰好有N位。

输出格式：

包含一行。在这一行中，应当包含唯一的那组解。解是这样表示的：输出N个数字，分别表示A，B，C……所代表的数字，相邻的两个数字用一个空格隔开，不能有多余的空格。

输入输出样例

输入样例#1：

5
ABCED
BDACE
EBBAA

输出样例#1：

1 0 3 4 2

说明

对于30％的数据，保证有N<＝10；

对于50％的数据，保证有N<＝15；

对于全部的数据，保证有N<＝26。

noip2004提高组第4题

传送门

先挂个题解当做自己ac了= =其实vijos上测试ac，而洛谷却tle了……

不想说什么了，洛谷的评测机有时候真的不行。。

这题官方正解是高斯消元……然而我被“搜索过此题”的思想感染……

然后考试的时候就打了个爆搜……

事实上有很多优化，从后往前一位位搜索，

比如当前位，已知两个数，已知一个数，或者有两个字母一样等等情况都可以用作剪枝。

然后如果人品ok答案不远的话还可以快速出解。

打得很长……剪枝很多但是常数自然也大了= =提交到洛谷上发现80

猜想了一下觉得肯定A最大，B次大这种数据在构造，所以就从n-1到0枚举这一位填什么，

然后就很快了，拿了90.。

剩下那个点弄下来却发现0.8s可以跑出答案……woc？？

……于是就去vijos上交然后发现ac了。不想说什么了= =

#include<bits/stdc++.h>
using namespace std;
int n,zh[30],p[30];
char A[30],B[30],C[30];
bool flag;
void print_ans(){
	flag=1;
	for (int i=0;i<n-1;i++) printf("%d ",zh[i]);
	printf("%d\n",zh[n-1]);
}
void dfs(int now,int plus){
	if (flag) return;
	if (!now){
		if (!plus) print_ans();
		return;
	}
	int t1=zh[A[now]-'A'],t2=zh[B[now]-'A'],t3=zh[C[now]-'A'];
	if (t1<0 && t2<0 && t3<0){
		if (A[now]==B[now]){
			if (A[now]!=C[now]){
				for (int i=n-1;i>=0;i--)
					if (!p[i]){
						p[i]=1,zh[A[now]-'A']=i,zh[C[now]-'A']=(i+i+plus)%n;
						dfs(now-1,(i+i+plus)/n);
						p[i]=0,zh[A[now]-'A']=t1,zh[C[now]-'A']=t2;
					}
			} else{
				for (int i=n-1;i>=0;i--)
					if (!p[i] && (i+i+plus)%n==i){
						p[i]=1,zh[A[now]-'A']=i;
						dfs(now-1,(i+i+plus)/n);
						p[i]=0,zh[A[now]-'A']=t1;
					}
			}
		} else{
			if (A[now]==C[now]){
				if (!plus){
					if (p[0]) return;
					p[0]=1,zh[B[now]-'A']=0;
					for (int i=n-1;i;i--)
						if (!p[i]){
							p[i]=1,zh[A[now]-'A']=i;
							dfs(now-1,0);
							p[i]=0,zh[A[now]-'A']=t1;
						}
					p[0]=0,zh[B[now]-'A']=t2;
				} else{
					if (p[n-plus]) return;
					p[n-plus]=1,zh[B[now]-'A']=n-plus;
					for (int i=n-1;i>=0;i--)
						if (!p[i]){
							p[i]=1,zh[A[now]-'A']=i;
							dfs(now-1,1);
							p[i]=0,zh[A[now]-'A']=t1;
						}
					p[n-plus]=0,zh[B[now]-'A']=t2;
				}
			} else
			if (B[now]==C[now]){
				if (!plus){
					if (p[0]) return;
					p[0]=1,zh[A[now]-'A']=0;
					for (int i=n-1;i;i--)
						if (!p[i]){
							p[i]=1,zh[B[now]-'A']=i;
							dfs(now-1,0);
							p[i]=0,zh[B[now]-'A']=t1;
						}
					p[0]=0,zh[A[now]-'A']=t2;
				} else{
					if (p[n-plus]) return;
					p[n-plus]=1,zh[A[now]-'A']=n-plus;
					for (int i=n-1;i>=0;i--)
						if (!p[i]){
							p[i]=1,zh[B[now]-'A']=i;
							dfs(now-1,1);
							p[i]=0,zh[B[now]-'A']=t1;
						}
					p[n-plus]=0,zh[A[now]-'A']=t2;
				}
			} else{
				for (int i=n-1;i>=0;i--)
					if (!p[i])
					for (int j=n-1;j>=0;j--)
						if (!p[j] && !p[(i+j+plus)%n]){
							if (i==j || i==((i+j+plus)%n) || j==((i+j+plus)%n)) continue;
							p[i]=p[j]=p[(i+j+plus)%n]=1;
							zh[A[now]-'A']=i,zh[B[now]-'A']=j,zh[C[now]-'A']=(i+j+plus)%n;
							dfs(now-1,(i+j+plus)/n);
							p[i]=p[j]=p[(i+j+plus)%n]=0;
							zh[A[now]-'A']=t1,zh[B[now]-'A']=t2,zh[C[now]-'A']=t3;						
						}
			}
		}
	}
	if (~t1){
		if (~t2){
			int t=(t1+t2+plus)%n;
			if (~t3){
				if (t3!=t) return;
				dfs(now-1,(t1+t2+plus)/n);
			} else{
				if (p[t]) return;
				p[t]=1,zh[C[now]-'A']=t;
				dfs(now-1,(t1+t2+plus)/n);
				p[t]=0,zh[C[now]-'A']=t3;
			}
			return;
		} else{
			if (~t3){
				if ((t1+plus)%n==t3){
					if (p[0]) return;
					p[0]=1,zh[B[now]-'A']=0;
					dfs(now-1,0);
					p[0]=0,zh[B[now]-'A']=t2;
				} else{
					int tmp=(t3-t1-plus+n)%n;
					if (p[tmp]) return;
					p[tmp]=1,zh[B[now]-'A']=tmp;
					dfs(now-1,(t1+tmp+plus)/n);
					p[tmp]=0,zh[B[now]-'A']=t2;
				}
				return;
			} else{
				for (int i=n-1;i>=0;i--)
					if (!p[i]){
						int tmp=(i+plus+t1)%n;
						if ((i==tmp && B[now]!=C[now]) || p[tmp]) continue;
						p[i]=p[tmp]=1,zh[B[now]-'A']=i,zh[C[now]-'A']=tmp;
						dfs(now-1,(t1+i+plus)/n);
						p[i]=p[tmp]=0,zh[B[now]-'A']=t2,zh[C[now]-'A']=t3;
					}
			}
		}
	}
	if (~t2){
		if (~t3){
			if ((t2+plus)%n==t3){
				if (p[0]) return;
				p[0]=1,zh[A[now]-'A']=0;
				dfs(now-1,0);
				p[0]=0,zh[A[now]-'A']=t1;
			} else{
				int tmp=(t3-t2-plus+n)%n;
				if (p[tmp]) return;
				p[tmp]=1,zh[A[now]-'A']=tmp;
				dfs(now-1,(t2+tmp+plus)/n);
				p[tmp]=0,zh[A[now]-'A']=t1;
			}
			return;
		} else{
			for (int i=n-1;i>=0;i--)
				if (!p[i]){
					int tmp=(i+plus+t2)%n;
					if ((i==tmp && A[now]!=C[now]) || p[tmp]) continue;
					p[i]=p[tmp]=1,zh[A[now]-'A']=i,zh[C[now]-'A']=tmp;
					dfs(now-1,(t2+i+plus)/n);
					p[i]=p[tmp]=0,zh[A[now]-'A']=t1,zh[C[now]-'A']=t3;
				}
		}
	}
	if (~t3){
		for (int i=n-1;i>=0;i--)
			if (!p[i]){
				int j=t3-i-plus;
				if (j<0) j+=n;
				if ((i==j && A[now]!=B[now]) || p[j]) continue;
				p[i]=p[j]=1,zh[A[now]-'A']=i,zh[B[now]-'A']=j;
				dfs(now-1,(i+j+plus)/n);
				p[i]=p[j]=0,zh[A[now]-'A']=t1,zh[B[now]-'A']=t2;
			}
	}
}
int main(){
	scanf("%d",&n);
	scanf("%s",A+1);
	scanf("%s",B+1);
	scanf("%s",C+1);
	memset(zh,255,sizeof(zh));
	flag=0;
	dfs(n,0);
	return 0;
}