poj1900 Game

Posted 2020-07-19 lyxxx

http://poj.org/problem?id=1900

一道有意思的很好的题目。

A，S，P三个人玩一个游戏。A先宣布一个正整数N(2<=N<=200)，然后选择两个不同的正整数x，y(1<=x,y<=N)，然后把x+y的值告诉S，把x*y的值告诉P。

接下来S和P轮流告诉对方自己现在是否知道这两个数是多少，S先开始。例如：

S和P都已知N=10，S知道x+y=9，P知道x*y=18。

S：我不知道x和y的值。

P：我也不知道。

S：我还不知道。

P：我仍然不知道。

S：我现在知道了，x=3，y=6.

给定N和M（两人总共说“我不知道”的次数，0<=M<=100），求所有可能的x和y。

 

乍一看，很神奇是吧？

我们来分析一下，3和6这个数据太大了（我刚开始分析了半天才看出了策略），先来看一些更简单的例子：

如果N=10，S知道x+y=3.可能的组合只有3=1+2. 毫无疑问x=1，y=2. 这样只要M=0次就能猜出来了。

如果N=10，x+y=8，x*y=7：

首先S知道所有可能的组合有8=2+6=3+5. S不知道到底是哪一种。

然后P知道所有可能的组合只有7=1*7. 因此x=1，y=7. 这样只要M=1次就猜出来了。

再如N=10，x+y=4，x*y=4.

首先S所有可能的组合为4=1+3=2+2. S还不知道。

然后P所有可能的组合为4=1*4=2*2. P也还不知道。

接下来又是S。他已经知道P在上一次没有猜出来，而如果x=1，y=3，那么P所知道的数应该为3，从而他应该在上一次就猜出来了。

所以4=1+3这种组合是错的。因此S知道，只有4=2+2一种可能。这样只要M=2次就可以猜出来了。

 

从以上例子我们可以得到一个结论：如果在第M次“不知道”后，某个人所知道的全部组合ALL中只有一种还没有猜出来（也就是说其余的所有组合形式都可以在M次之内猜出来），那么这一种就是答案无疑了。相反，如果在M次之后，仍然有多于一种情况都没有猜出来，显然也就无法决定到底是这些情况中的哪一种了。

 

根据这种策略，只要从m=0从小到大枚举，对每种x，y，根据M的奇偶性把x+y或x*y进行分解，看是否只有x，y这一种情况没有还猜出来。然后再M次之后，看所有xy的组合，那种可以在M次之后猜出来即可。

这道题的时限卡得很紧，有两点需要注意：

第一，如果某次在m次之后没有任何一个xy被更新的话，接下来m再怎么变大也不会有变动了（事实上有的状态无论多少次都猜不出来）。直接break

第二，x*y如果每次都分解一遍的话会被模和除运算卡掉。在开始预处理一下才行。

 

注：我这道题第一天看想出了正确解法，但是没注意到x和y不相等这个条件，得到的答案不对，所以搁置了两个星期。。今天拾起来再提。

刚开始WA掉了因为以为只有x=1，y=2这种情况是M=0的（其实N=10，x+y=19这种也是啊）。

 

 1 #include<iostream>
 2 #include<cstdio>
 3 #include<string>
 4 #include<cstring>
 5 #include<algorithm>
 6 #include<cmath>
 7 using namespace std;
 8 
 9 const int N = 205;
10 int n,f[N][N],tot[N*N];
11 pair<int,int> fac[N*N][15]/*15试出来的，否则MLE*/;
12 
13 int search(int a,int b,int dep)
14 {//搜索j，k是否要i次才能猜出来
15     int s = a+b,p = a*b,i,flag = 1;
16     if(dep%2==0) {
17         for(i = 1; i+i<s; i++)
18             if(i!=a && s-i<=n && f[i][s-i]>=dep) {
19                 flag = 0; break;
20             }
21         if(flag) f[a][b] = dep;
22     }
23     else {
24         for(i = 1; i<=tot[p]; i++)
25             if(fac[p][i].first!=a && f[fac[p][i].first][fac[p][i].second]>=dep) {
26                 flag = 0; break;
27             }
28         if(flag) f[a][b] = dep;
29     }
30     return f[a][b]==dep;
31 }
32 
33 void preset() 
34 {
35     for(int i = 1; i<n; i++)
36         for(int j = i+1; j<=n; j++) {
37             tot[i*j]++;
38             fac[i*j][tot[i*j]] = make_pair(i,j);
39         }
40 }
41 
42 int main()
43 {
44     int i,j,k,m,ans = 0;
45     cin >> n >> m;
46     memset(f, 0x3f, sizeof(f));
47     preset();
48     for(i = 0; i<=m; i++) {
49         int flag = 0;
50         for(j = 1; j<=n; j++)
51             for(k = j+1; k<=n; k++)
52                 if(f[j][k]>m)
53                     flag |= search(j,k,i);
54         if(!flag) break;
55     }
56     for(i = 1; i<=n; i++)
57         for(j = i+1; j<=n; j++)
58             if(f[i][j]==m) ans++;
59     printf("%d\\n", ans);
60     for(i = 1; i<=n; i++)
61         for(j = i+1; j<=n; j++)
62             if(f[i][j]==m) printf("%d %d\\n", i,j);
63     return 0;
64 }

View Code