BZOJ3659 : Which Dreamed It
2016-04-03 00:37
295 查看
首先判断一下是否无解,并剔除孤立点。
根据best theorem,有向图中以$i$为起点的欧拉回路个数为:
以$i$为根的树形图个数$\times\prod_{i=1}^n (deg(i)-1)!$。
根据matrix tree theorem,以$i$为根的树形图个数$=$基尔霍夫矩阵去掉第$i$行第$i$列的行列式。
$ans=以1为起点的欧拉回路个数\times 1的度数$。
高斯消元即可,时间复杂度$O(n^3\log P)$。
根据best theorem,有向图中以$i$为起点的欧拉回路个数为:
以$i$为根的树形图个数$\times\prod_{i=1}^n (deg(i)-1)!$。
根据matrix tree theorem,以$i$为根的树形图个数$=$基尔霍夫矩阵去掉第$i$行第$i$列的行列式。
$ans=以1为起点的欧拉回路个数\times 1的度数$。
高斯消元即可,时间复杂度$O(n^3\log P)$。
#include<cstdio> typedef long long ll; const int N=110,M=200010,P=1000003; int n,m,i,j,k,x,y,in ,ou ,vis ,g ;ll f[M],a ,ans; inline void read(int&a){char c;while(!(((c=getchar())>='0')&&(c<='9')));a=c-'0';while(((c=getchar())>='0')&&(c<='9'))(a*=10)+=c-'0';} void dfs(int x){ vis[x]=++m; for(int i=1;i<=n;i++)if(g[x][i]&&!vis[i])dfs(i); } inline void swap(ll&a,ll&b){ll c=a;a=b;b=c;} ll det(int n){ ll ans=1;bool flag=1; for(i=1;i<=n;i++)for(j=1;j<=n;j++)a[i][j]=(a[i][j]%P+P)%P; for(i=1;i<=n;i++){ for(j=i+1;j<=n;j++)while(a[j][i]){ ll t=a[i][i]/a[j][i]; for(k=i;k<=n;k++)a[i][k]=(a[i][k]+P-t*a[j][k]%P)%P; for(k=i;k<=n;k++)swap(a[i][k],a[j][k]); flag^=1; } ans=ans*a[i][i]%P; if(!ans)return 0; } if(!flag)ans=P-ans; return ans; } int solve(){ for(m=0,i=1;i<=n;i++)in[i]=ou[i]=vis[i]=0; for(i=0;i<=n;i++)for(j=0;j<=n;j++)a[i][j]=g[i][j]=0; int ed=0; for(i=1;i<=n;i++)for(read(k);k--;g[i][j]++)read(j),ed++; for(dfs(i=1);i<=n;i++)if(!vis[i]&&g[i])return 0; for(i=1;i<=n;i++)for(j=1;j<=n;j++)if(g[i][j]){ x=vis[i],y=vis[j]; ou[x]+=g[i][j];in[y]+=g[i][j]; a[x-1][y-1]-=g[i][j],a[x-1][x-1]+=g[i][j]; } for(i=1;i<=m;i++)if(in[i]!=ou[i])return 0; if(m==1)return f[g[1][1]]; ans=det(m-1)*in[1]; for(i=1;i<=m;i++)ans=ans*f[in[i]-1]%P; return ans; } int main(){ for(f[0]=i=1;i<M;i++)f[i]=f[i-1]*i%P; while(1){ read(n); if(!n)return 0; printf("%d\n",solve()); } }
相关文章推荐
- css
- <css 五>list
- Windows下的bat中文乱码问题
- 归并排序(MergeSort)
- html
- Android studio出现Error:Unable to tunnel through proxy. Proxy returns "HTTP/1.1 400 Bad Request"的解决办法
- 数据加密原理
- 0.Win8.1,Win10,Windows Server 2012 安装 Net Framework 3.5
- 关于sublime的一些使用方式
- CTF——你关注最新的漏洞吗
- 一步一步利用django创建博客应用(三)
- ssh整合context:component-scan包名写了*号:Failed to parse configuration class [org.springframework.cache.aspectj.AspectJJCacheConfiguration]
- swift系统自带二维码的生成
- 数据库原理(一)—— 关系代数(二)
- apache版本的hadoop源码编译
- C#-变量类型(值类型、引用类型)
- MindManger 2016 64位 破解版
- u-boot-2016.03在mini2440移植 之DM9000
- VS
- MySQL 数据库学习笔记1