NOIP2009 题解 潜伏者 Hankson的趣味题 最优贸易 靶形数独

大白兔的奶糖’s T解

T1 潜伏者

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：通过内线掌握的信息，尝试破译密码。然后利用破译的密码，翻译电报中的加密信息。

【输入输出样例 1】
AA
AB
EOWIE

【输入输出样例 2】
QWERTYUIOPLKJHGFDSAZXCVBN
ABCDEFGHIJKLMNOPQRSTUVWXY
DSLIEWO

【输入输出样例 3】
MSRTZCJKPFLQYVAWBINXUEDGHOOILSMIJFRCOPPQCEUNYDUMPP
YIZSDWAHLNOVFUCERKJXQMGTBPPKOIYKANZWPLLVWMQJFGQYLL
FLSO

【输入输出样例 1】
Failed

【输入输出样例 2】
Failed

【输入输出样例 3】
NOIP

输入

共3 行，每行为一个长度在1 到100 之间的字符串。
第 1 行为小C 掌握的一条加密信息。
第 2 行为第1 行的加密信息所对应的原信息。
第 3 行为R 国司令部要求小C 翻译的加密信息。
输入数据保证所有字符串仅由大写字母‘A’-‘Z’构成，且第1 行长度与第2 行相等。

输出

共1 行。
若破译密码停止时出现 2，3 两种情况，请你输出“Failed”（不含引号，注意首字母大
写，其它小写）。
否则请输出利用密码翻译电报中加密信息后得到的原信息。

题解：
模拟题，然而我最开始只过了两个点，非常遗憾（只得了“Failed”的分），题意应该为第一行的字母对应第二行的字母，我没有考虑到的是两个字母对应一个的情况，此外为了优化空间时间单个字符读入。代码中用了桶排来考虑两个字母对应一个的情况。sec数组为对应情况。

#include "cstdio"
#include "cctype"
#include "cstring"

int inf[101], sec[26], len1, len2, len3;
int bowl[26];
char ch;

int main()  {
    freopen("spy.in", "r", stdin);
    freopen("spy.out", "w", stdout);
    memset(sec, -1, sizeof(sec));
    for( ch = (char)getchar(); isalpha(ch); inf[++len1] = ch - 'A', ch = (char)getchar());
    for( ch = (char)getchar(); isalpha(ch), ++len2;  ch = (char)getchar() )  {
        if( ~sec[inf[len2]] && sec[inf[len2]] != ch - 'A' )  break;
        sec[inf[len2]] = ch - 'A';
    }
    if( len2 <= len1 )  {// <= !
        printf("Failed\n");
        return 0;
    }
    for( int i = 0; i <=25; ++i )  
      if( sec[i] == -1 )  {
        printf("Failed\n");
        return 0;
      }  else  {
        ++bowl[sec[i]];
        if(bowl[i] > 1)  {
            printf("Failed\n");
            return 0;
        }
      }
    for( ch = (char)getchar(); isalpha(ch); inf[++len3] = ch - 'A', ch = (char)getchar());// reuse inf[];
    for( int i = 1; i <= len3; ++i )
      putchar(sec[inf[i]] + 'A');
    putchar('\n');
    return 0;
}

---

T2 Hankson的趣味题

题目描述
　　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 的个数；

样例数据 1
输入

2
41 1 96 288
95 1 37 1776

输出

6
2

【说明】第一组输入数据，x 可以是9、18、36、72、144、288，共有6 个。第二组输入数据，x 可以是48、1776，共有2 个。

备注
【数据范围】
对于 50%的数据，保证有1≤a0，a1，b0，b1≤10000 且n≤100。
对于 100%的数据，保证有1≤a0，a1，b0，b1≤2,000,000,000 且n≤2000。

题解：
使用质数唯一分解优化，用乘法原理计算方案。对于某个质因数x
设a0分解完得到c0个x，a1分解得c1个x，b0分解得c2个x，b1分解得c3个x，ans分解得t个x
那么依照题意，显然c0>=c1,c2<=c3
若c0>c1则t=c1
若c2

#include "cstdio"

int a0, b0, a1, b1, ans, T, prime[10000], mark[50001], tot;
inline void  getprime()  {
    for( register int i = 2; i <= 50000; ++i )  {
      if( !mark[i] )  prime[++tot] = i;
      for (register int j = 1; j<= tot && prime[j] * i <= 50000; ++j)  {
        mark[prime[j]*i] = 1;
        if( prime[j] % i == 0 )  break;
      }
    }
}
void getans( int x )  {
    int g[3];
    g[0] = g[1] = g[2] = g[3] = 0;
    while(!(a0%x))  { a0 /= x; ++g[0]; }
    while(!(a1%x))  { a1 /= x; ++g[1]; }
    while(!(b0%x))  { b0 /= x; ++g[2]; }
    while(!(b1%x))  { b1 /= x; ++g[3]; }
    if( g[0] == g[1] && g[2] == g[3] )  {
        if( g[1] <= g[3])  ans *= g[3] - g[1] + 1;
        else ans = 0;
    }
    else if( g[0] != g[1] && g[2] != g[3] && g[1] != g[3] )  ans = 0;
}

int main()  {
    freopen("son.in", "r", stdin);
    freopen("son.out", "w", stdout);
    getprime();
    scanf("%d", &T);
    while( T-- )  {
        ans = 1;
        scanf("%d%d%d%d", &a0, &a1, &b0, &b1);
        for( register int i = 1; i <= tot; ++i )  
            getans(prime[i]);
        if( b1 != 1 )  getans(b1);
        printf("%d\n",ans);
    }
    return 0;
}

T3 最优贸易

题目描述 Description

【问题描述】
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 条道路的信息（每条道路所连接的两个城市的编号
以及该条道路的通行情况）。请你告诉阿龙，他最多能赚取多少旅费。

输入描述 Input Description

第一行包含 2 个正整数n 和m，中间用一个空格隔开，分别表示城市的数目和道路的
数目。
第二行 n 个正整数，每两个整数之间用一个空格隔开，按标号顺序分别表示这n 个城
市的商品价格。
接下来 m 行，每行有3 个正整数，x，y，z，每两个整数之间用一个空格隔开。如果z=1，
表示这条道路是城市x 到城市y 之间的单向道路；如果z=2，表示这条道路为城市x 和城市
y 之间的双向道路。

输出描述 Output Description

包含1 个整数，表示最多能赚取的旅费。如果没有进行贸易，
则输出0。

样例输入 Sample Input

5 5
4 3 5 6 1
1 2 1
1 4 1
2 3 2
3 5 1
4 5 2

样例输出 Sample Output

5

数据范围及提示 Data Size & Hint

【数据范围】
输入数据保证 1 号城市可以到达n 号城市。
对于 10%的数据，1≤n≤6。
对于 30%的数据，1≤n≤100。
对于 50%的数据，不存在一条旅游路线，可以从一个城市出发，再回到这个城市。
对于 100%的数据，1≤n≤100000，1≤m≤500000，1≤x，y≤n，1≤z≤2，1≤各城市水晶球价格≤100。

题解：听他们说用SPFA，可是我觉得明显是Tarjan加dfs啊。先用Tarjan求强连通分量缩点，Tarjan时顺便求出该强连通分量里面的最大最小费用（用bmx[]与bmn[]求）。然后重新建图，之后dfs就好。比一般的双向SPFA快一点（0.2s-0.3s左右）

#include "cstdio"
#include "cctype"
#include "cstring"
#include "algorithm"
using namespace std;

template <class T>
inline bool readIn( T &x )  {
    x = 0;
    char ch = (char) getchar();
    while ( !isdigit(ch) )  ch = (char) getchar();
    while ( isdigit(ch) )  {
        x = ( x << 1 ) + ( x << 3 ) + ch - 48;
        ch = (char) getchar();
    }
}

const int MAXN = 100005, MAXM = 500005;

int a[MAXN], place[MAXN], head[MAXN], dfn[MAXN], low[MAXN], stack[MAXN], u[MAXN], v[MAXN], top, cnt, n, m, z, ne, idc;
int bmx[MAXN], bmn[MAXN], f[MAXN], ans = -1;
bool vis[MAXN];  


struct edge  {
    int to, pre;
    edge()  {
        to = pre = 0;
    }
    inline bool adde( int to, int pre )  {
        this -> to = to, this -> pre = pre;
    }
}g[MAXM];

inline bool adde( int u, int v )  {
    g[++ne].adde(v, head[u]), head[u] = ne;
}

void Tarjan( int u )  {
    dfn[u] = low[u] = ++idc;  vis[u] = true;  stack[++top] = u;
    for( int i = head[u]; i; i = g[i].pre )  {
        int v = g[i].to;
        if( !dfn[v] )  {
            Tarjan(v);
            low[u] = min(low[u], low[v]);
        }  else  if( vis[v] )
           low[u] = min( dfn[v], low[u] );
    }

    if( low[u] == dfn[u]  )  {
        vis[u] = false;
        place[u] = ++cnt;
        bmx[cnt] = max(bmx[cnt], a[u]);  
        bmn[cnt] = min(bmn[cnt], a[u]); 
        while( stack[top] != u )  {
            place[stack[top]] = cnt;
            bmx[cnt] = max(bmx[cnt], a[stack[top]]);  
            bmn[cnt] = min(bmn[cnt], a[stack[top]]); 
            vis[stack[top--]] = false;
        }
        --top;
    }
}

void dfs(int u) {  
    vis[u] = true;  
    if (u == place[n]) f[u] = max(f[u], bmx[u]);  
    for (int i = head[u]; i; i = g[i].pre) {  
        int v = g[i].to;
        if (!vis[v]) dfs(v);  
        f[u] = max(f[u], f[v]);  
    }  
    if (f[u]) f[u] = max(f[u], bmx[u]); 
    ans = max(ans, f[u] - bmn[u]);
} 

int main()  {
    freopen("trade.in", "r", stdin);
    freopen("trade.out", "w", stdout);
    readIn(n), readIn(m);
    for( register int i = 1; i <= n; ++i )  readIn(a[i]);
    for( register int i = 1; i <= m; ++i )  {
        readIn(u[i]), readIn(v[i]), readIn(z);
        adde( u[i], v[i] );  if(z == 2) adde(v[i], u[i]);
    }
    memset(bmn, 0x3f, sizeof(bmn));
    Tarjan(1);
    ne = 0;
    memset(head, 0, sizeof(head));
    for (register int i = 1; i <= m; ++i)  
        if (place[u[i]] != place[v[i]])  
            adde(place[u[i]], place[v[i]]);  
    dfs(place[1]);  
    printf("%d\n", ans);  
    return 0;  
}

T4 靶形数独

题目描述 Description

小城和小华都是热爱数学的好学生，最近，他们不约而同地迷上了数独游戏，好胜的他
们想用数独来一比高低。但普通的数独对他们来说都过于简单了，于是他们向Z 博士请教，
Z 博士拿出了他最近发明的“靶形数独”，作为这两个孩子比试的题目。
靶形数独的方格同普通数独一样，在 9 格宽×9 格高的大九宫格中有9 个3 格宽×3 格
高的小九宫格（用粗黑色线隔开的）。在这个大九宫格中，有一些数字是已知的，根据这些

数字，利用逻辑推理，在其他的空格上填入1 到9 的数字。每个数字在每个小九宫格内不能
重复出现，每个数字在每行、每列也不能重复出现。但靶形数独有一点和普通数独不同，即
每一个方格都有一个分值，而且如同一个靶子一样，离中心越近则分值越高。

上图具体的分值分布是：最里面一格（黄色区域）为 10 分，黄色区域外面的一圈（红
色区域）每个格子为9 分，再外面一圈（蓝色区域）每个格子为8 分，蓝色区域外面一圈（棕
色区域）每个格子为7 分，最外面一圈（白色区域）每个格子为6 分，如上图所示。比赛的
要求是：每个人必须完成一个给定的数独（每个给定数独可能有不同的填法），而且要争取
更高的总分数。而这个总分数即每个方格上的分值和完成这个数独时填在相应格上的数字
的乘积的总和。如图，在以下的这个已经填完数字的靶形数独游戏中，总分数为2829。游
戏规定，将以总分数的高低决出胜负。

由于求胜心切，小城找到了善于编程的你，让你帮他求出，对于给定的靶形数独，能
够得到的最高分数。

输入描述 Input Description

一共 9 行。每行9 个整数（每个数都在0—9 的范围内），表示一个尚未填满的数独方
格，未填的空格用“0”表示。每两个数字之间用一个空格隔开。

输出描述 Output Description

输出可以得到的靶形数独的最高分数。如果这个数独无解，则输出整数-1。

样例输入 Sample Input

【输入输出样例 1】

7 0 0 9 0 0 0 0 1
1 0 0 0 0 5 9 0 0
0 0 0 2 0 0 0 8 0
0 0 5 0 2 0 0 0 3
0 0 0 0 0 0 6 4 8
4 1 3 0 0 0 0 0 0
0 0 7 0 0 2 0 9 0
2 0 1 0 6 0 8 0 4
0 8 0 5 0 4 0 1 2

【输入输出样例 2】

0 0 0 7 0 2 4 5 3
9 0 0 0 0 8 0 0 0
7 4 0 0 0 5 0 1 0
1 9 5 0 8 0 0 0 0
0 7 0 0 0 0 0 2 5
0 3 0 5 7 9 1 0 8
0 0 0 6 0 1 0 0 0
0 6 0 9 0 0 0 0 1
0 0 0 0 0 0 0 0 6

样例输出 Sample Output

【输入输出样例 1】

2829

【输入输出样例 1】

2852

数据范围及提示 Data Size & Hint

【数据范围】
40%的数据，数独中非0 数的个数不少于30。
80%的数据，数独中非0 数的个数不少于26。
100%的数据，数独中非0 数的个数不少于24。

题解：

纯粹就是大暴力嘛，注意如果输入的时候反向要快得多，哈哈哈。注意要卡时间，不然倒数第二个点过不了。

#include "ctime"
#include "cstdio"
#include "cstdlib"
#define max(a, b) (a) > (b) ? (a) : (b)

bool r[10][10], l[10][10], s[10][10];
int a[10][10], b[10][10];
int getscore[10][10][10];
int ans = -1, score, t;

int culscore(int x,int y,int k)  {
    if(x==5&&y==5)return 10*k;
    else if(x>=4&&x<=6&&y>=4&&y<=6)return 9*k;
    else if(x>=3&&x<=7&&y>=3&&y<=7)return 8*k;
    else if(x>=2&&x<=8&&y>=2&&y<=8)return 7*k;
    else return 6*k;
}

inline bool fillin(int x,int y,int k)  {
    if(r[x][k]||l[y][k]||s[(x-1)/3*3+(y-1)/3][k])return false;
    b[x][y]=k;
    r[x][k]=l[y][k]=s[(x-1)/3*3+(y-1)/3][k]=1;
    score+=getscore[x][y][k];
    return true;
}

void search(int x,int y)  {
    if( (double) (clock() - t) / CLOCKS_PER_SEC > 1.95)  {
        printf("%d\n",ans);
        exit(0);
    }
    if(x == 10 && y == 1)  {
        ans = max(ans, score);
        return;
    }
    if(b[x][y]) y == 9? search(x + 1, 1) : search(x, y + 1);
    else   {
      for(int i = 1; i <= 9; ++i)  {
          int t = score;
          if(fillin(x, y, i))  {
              y == 9? search(x + 1, 1) : search(x, y + 1);
              r[x][i] = l[y][i] = s[(x-1)/3*3+(y-1)/3][i] = 0;
              score = t;
          }
      }
      b[x][y] = 0;
    } 
}

int main()  {
    t = clock();
    freopen("sudoku.in", "r", stdin);
    freopen("sudoku.out", "w", stdout);
    for(register int i = 1; i <= 9; ++i)
       for(register int j = 1; j <= 9; ++j)
          for(register int k = 1; k <= 9; ++k)
             getscore[i][j][k] = culscore(i,j,k);
    for(int i=9;i>0;i--)
       for(int j=9;j>0;j--)  {
          scanf("%d",&a[i][j]);
          if(a[i][j])  fillin(i, j, a[i][j]);
      }
    search(1, 1);
    printf("%d\n",ans);
    return 0;
}

大家都很强，可与之共勉