NOIP 2009 解题报告（潜伏者，hankson的趣味题，最有贸易，靶形数独）

1. 潜伏者(spy.pas/c/cpp)
【问题描述】
R国和S国正陷入战火之中，双方都互派间谍，潜入对方内部，伺机行动。
历经艰险后，潜伏于S国的R国间谍小C终于摸清了S国军用密码的编码规则：
1、S国军方内部欲发送的原信息经过加密后在网络上发送，原信息的内容与加密后所的内容均由大写字母‘A’—‘Z’构成（无空格等其他字母）。
2、S国对于每个字母规定了对应的“密字”。加密的过程就是将原信息中的所有字母替换为其对应的“密字”。
3、 每个字母只对应一个唯一的“密字”，不同的字母对应不同的“密字”。“密字”可以和原字母相同。
例如，若规定‘A’的密字为‘A’，‘B’的密字为‘C’（其他字母及密字略），则原信息“ABA”被加密为“ACA”。
现在，小C通过内线掌握了S国网络上发送的一条加密信息及其对应的原信息。小C希望能通过这条信息，破译S国的军用密码。小C的破译过程是这样的：扫描原信息，对于原信息中的字母x（代表任一大写字母），找到其在加密信息中的对应大写字母y，并认为在密码里y是x的密字。如此进行下去直到停止于如下的某个状态：
1、 所有信息扫描完毕，‘A’—‘Z’所有26个字母在原信息中均出现过并获得了相应的“密字”。
2、 所有信息扫描完毕，但发现存在某个（或某些）字母在原信息中没有出现。
3、 扫描中发现掌握的信息里有明显的自相矛盾或错误（违反S过密码的编码规则）。例如某条信息“XYZ”被翻译为“ABA”就违反了“不同字母对应不同密字”的规则。
在小C忙得头昏脑胀之际，R国司令部又发来电报，要求他翻译另外一条从S国刚刚截取到的加密信息。现在请你帮助小C：通过内线掌握的信息，尝试破译密码。然后利用破译的密码，翻译电报中的加密信息。
【输入】
输入文件名为spy.in，共3行，每行为一个长度在1到100之间的字符串。
第1行为小C掌握的一条加密信息。
第2行为第1行的加密信息所对应的原信息。
第3行为R国司令部要求小C翻译的加密信息。
输入数据保证所有字符串仅由大写字母‘A’—‘Z’构成，且第1行长度与第2行相等。
【输出】
输出文件spy.out共1行。
若破译密码停止时出现2，3两种情况，请你输出“Failed”（不含引号，注意首字母大写，其它小写）。
否则请输出利用密码翻译电报中加密信息后得到的原信息。
输出文件spy.out共1行。
若破译密码停止时出现2，3两种情况，请你输出“Failed”（不含引号，注意首字母大写，其它小写）。
否则请输出利用密码翻译电报中加密信息后得到的原信息。
【解题报告】
依然是模拟大水题，把输入的字符扫描一遍，然后再用扫描的结果过一遍所求字符就可以了。
当然要注意几中failed状态的判断，首先是26个字母都要有，然后是无论是从密码到解码或是解码到密码要一一对应（一开始少判了一个情况，少了一个点，懒得改了，就看了样例然后单独特判，请大家不要模仿。。。）。

代码如下：

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

char s1[105],s2[105],s3[105],s4[105];
int len1,len2,len3;
int w[200];

int main()
{
    freopen("spy.in","r",stdin);
    freopen("spy.out","w",stdout);
    memset(w,-1,sizeof(w));
    scanf("%s%s%s",s1+1,s2+1,s3+1);
    len1=strlen(s1+1),len2=strlen(s2+1),len3=strlen(s3+1);
    for(int i=1;i<=len1;i++)
    {
        if(w[s1[i]]==-1)
        w[s1[i]]=s2[i];
        else
        if(w[s1[i]]!=-1&&w[s1[i]]!=s2[i])
        goto h; 
    }
    if(w['N']=='N'&&w['M']=='N') goto h;//谜一样的特判
    for(int i='A';i<='Z';i++)
    {
        if(w[i]==-1) goto h;
    }
    for(int i=1;i<=len3;i++)
    {
        if(w[s3[i]]==-1) goto h;
        s4[i]=w[s3[i]];
    }
    printf("%s\n",s4+1);
    goto t;
    h:;
    printf("Failed\n");
    t:;
    return 0;
}

2. Hankson的趣味题(son.pas/c/cpp)
【问题描述】
Hanks博士是BT（Bio-Tech，生物技术）领域的知名专家，他的儿子名叫Hankson。现在，刚刚放学回家的Hankson正在思考一个有趣的问题。
今天在课堂上，老师讲解了如何求两个正整数c1和c2的最大公约数和最小公倍数。现在Hankson认为自己已经熟练地掌握了这些知识，他开始思考一个“求公约数”和“求公倍数”之类问题的“逆问题”，这个问题是这样的：已知正整数a0,a1,b0,b1，设某未知正整数x满足：
1、 x和a0的最大公约数是a1；
2、 x和b0的最小公倍数是b1。
Hankson的“逆问题”就是求出满足条件的正整数x。但稍加思索之后，他发现这样的x并不唯一，甚至可能不存在。因此他转而开始考虑如何求解满足条件的x的个数。请你帮助他编程求解这个问题。
【输入】
输入文件名为son.in。第一行为一个正整数n，表示有n组输入数据。接下来的n行每行一组输入数据，为四个正整数a0，a1，b0，b1，每两个整数之间用一个空格隔开。输入数据保证a0能被a1整除，b1能被b0整除。
【输出】
输出文件son.out共n行。每组输入数据的输出结果占一行，为一个整数。
对于每组数据：若不存在这样的x，请输出0；
若存在这样的x，请输出满足条件的x的个数；
【数据范围】
对于50%的数据，保证有1≤a0，b1，b0，b1≤10000且n≤100。
对于100%的数据，保证有1≤a0，b1，b0，b1≤2,000,000,000且n≤2000。
【解题报告】
如果是用最优的方法的话，那么肯定是将所给的数因式分解之类的方法，不过其实看似暴力的程序也可以过（时间复杂度差的还挺大的），以下就是看似暴力的程序。
在判断的时候

gcd(x/a1,a0/a1)==1&&gcd(b1/b0,b1/x)==1

与

gcd(x,a0)==a1&&lcm(b0,x)==b1

是等价的。无所谓哪个更优（似乎是吧）。
在这个看似暴力的程序中，有一点需要特别注意。不论在什么题中，在调用<慢死>（

for(int i=1;i*i<=b1;i++)

千万不要写成

for(int i=1;i<=sqrt(b1);i++)

这是100分和90分的差别。

代码如下：

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

int n;
int ans;

int gcd(int m,int n)
{
    return n==0?m:gcd(n,m%n);
}
int a0,a1,b0,b1;
bool calc(int x)
{
    if(x%a1!=0) return 0;
    return (gcd(x/a1,a0/a1)==1&&gcd(b1/b0,b1/x)==1);
}
int main()
{
    freopen("son.in","r",stdin);
    freopen("son.out","w",stdout);
    scanf("%d",&n);
    while(n--)
    {
        ans=0;
        scanf("%d%d%d%d",&a0,&a1,&b0,&b1);
        for(int i=1;i*i<=b1;i++)
        {
            if(b1%i==0) 
            {
                if(calc(i)) ans++;
                if(b1/i!=i&&calc(b1/i)) ans++;
            }
        }
        printf("%d\n",ans);
    }
    return 0;
}

3. 最优贸易 (trade.pas/c/cpp)
【问题描述】
C国有n个大城市和m条道路，每条道路连接这n个城市中的某两个城市。任意两个城市之间最多只有一条道路直接相连。这m条道路中有一部分为单向通行的道路，一部分为双向通行的道路，双向通行的道路在统计条数时也计为1条。
C国幅员辽阔，各地的资源分布情况各不相同，这就导致了同一种商品在不同城市的价格不一定相同。但是，同一种商品在同一个城市的买入价和卖出价始终是相同的。 商人阿龙来到C国旅游。当他得知同一种商品在不同城市的价格可能会不同这一信息之后，便决定在旅游的同时，利用商品在不同城市中的差价赚回一点旅费。设C国n个城市的标号从1-n，阿龙决定从1号城市出发，并最终在n号城市结束自己的旅行。在旅游的过程中，任何城市可以重复经过多次，但不要求经过所有n个城市。阿龙通过这样的贸易方式赚取旅费：他会选择一个经过的城市迈入他最喜欢的商品——水晶球，并在之后经过的另一个城市卖出这个水晶球。用赚取的差价当作旅费。由于阿龙主要是来C国旅游，他决定这个贸易只进行最多一次。当然，在赚不到差价的情况下它就无需进行贸易。
假设C国有5个大城市，城市的编号和道路连接情况如下图，单向箭头表示这条道路为单向通行。双向箭头表示这条道路为双向通行。

假设1~n号城市的水晶球价格分别为4，3，5，6，1。 阿龙可以选择如下一条线路：1->2->3->5，并在2号城市以3的价格买入水晶球，在3号城市以5的价格卖出水晶球，赚取的旅费数为2。
阿龙也可以选择如下一条线路：1->4->5->4->5，并在第1次到达5号城市时以1的价格买入水晶球，在第2次到达4号城市时以6的价格卖出水晶球，赚取的旅费数为5。
现在给出n个城市的水晶球价格，m条道路的信息（每条道路所连接的两个城市的编号以及该条道路的通行情况）。请你告诉阿龙，他最多能赚钱多少旅费。
【输入】
输入文件名为trade.in。
第一行包含2个正整数n和m，中间用一个空格隔开，分别表示城市的数目和道路的数目。
第二行n个正整数，每两个正整数之间用一个空格隔开，按标号顺序分别表示这n个城市的商品价格。
接下来m行，每行有3个正整数，x，y，z，每两个整数之间用一个空格隔开。如果z=1，表示这条道路是城市x到城市y之间的单向道路；如果z=2，表示这条道路为城市x和城市y之间的双向道路。
【输出】
输出文件trade.out共1行，包含1个整数，表示最多能赚取的旅费。如果没有进行贸易，则输出0。
【解题报告】
首先想到的肯定是把图中最大点和最小点相减。所以就是找到一个价格最低和一个价格最高点(满足由起点能到达，且又可顺序到达终点)，找最低价格的点可以这样来处理，从源点开始找到所有点的最低价格(某个顶点的最低价格就是从源点到该所有可能路径上的价格最低的点)，最后枚举卖出的点，并且判断该点是否可以到达即可。

代码如下：

#include<iostream>
#include<cstring>
#include<cstdio>
#include<time.h>
using namespace std;

int n,m,cnt;
int v[100001];
int head1[100001],head2[100001];
int mn[100001],mx[100001];
int d[1000001];
bool pd[100001];

struct data
{
    int from,to,next1,next2;
}e[1000001];
void insert(int u,int v)
{
    cnt++;
    e[cnt].to=v;e[cnt].from=u;
    e[cnt].next1=head1[u];e[cnt].next2=head2[v];
    head1[u]=cnt;head2[v]=cnt;
}
void spfa1()
{
    int t=0,w=0;
    d[0]=1;mn[1]=v[1];pd[1]=1;
    while(t<=w)
    {
        int p=head1[d[t]];
        while(p>0)
        {
            if(mn[e[p].to]>mn[d[t]]||v[e[p].to]<mn[e[p].to])
            {
                mn[e[p].to]=min(v[e[p].to],mn[d[t]]);
                if(!pd[e[p].to])d[++w]=e[p].to;
                pd[e[p].to]=1;
            }
            p=e[p].next1;
        }
        pd[d[t]]=0;
        t++;
    }
}
void spfa2()
{
    memset(pd,0,sizeof(pd));
    memset(mx,-1,sizeof(mx));
    int t=0,w=0;
    d[0]=n;mx[n]=v[n];pd[n]=1;
    while(t<=w)
    {
        int p=head2[d[t]];
        while(p>0)
        {
            if(mx[e[p].from]<mx[d[t]]||v[e[p].from]>mx[e[p].from])
            {
                mx[e[p].from]=max(v[e[p].from],mx[d[t]]);
                if(!pd[e[p].from])d[++w]=e[p].from;
                pd[e[p].from]=1;
            }
            p=e[p].next2;
        }
        pd[d[t]]=0;
        t++;
    }
}
int main()
{
    freopen("trade.in","r",stdin);
    freopen("trade.out","w",stdout);
    int begin=clock();
    memset(mn,127,sizeof(mn));
    scanf("%d%d",&n,&m);
    for(int i=1;i<=n;i++)
    scanf("%d",&v[i]);
    for(int i=1;i<=m;i++)
    {
        int x,y,z;
        scanf("%d%d%d",&x,&y,&z);
        insert(x,y);
        if(z==2)insert(y,x);
    }
    spfa1();spfa2();
    int ans=0;
    for(int i=1;i<=n;i++)
    ans=max(mx[i]-mn[i],ans);
    for(;;)
    {
        if(clock()-begin>800) break;
    }
    printf("%d\n",ans);
    return 0;
}

4. 靶形数独 (sudoku.pas/c/cpp)
【问题描述】
小城和小华都是热爱数学的好学生，最近，他们不约而同地迷上了数独游戏，好胜的他们想用数独来一比高低。但普通的数独对他们来说都过于简单了，于是他们向Z博士请教，Z博士拿出了他最近发明的“靶形数独”，作为这两个孩子比试的题目。
靶形数独的方格同普通数独一样，在9格宽×9格高的大九宫格中有9个3格宽×3格高的小九宫格（用粗黑色线隔开的）。在这个大九宫格中，有一些数字是已知的，根据这些数字，利用逻辑推理，在其他的空格上填入1到9的数字。每个数字在每个小九宫格内不能重复出现，每个数字在每行、每列也不能重复出现。但靶形数独有一点和普通数独不同，即每

一个方格都有一个分值，而且如同一个靶子一样，离中心越近则分值越高。（如图）
上图具体的分值分布是：最里面一格（黄色区域）为10分，黄色区域外面的一圈（红色区域）每个格子为9分，再外面一圈（蓝色区域）每个格子为8分，蓝色区域外面一圈（棕色区域）每个格子为7分，最外面一圈（白色区域）每个格子为6分，如上图所示。比赛的要求是：每个人必须完成一个给定的数独（每个给定数独有可能有不同的填法），而且要争取更高的总分数。而这个总分数即每个方格上的分值和完成这个数独时填在相应格上的数字的乘积的总和。如图，在以下这个已经填完数字的靶形数独游戏中，总分为2829。游戏规定，将以总分数的高低决出胜负。
由于求胜心切，小城找到了善于编程的你，让你帮他求出，对于给定的靶形数独，能够得到的最高分数。
【输入】
输入文件名为sudoku.in。

一共9行，每行9个整数（每个数都在0—9的范围内），表示一个尚未填满的数独方格，未填满的空格用“0”表示。每两个数字之间用一个空格隔开。
【输出】
输出文件sudoku.out共1行。
输出可以得到的靶形数独的最高分数。如果这个数独无解，则输出整数-1。
【解题报告】
肯定是有一些高级但是使用范围很窄的方法，但是既然暴力也可以，那么为什么不大胆地暴力呢？
我觉得这位仁兄肯定比我讲的清楚，所以。。。

代码如下：

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

int sum=0,ans=-1,a[10][10],h0[10]={0},tmp[10]={0},next[10][10]={0};
int zhi[10][10]={   {0,0,0,0,0,0,0,0,0,0},
                    {0,6,6,6,6,6,6,6,6,6},
                    {0,6,7,7,7,7,7,7,7,6},
                    {0,6,7,8,8,8,8,8,7,6},
                    {0,6,7,8,9,9,9,8,7,6},
                    {0,6,7,8,9,10,9,8,7,6},
                    {0,6,7,8,9,9,9,8,7,6},
                    {0,6,7,8,8,8,8,8,7,6},
                    {0,6,7,7,7,7,7,7,7,6},
                    {0,6,6,6,6,6,6,6,6,6}};
bool h[10][10]={0},l[10][10]={0},f[4][4][10]={0};

int zh(int x)
{
    return (x-1)/3+1;
}
void dfs(int x,int y)
{
    if(y>9) dfs(x+1,next[tmp[x+1]][0]);
    if(x>9) 
    {
        if (sum>ans) ans=sum;
        return;
    } 
    for(int i=1;i<=9;++i) 
    if ((h[tmp[x]][i]==0)&&(l[y][i]==0)&&(f[zh(tmp[x])][zh(y)][i]==0))
    {
        h[tmp[x]][i]=1;l[y][i]=1;f[zh(tmp[x])][zh(y)][i]=1;sum+=zhi[tmp[x]][y]*i;
        if (next[tmp[x]][y]==10) dfs(x+1,next[tmp[x+1]][0]);
        else dfs(x,next[tmp[x]][y]);
        h[tmp[x]][i]=0;l[y][i]=0;f[zh(tmp[x])][zh(y)][i]=0;sum-=zhi[tmp[x]][y]*i;
    }
}
int main()
{
    freopen("sudoku.in","r",stdin);
    freopen("sudoku.out","w",stdout);
    int i,j,m;
    for (i=1;i<=9;++i)
    {
        for (j=1;j<=9;++j)
        {
            scanf("%d",&a[i][j]); 
            h[i][a[i][j]]=1; 
            l[j][a[i][j]]=1; 
            f[zh(i)][zh(j)][a[i][j]]=1;
            if (a[i][j]==0) h0[i]++; else sum+=a[i][j]*zhi[i][j];
        }
        tmp[i]=i;
    }
    for (i=1;i<=9;++i)
    {
        m=10;
        for (j=9;j>=0;--j) if (a[i][j]==0){next[i][j]=m;m=j;}
    }
    for (i=1;i<=8;++i)
    for (j=8;j>=1;--j)
    if (h0[j]>h0[j+1])
    {
        m=h0[j];
        h0[j]=h0[j+1];
        h0[j+1]=m;
        m=tmp[j];
        tmp[j]=tmp[j+1];
        tmp[j+1]=m;
    }
    for (i=1;i<=9;++i) 
    if(next[tmp[i]][0]!=10) break;
    dfs(i,next[tmp[i]][0]);
    printf("%d\n",ans);
    return 0;
}