[FT][10+]NOIP 2004 虫食算

题目描述

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

43#9865#045

+8468#6633

44445509678

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

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

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

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

BADC

CBDA
DCCC 上面的算式是一个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题

---

【分析】
搞了将近三节课…下了一堆数据终于骗过去了…
QwQQQQ真的好难
思路是这样的：

从右往左，一列列搜索。搜索时传递第几列和是否进位。

对于两个加数，每当遇到一个没有匹配到的字母，就找一个数字塞进去；对于和，计算两个加数之和并mod n，如果与已有冲突就返回。

简单地这么处理就已经有80分了；如果找数字的时候从后往前找则能够到达90分。

有一个比较好想到的剪枝：如果 其中一个加数 与 和 已经确定，则可以算出 另外一个加数 。（同样，已经冲突就返回）。

但是，实际上这并没有什么卵用，大多数时候看不出效率差。（我后来把这个剪枝删掉，仍然跑了16ms）。

在上述思路的基础上，有个能大幅提高效率的剪枝：搜索当前列左边的所有列，如果三个数都已经确定（A+B=C），且无论是否进位都无法满足这个等式，就可以退掉当前搜索。

---

【代码】

//洛谷 P1092 虫食算
#include<iostream>
#include<cmath>
#include<cstdio>
#include<cstring>
#include<algorithm>
#define M(a) memset(a,0,sizeof a)
#define fo(i,j,k) for(i=j;i<=k;i++)
using namespace std;
const int mxn=1e5+5;
char a[30],b[30],c[30],x[30],y[30],z[30];
int n;
int zimu[100];
bool num[100];
inline void print()
{
    int i;
    fo(i,0,n-1) if(!num[i]) return;
    for(i=65;i<=65+n-1;i++)
      printf("%d ",zimu[i]);
    printf("\n");
    exit(0);
}
inline void dfs(int dep,int ad)
{
    if(dep>n) print();
    int i,j,jia=0,flag=zimu[c[dep]];
    fo(i,dep,n)
      if(zimu[a[i]]!=-1 && zimu[b[i]]!=-1 && zimu[c[i]]!=-1)
        if((zimu[a[i]]+zimu[b[i]])%n!=zimu[c[i]] && (zimu[a[i]]+zimu[b[i]]+1)%n!=zimu[c[i]])
          return;  //剪枝 
    if(zimu[a[dep]]!=-1 && zimu[b[dep]]!=-1 && zimu[c[dep]]!=-1)
    {
        jia=0;
        int tmp=zimu[a[dep]]+zimu[b[dep]]+ad; 
        if(tmp>=n) jia=1;
        dfs(dep+1,jia);
    }
    else if(zimu[a[dep]]!=-1 && zimu[b[dep]]!=-1)  //ab均有,c空 
    {
        jia=0;
        int tmp=zimu[a[dep]]+zimu[b[dep]]+ad;
        if(num[tmp%n]) return;
        if(tmp>=n) jia=1;
        zimu[c[dep]]=tmp%n;
        num[tmp%n]=1;
        dfs(dep+1,jia);
        num[tmp%n]=0;
        zimu[c[dep]]=flag;
        num[flag]=1;zimu[c[dep]]=flag;
    }
    else if(zimu[a[dep]]!=-1)  //b空 
    {
        fo(i,0,n-1)
          if(!num[i])
          {
              jia=0;
              int tmp=zimu[a[dep]]+i+ad;
              if(zimu[c[dep]]>=0 && tmp%n!=zimu[c[dep]]) continue;
              if(num[tmp%n] && zimu[c[dep]]!=(tmp%n)) continue;
              if(tmp>=n) jia=1;
              zimu[b[dep]]=i;
              zimu[c[dep]]=tmp%n;
              num[i]=num[tmp%n]=1;
              dfs(dep+1,jia);
              num[i]=num[tmp%n]=0;
              num[flag]=1;zimu[c[dep]]=flag;
          }
        zimu[b[dep]]=-1;
    }
    else if(zimu[b[dep]]!=-1)  //a空 
    {
        fo(i,0,n-1)
          if(!num[i])
          {
              jia=0;
              int tmp=zimu[b[dep]]+i+ad;
              if(zimu[c[dep]]>=0 && tmp%n!=zimu[c[dep]]) continue;
              if(num[tmp%n] && zimu[c[dep]]!=(tmp%n)) continue;
              if(tmp>=n) jia=1;
              zimu[a[dep]]=i;
              zimu[c[dep]]=tmp%n;
              num[i]=num[tmp%n]=1;
              dfs(dep+1,jia);
              num[i]=num[tmp%n]=0;
              num[flag]=1;zimu[c[dep]]=flag;
          }
        zimu[a[dep]]=-1;
    }
    else
    {
        fo(i,0,n-1) if(!num[i])
        {
            fo(j,0,n-1)
            if((!num[j] && i!=j && a[dep]!=b[dep]) || (!num[j] && i==j && a[dep]==b[dep]))
            {
                jia=0;
                int tmp=i+j+ad;
                if(zimu[c[dep]]>=0 && tmp%n!=zimu[c[dep]]) continue;
                if(num[tmp%n] && zimu[c[dep]]!=(tmp%n)) continue;
                zimu[a[dep]]=i,zimu[b[dep]]=j,zimu[c[dep]]=tmp%n;
                if(tmp>=n) jia=1;
                num[i]=num[j]=num[tmp%n]=1;
                dfs(dep+1,jia);
                num[i]=num[j]=num[tmp%n]=0;
                num[flag]=1;zimu[c[dep]]=flag;
            }
        }
        zimu[a[dep]]=zimu[b[dep]]=-1;
    } 
}
int main()
{
    int i,j;
    scanf("%d",&n);
    scanf("%s%s%s",x+1,y+1,z+1);
    fo(i,1,n)
    {
        a[i]=x[n-i+1];
        b[i]=y[n-i+1];
        c[i]=z[n-i+1]; 
    }
    memset(zimu,-1,sizeof zimu);
    dfs(1,0);
}