ZOJ 1003 Crashing Ballon 解题思路

前言

此题的思路参考了 成鹏致远 的博客 http://infohacker.blog.51cto.com/6751239/1195118，仅对原博的解题思路进行一些补充。

此篇思路亦是我的真实思路不免有些地方绕了圈子毕竟是菜鸡 希望以此为警惕，在之后的一些情况里尽量避免。

原题题目及翻译省略。

抽象过程

先构造数学模型。100个气球可以看成一个集合

M={1，2，3，…，98，99，100}，

两位参赛者分别用X₁和X₂表示，则

X₁=M₁*M₂*M₃ …*M_n($n\ge 100$),

X₂=M^’₁*M^’₂*M^’₃…*M^’_m($m\le 100-n$)，

并且

X₁$\ge$X₂。

列出这几个式子并不能帮助我们解决问题，所以我们只能再回到题目中去。
人家已经提醒我们有几种情况了，分别是
一、X₁为真同时X₂为假
二、X₁为假同时X₂为真
三、X₁为假且X₂亦假
而此题的解题目标就是在X₁和X₂为任意符合题意的取值下判断这一组数属于哪种情况。那么我们结合上面所列的式子带入这几种情况。

X₁为真X₂为假

假设X₁为真同时X₂为假，则X₁和X₂的所有因子M₁、M₂、M₃ …、M_n、M^’₁、M^’₂、M^’₃…、M^’_m各不相同。那么此时将此组数的乘积除以1~100中的某些个数的结果应该可以1。即
$$(\prod _{n=1}^n{M}_n\cdot \prod _{m=1}^n{M}_m)÷(x_1{x}_2{x}_3...)=1$$
其中x₁、x₂…皆属于[1,100]。

到这里，应该就可以利用循环遍历对此种情况进行判断了，因为X₁和X₂可以看成是各个因子相乘得到的，直接从2开始对X₁、X₂进行遍历做除法就可以了。
这样真的可行吗？当一个中间思路完成时，一定要三思而后行，对算法进行整理优化，否则会加大后续工作的难度。
若按上述流程进行判断，我们的确可以将第一种情况判断出来，但是若最后结果不为1，那么情况二和情况三会发生混淆，无法进行进一步判断。
那么该如何改进呢？其实稍微用心观察情况会发现X₁、X₂是相互独立的，那么不妨将上面的那一串式子改成$$\prod _{n=1}^n{M}_n÷(\prod P)=1$$ &&$$\prod _{m=1}^n{M}_m÷(\prod Q)=1$$ (不要在意格式 )
其中P∪Q∈{1、2、3…、100}&&P∩Q=$\varphi$。
那么在改动之后，代码实现时就需要在遍历所有因子时尝试分别除以X₁和X₂直至X₁和X₂都为1。这样就可以证明声称自己获得较高分数的X₁的确说的是真话。
第一种情况的思考没有结束，有一个疑问出现了，那就是为什么情况一成立时结果一定还要X₂为1呢？首先题目中已经规定得分较小者必说实话，并且可以稍稍思考一下——若X₂不为1，则X₂在提出挑战时说了谎，那为啥你说了谎还要比X₁小呢？（不从实际情况考虑）从题意和逻辑都说明此种情况出现了矛盾，因此结果都必须为1。

X₁为假X₂为真

若 X₁为假X₂为真，仿照已有思路，则
$$\prod _{n=1}^n{M}_n÷(\prod P)！=1$$ &&$$\prod _{m=1}^n{M}_m÷(\prod Q)=1$$
其中P∪Q∈{1、2、3…、100}&&P∩Q=$\varphi$。

X₁为假X₂为假

$$\prod _{n=1}^n{M}_n÷(\prod P)！=1$$ &&$$\prod _{m=1}^n{M}_m÷(\prod Q)！=1$$
其中P∪Q∈{1、2、3…、100}&&P∩Q=$\varphi$。

思路完善

三种情况的思路已然明晰，但是在实现过程中我们可以想象到——三种情况是在同时进行判断的。那就会产生一个新问题——优先级谁最高？只要X₁和X₂的因子存在各不相同的可能，就可判定X₁说了实话，反之，若情况一不成立，则无法立刻判断是否为情况二或情况三，需要遍历所有情况然后才能进行判断。
关于选定X₁、X₂因子的方法，因为是“一对多”的树状结构，首先联想到的实现方法应该是DFS和BFS，又并且有情况一的存在，所以用DFS较快。

代码

#include <iostream>
#include<algorithm>
using namespace std;
int a_iftrue=0,b_iftrue=0;
void dfs(int x,int y,int i)
{
    if(x==1&&y==1)
    {
        a_iftrue=1;
        return;
    }
    if(y==1)
        b_iftrue=1;
    while(i<=100&&(i<=x||i<=y))
    {
        if(x%i==0)
        {
            dfs(x/i,y,i+1);
            if(a_iftrue==1)
                return;
        }
        if(y%i==0)
        {
            dfs(x,y/i,i+1);
            if(a_iftrue==1)
                return;
        }
        i++;
    }
}
int main()
{
    int a=0,b=0;
    while(cin>>a>>b)
    {
        if(a<b)
        {
            int c=0;
            c=a;
            a=b;
            b=c;
        }
        a_iftrue=0,b_iftrue=0;
        dfs(a,b,2);
        if(a_iftrue==0&&b_iftrue==1)
            cout<<b<<endl;
        else
            cout<<a<<endl;
    }
    return 0;
}

代码注释就直接看成鹏致远大佬的博客就行了，我写的凑活AC，没啥参考价值。