cogs2060 除法表达式【解法二】

【题目描述】

经过无尽的蘑菇和无尽的不知道错在哪里、不知道发生了什么只下，恭喜KZ成功的造出了数据！！可喜可贺可喜可贺

给出一个这样的除法表达式：X1/X2/X3/X4/X5/X6…./Xk,

正常的除法表达式是从左往右计算的，但是我们可以向其中添加括号，以改变它的运算顺序，如1/2/1/2=1/4 ,
(1/2)/(1/2)=1 ;

那么， 给定一个除法表达式，能否通过添加括号使它的值为正整数？

【输入格式】

多行

每行一个 除法表达式 【输出格式】

多行

对应输入，分别输出“YES”（能为正整数）或“NO”（不能为正整数）

解法一见【这里】。
和解法一的结论一样，但是不直接求gcd约分，而是对x2进行质因数分解，用其他数消去对应因数。

#include<cstdio>
#include<cstring>
int prm[100010],cnt[100010];
char s[10000010];
int main()
{
    freopen("baoquansl.in","r",stdin);
    freopen("baoquansl.out","w",stdout);
    int i,j,k,l,m,n,p,q,x,y,z;
    bool ans;
    while (scanf("%s",s+1)==1)
    {
        l=strlen(s+1);
        y=0;
        for (i=1;s[i]>='0'&&s[i]<='9';i++)
          y=y*10+s[i]-'0';
        x=0;
        for (i++;s[i]>='0'&&s[i]<='9';i++)
          x=x*10+s[i]-'0';
        k=0;
        for (j=2;x>1;j++)
          if (x%j==0)
          {
            prm[++k]=j;
            cnt[k]=0;
            while (x%j==0)
            {
                cnt[k]++;
                x/=j;
            }
          }
        for (j=1;j<=k;j++)
          while (y%prm[j]==0&&cnt[j])
          {
              y/=prm[j];
              cnt[j]--;
          }
        y=0;
        for (i++;i<=l;i++)
          if (s[i]>='0'&&s[i]<='9')
            y=y*10+s[i]-'0';
          else
          {
            for (j=1;j<=k;j++)
              while (y%prm[j]==0&&cnt[j])
              {
                y/=prm[j];
                cnt[j]--;
              }
            y=0;
          }
        for (j=1;j<=k;j++)
          while (y%prm[j]==0&&cnt[j])
          {
              y/=prm[j];
              cnt[j]--;
          }
        ans=1;
        for (i=1;i<=k;i++)
          if (cnt[i])
          {
            ans=0;
            break;
          }
        if (ans) printf("YES\n");
        else printf("NO\n");
    }
}