51nod 1318 最大公约数与最小公倍数方程组 2-SAT+数学

题意

求解一个有趣的方程组，方程组有N个未知正整数，x[0],x[1],x[2],…,x[N-1]。
方程组由M个方程组成，方程只有两种类型：
1）GCD(x[i],x[j]) = G ；其中，i != j，且GCD(a,b)为正整数a与b的最大公约数的值
2）LCM(x[i],x[j]) = L ；其中，i != j，且LCM(c,d)为正整数c与d的最小公倍数的值
你需要判断这样一个方程组是否存在解，且N个未知数都是正整数。
问题包含多组测试数据。
输入数据的第一行是一个整数T，表示有T组测试数据，1<=T<=10.
之后有T组相同结构的数据。
每组数据的第一行有两个整数N、M，分别表示未知数个数与方程个数，其中1<=N<=200,1<=M<=200.
之后M行，每行表示一个方程。
第i个方程格式为：”ti ai bi ci”,其中，ti为一个字符，ai，bi为[0,N-1]上的整数，1<=ci<=1,000,000,000。
如果ti=’L’,表示该方程为 LCM(x[ai],x[bi])=ci；
如果ti=’G’,表示该方程为 GCD(x[ai],x[bi])=ci。

分析

这题是真的难打。。。
思路其实并不难。我们对每一种素数分开处理。设当前处理到素数p，p在x中的指数为p1，在y中的指数为p2，LCM(x,y)则表示max(p1,p2)=c，GCD(x,y)则表示min(p1,p2)=c。
可以把每个变量拆成k个点对，其中的第i个点对(xi,yi)，xi表示该变量选p^i，yi则表示不选。2-SAT建图之后跑一下即可。

代码

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

const int N=205;

int n,m,zhi,tot,tot_block,dfn[N*N],low[N*N],stack[N*N],top,bel[N*N],last[N*N],prime[N*N],tim,cnt,mx;
bool del[N*N],ins[N*N],val[N*N],vis[N*N];
struct edge{
  int to,next;}e[N*N*N];
struct data{
  int op,a,b,c;}pro[N];

int read()
{
    int x=0,f=1;char ch=getchar();
    while (ch<'0'||ch>'9'){
  if(ch=='-')f=-1;ch=getchar();}
    while (ch>='0'&&ch<='9'){
  x=x*10+ch-