某考试 T2 orzcyr
2024-09-08 15:00:45
非常nice的一道行列式的题目。
考虑如果没有路径不相交这个限制的话,那么这个题就是一个行列式:设 a[i][j] 为从编号第i小的源点到编号第j小的汇点的路径条数,那么矩阵a[][]的行列式就是的答案,因为行列式的定义就是给行一个列的排列,贡献就是所有a[i][p[i]]再乘上 (-1)^(p[] 这个排列的逆序对数).
但是路径不相交就很恶心。。。。根本没法分开算嘛。。。。
不过逆序对可是有一个特殊性质的: 如果把 p[i] 和 p[j] swap一下,那么这个排列的逆序对数的变化值一定是奇数。
这个不难证明,因为仅有权值和下标都在交换的两个数的中间的那些数会产生逆序对变化,但是变化都是双倍的,所以仅有 (p[i].p[j]) 造成了 +/- 1的影响是奇数。
然后我们可以发现两条相交的路径 (a,b) , (c,d) 我们把交点后面的路径 swap 一下,那么就是 (a,d) , (c,b)了,原理就是上述了。
#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cmath>
#include<cstring>
#include<cstdlib>
#include<queue>
#include<ctime>
#define ll long long
using namespace std;
const int maxn=605;
int a[maxn][maxn],hd[maxn],num,f[maxn];
int to[maxn*100],ne[maxn*100],X,Y,d[maxn];
int id[maxn],od[maxn],n,m,p,dy[maxn];
int ans=1; inline void addline(int x,int y){ id[y]++,od[x]++,to[++num]=y,ne[num]=hd[x],hd[x]=num;}
inline int add(int x,int y){ x+=y; return x>=p?x-p:x;}
inline void ADD(int &x,int y){ x+=y; if(x>=p) x-=p;}
inline int mul(int x,int y,const int ha){ return x*(ll)y%ha;}
inline int ksm(int x,int y){
int an=1;
for(;y;y>>=1,x=mul(x,x,p)) if(y&1) an=mul(an,x,p);
return an;
} inline void build(int u){
queue<int> q; int x;
for(int i=1;i<=n;i++) if(!id[i]) q.push(i); memcpy(d,id,sizeof(id)); while(!q.empty()){
x=q.front(),q.pop();
// cout<<x<<' '<<f[x]<<endl;
for(int i=hd[x];i;i=ne[i]){
ADD(f[to[i]],f[x]);
if(!(--d[to[i]])) q.push(to[i]);
}
} for(int i=1;i<=n;i++) if(!od[i]) a[u][dy[i]]=f[i];
} inline void xy(){
/*
for(int i=1;i<=n;i++){
for(int j=1;j<=n;j++) printf("%d ",a[i][j]);
puts("");
}
*/ for(int i=1,inv,tmp;i<=n;i++){
if(!a[i][i]){
ans=p-ans;
for(int j=i+1;j<=n;j++) if(a[j][i]){
for(int k=i;k<=n;k++) swap(a[i][k],a[j][k]);
break;
}
} ans=mul(ans,a[i][i],p);
inv=ksm(a[i][i],p-2); for(int j=i+1;j<=n;j++) if(a[j][i]){
tmp=a[j][i]*(ll)inv%(const int)p;
for(int k=i;k<=n;k++) ADD(a[j][k],p-mul(a[i][k],tmp,p));
}
}
} inline void solve(){
for(int i=1;i<=n;i++) if(!od[i]) dy[i]=++Y;
for(int i=1;i<=n;i++) if(!id[i]){
memset(f,0,sizeof(f)),f[i]=1;
X++,build(X);
} /*
for(int i=1;i<=X;i++){
for(int j=1;j<=X;j++) printf("%d ",a[i][j]);
puts("");
}
*/ n=X,xy();
} int main(){
freopen("orzcyr.in","r",stdin);
freopen("orzcyr.out","w",stdout); scanf("%d%d%d",&n,&m,&p); int uu,vv;
while(m--) scanf("%d%d",&uu,&vv),addline(uu,vv); solve(); printf("%d\n",ans);
// cout<<X<<' '<<Y<<endl;
return 0;
}
最新文章
- Android系统build.prop文件
- MongoDB学习笔记~MongoVUE对数据进行查询,排序和按需显示
- spring功能总结
- 12款界面精美的 HTML5 &; CSS3 网站模板
- Linux变量
- Linux学习 : 裸板调试 之 使用MMU
- docker入门(一)
- Android --- 字符串\n的换行问题
- Linux内核——定时器和时间管理
- .NET反编译之Reflector基础示例
- Android 开发笔记___图像按钮__imageButton
- TypeScript系列 - 什么是TypeScript
- C++(实验三)
- redis 的使用,及如何使用redis维护数亿人的登录状态
- 2017-12-15python全栈9期第二天第五节之while else的用法二当不被break打断时else内容的结果会被打印
- Android生成自定义二维码
- Retrofit 代理模式
- Collection集合总结,List和set集合的用法,HashSet和LinkedHashSetde用法
- SpringMvc+hibernate+easyui简单的权限管理系统
- HDU2157(SummerTrainingDay05-F dp)