BZOJ2570：算符破译题解

转载自Array98大大的博客

【题目描述】 
考古学发现，几千年前古梅文明时期的数学非常的发达，他们懂得多位数的加法和乘法，其表达式和运算规则等都与现在通常所用的方式完全相同（如整数是十进制，左边是高位，最高位不能为零；表达式为中缀运算，先乘后加等），唯一的区别是其符号的写法与现在不同。有充分的证据表明，古梅文明的数学文字一共有13个符号，与 0,1,2,3,4,5,6,7,8,9,+,*,= 这13个数字和符号（称为现代算符）一一对应。为了便于标记，我们用13个小写英文字母a,b,…m代替这些符号（称为古梅算符）。但是，还没有人知道这些古梅算符和现代算符之间的具体对应关系。 
在一个石壁上，考古学家发现了一组用古梅算符表示的等式，根据推断，每行有且仅有一个等号，等号左右两边为运算表达式（只含有数字和符号），并且等号两边的计算结果相等。 
假设这组等式是成立的，请编程序破译古梅算符和现代算符之间的对应关系。

【输入格式】 
输入文件的第一行为等式的个数N（1<=N<=1000），以下N行每行为一个等式。 
每个等式的长度为5个字符到11个字符。

【输出格式】 
如果不存在对应关系能够满足这组等式，输出“noway”和一个换行/回车符。 
如果有对应关系能够满足这组等式，输出所有能够确定的古梅算符和现代算符的对应关系。每一行有两个字符，其中第一个字符是古梅算符，第二个字符是对应的现代算符。输出按照字典顺序排序。

【样例输入】 
2 
abcdec 
cdefe

【样例输出】 
a6 
b* 
d= 
f+

【样例说明】 
在上例中，可能对应的现代表达式为{6*2=12，2=1+1}，{6*4=24，4=2+2}，{6*8=48，8=4+4}。可见，能够确定的对应关系只有a对应6，b对应*，d对应=，f对应+，应该输出；而{c,e}虽然能够找到对应的现代算符使得等式成立，但没有唯一的对应关系，不能输出。其他古梅算符{g,h…m}完全不能确定，也不能输出。

【分析】 
相比消棋子我觉得还是这道题酸爽一点……毕竟不是单纯的模拟还要敲搜索疯狂的剪枝…………总之我讨厌写搜索，最让人痛恨的是调试搜索TAT 
没有充足心里准备不要敲这道题= =

反正终归是水过了嘛……下面来讲讲怎么做的。 
首先，裸的搜索要敲对吧（PS：爆搜连样例都过不去调试毛线啊……） 
这里无视Std各种各样的优化（标程一共敲了600行）

然后编写是有一点要注意我们可以优先确定三个操作符的对应关系（体现在代码中的三个for循环，这样搜索就只需要搜数字）

第一个优化：搜索时位运算压位，方便编写，加快速度，这个看代码应该都能理解。 
第二个优化：其实算不上优化，就是我开始比较蠢……没写位运算，然后就是说把判断放到了搜索完之后……显而易见 O(N!) 连样例都过不了。 
第三个优化：确定了运算符之后我们就可以大致确定出等号两边的数的位数，也就是大概的范围，这样可以直接跳过大量不可能的情况。最重要的优化。（对应代码中的Pre_Dfs）

/**************************************************************
    Problem: 2570
    User: Array98
    Language: C++
    Result: Accepted
    Time:136 ms
    Memory:1060 kb
****************************************************************/

#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;

const int INF=int(1e9);
const int MaxN=1010;
int T;
char a[MaxN][15]; 
int len[MaxN],lg2[(1<<10)+10];
int w[15],Ans[15],st[15],cnt[15][MaxN];
bool c1[15],c2[15],mark[15],Map[150][1150];
char Equ,Mul,Add;
bool Solved;

inline int Id(char x) { return x-'a'; }

int calc()
{
    int i,sum=0;
    for(;mark[st[0]-1];st[0]--)
        st[st[0]-1]=st[st[0]-1]*st[st[0]];
    for (int i=1; i<=st[0]; i++) sum+=st[i]; st[0]=0;
    return sum;
}

bool check(int d)
{
    int Last=-1,L=-1,R=-1; st[0]=0;
    for (int i=1; i<=len[d]; i++)
    {
        int c=w[Id(a[d][i])];
        if(c>=0)
            if(Last>=0)st[st[0]]=st[st[0]]*10+c;
                else st[++st[0]]=c;
        else
        {
            if(st[0]==0) return 0;
            for(;mark[st[0]-1];st[0]--)
                st[st[0]-1]=st[st[0]-1]*st[st[0]];

            if(c==-3) mark[st[0]]=1;
                else if(c==-1)L=calc();
                    else mark[st[0]]=0;
        }
        Last=c;
    }
    R=calc();
    return L==R;
}

void Dfs(int d,int Now,int S)
{
    bool Flag=1;
    for (int i=0; i<=12; i++) 
        if(!(Ans[i]==w[i] || Ans[i]==-200)){ Flag=0; break; }
    if(Flag) return;
    if(d>T)
    {   
        Solved=1;
        for (int i=0; i<=12; i++) if(w[i]!=INF)
        {
              if(Ans[i]==-100) Ans[i]=w[i];
              else if(Ans[i]!=w[i])Ans[i]=-200;
        }
        else Ans[i]=-200;
        return;
    }
    else
    {
        if(Now>len[d])
        {
            if(check(d)) Dfs(d+1,1,S);
            return;
        }
        int ch=Id(a[d][Now]);
        if(w[ch]==INF)
        {
            for(int S2=S; S2; S2-=S2&-S2)
            {
                int t=lg2[S2&-S2];
                if(Now==1 && t==0 && w[Id(a[d][2])]>0)continue;
                w[ch]=t;
                Dfs(d,Now+1,S-(S2&-S2));
                w[ch]=INF;
            }
        }
        else Dfs(d,Now+1,S);
    }
}

int M,lw[15],rw[15];

void Get_Range(int &L, int &R)
{ 
    L=R=0; 
    for (;M;M--) L=max(L,lw[M]),R=max(rw[M],R)+(R>0); 
}

bool Pre_Dfs()
{
    int Last,L_l,L_r,R_l,R_r;
    for (int d=1; d<=T; d++)
    {
        M=0;Last=-1;
        for (int i=1; i<=len[d]; i++)
        {
            int c=w[Id(a[d][i])];
            if(c==INF)
            { 
                if(Last<0) M++,lw[M]=rw[M]=0; 
                lw[M]++,rw[M]++;
            }
            else
            {
                if(M==0)return 0;
                for(; mark[M-1]; M--) lw[M-1]+=lw[M]-1,rw[M-1]+=rw[M];

                     if(c==-3) mark[M]=1;
                else if(c==-2) mark[M]=0;
                else if(c==-1) Get_Range(L_l,L_r);
            }
            Last=c;
        }
        Get_Range(R_l,R_r);
        if(L_r<R_l || R_r<L_l)
        {
            return 0;
        }
    }
    return 1;
}

int main()
{
    scanf("%d",&T);
    for (int i=1; i<=T; i++) scanf("%s",a[i]+1),len[i]=strlen(a[i]+1);
    for (int i=0; i<=15; i++) lg2[1<<i]=i,Ans[i]=-100;     
    memset(c1,1,sizeof(c1));

    for(char i='a';i<='m';i++)
        for (int j=1; j<=T; j++)
        {
            bool Flag=1;
            for (int k=1; k<=len[j]; k++)
            {
                if(k<len[j]) Map[a[j][k]][a[j][k+1]]=Map[a[j][k+1]][a[j][k]]=1;
                if(a[j][k]==i) cnt[Id(i)][j]++;
            }
            if(cnt[Id(i)][j]!=1) Flag=0;
            if(a[j][1]==i || a[j][len[j]]==i){ c1[i-'a']=0; Flag=0; }
            c2[Id(i)]=Flag;
        }

    for(Equ='a'; Equ<='m'; Equ++) if(c2[Id(Equ)])
        for(Mul='a'; Mul<='m'; Mul++) if(c1[Id(Mul)])
            for(Add='a'; Add<='m'; Add++) if(c1[Id(Add)])
                if(Equ!=Mul && Add!=Mul && Equ!=Add)
                    if(!Map[Equ][Mul] && !Map[Equ][Add] && !Map[Mul][Add])
                    {                 
                        for (int i=0; i<=12; i++) w[i]=INF;
                        w[Id(Equ)]=-1; w[Id(Add)]=-2; w[Id(Mul)]=-3;
                        if(!Pre_Dfs()) continue;
                        Dfs(1,1,(1<<10)-1);
                    }

    for (int i=0; i<=12; i++)
        if(Ans[i]>=-3)
        {
            putchar('a'+i);
                 if(Ans[i]==-3)putchar('*');
            else if(Ans[i]==-2)putchar('+');
            else if(Ans[i]==-1)putchar('=');
            else putchar('0'+Ans[i]);
            putchar('\n');
        }
    if(!Solved) printf("noway\n");

    return 0;
}

---

我自己打了5h，打了6k+，
最后WA了，是在下输了
蒟蒻就是蒟蒻