文章目錄

• 解析
• 代碼

傳送門

解析

關(guān)鍵性質(zhì)是一個定理：若m>3*n-6，必然不存在合法的平面圖
這誰知道啊
不過這題應(yīng)該往也許圖過于稠密時必然無解這方面想
所以我們只需要考慮m、n同階的情況就行了
這個時候我們直接暴力判斷跑2-SAT就行了

代碼

#include<bits/stdc++.h> using namespace std; #define ll long long const int N=35605; const int M=30005; const int mod=998244353; const double eps=1e-6; ll read(){ll x=0,f=1;char c=getchar();while(!isdigit(c)){if(c=='-') f=-1;c=getchar();}while(isdigit(c)){x=x*10+c-'0';c=getchar();}return x*f; } int n,m; struct node{int to,nxt; }p[2005000]; int fi[N],cnt; void addline(int x,int y){p[++cnt]=(node){y,fi[x]};fi[x]=cnt;//printf("x=%d y=%d\n",x,y); } int u[M],v[M],pos[N]; bool judge(int i,int j){return (u[i]<u[j]&&v[i]>u[j]&&v[i]<v[j])||(v[i]>v[j]&&u[i]>u[j]&&u[i]<v[j]); } int dfn[N],low[N],zhan[N],col[N],tim,tot,top; void tarjan(int x){dfn[x]=low[x]=++tim;zhan[++top]=x;for(int i=fi[x];~i;i=p[i].nxt){int to=p[i].to;if(!dfn[to]){tarjan(to);low[x]=min(low[x],low[to]);}else if(!col[to]) low[x]=min(low[x],dfn[to]);}if(low[x]==dfn[x]){col[x]=++tot;while(zhan[top]!=x){col[zhan[top--]]=tot;}top--;}return; } void init(){memset(fi,-1,sizeof(fi));cnt=-1;memset(dfn,0,sizeof(dfn));memset(low,0,sizeof(low));memset(col,0,sizeof(col));memset(zhan,0,sizeof(zhan));tim=top=0;tot=2*m; } int main(){int T=read();while(T--){n=read();m=read();init();for(int i=1;i<=m;i++){u[i]=read(),v[i]=read();}for(int i=1;i<=n;i++){int x=read();pos[x]=i;}if(m>3*n-6){printf("NO\n");continue;}for(int i=1;i<=m;i++){u[i]=pos[u[i]];v[i]=pos[v[i]];if(u[i]>v[i]) swap(u[i],v[i]);//printf("%d: (%d %d)\n",i,u[i],v[i]);}for(int i=1;i<=m;i++){//if(v[i]==u[i]+1) continue;for(int j=i+1;j<=m;j++){//if(v[j]==u[j]+1) continue;if(judge(i,j)){addline(i,j+m);addline(j,i+m);addline(j+m,i);addline(i+m,j);//printf("i=%d j=%d ok\n",i,j);}}}for(int i=1;i<=m*2;i++){if(!dfn[i]) tarjan(i);}int flag=1;for(int i=1;i<=m;i++){if(col[i]==col[i+m]){flag=0;break;}}if(flag) printf("YES\n");else printf("NO\n");} } /* 2 1 4 */

總結(jié)

以上是生活随笔為你收集整理的YBTOJ洛谷P3209：平面图判定（2-SAT）的全部內(nèi)容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網(wǎng)站內(nèi)容還不錯，歡迎將生活随笔推薦給好友。