P4100-[HEOI2013]钙铁锌硒维生素【矩阵求逆,最大匹配】

正题

题目链接:https://www.luogu.com.cn/problem/P4100


题目大意

给出nnn个线性无关的向量AiA_iAi,然后给出nnn个向量BiB_iBi,求一个字典序最小的排列ppp使得将任意的AiA_iAi替换为BpiB_{p_i}Bpi后依旧线性无关。

1≤n≤3001\leq n\leq 3001n300


解题思路

首先因为我们有nnn个向量AAA线性无关,那么显然这nnn个向量能表示任意向量,如果对于一个BpiB_{p_i}Bpi替换为AiA_iAi后依旧线性无关,那么BpiB_{p_i}BpiAiA_iAi是等价的(因为BpiB_{p_i}BpiAiA_iAi都代表了剩下n−1n-1n1个无法表示的部分)。

所以只需考虑每个BjB_jBj能否换到AiA_iAi即可,构建出矩阵A=[A1,A2...An]A=[A_1,A_2...A_n]A=[A1,A2...An]B=[B1,B2...Bn]B=[B_1,B_2...B_n]B=[B1,B2...Bn],考虑一个置换矩阵使得AR=BAR=BAR=B,那么就是对于每个BBB如何用AAA进行表示。

那么如果Ri,j=0R_{i,j}=0Ri,j=0也就是说BBB可以用AjA_jAj以外的其他AAA表示所以BBB替换到AjA_jAj之后肯定线性有关了,所以不行。

R=BAR=\frac{B}{A}R=AB,求逆得到RRR,这样我们就知道哪些AAA可以替换哪些BBB了,问题就变成了最小字典序匹配。对于这个问题我们可以考虑找一条增广环就好了。

时间复杂度O(n3)O(n^3)O(n3)


code

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
using namespace std;
const int N=310;
const double eps=1e-8;
int n,v[N],link[N];
double A[N][N],B[N][N];
bool GetInv(){for(int i=1;i<=n;i++){int z=i;for(int j=i+1;j<=n;j++)if(fabs(A[j][i])>fabs(A[z][i]))z=j;swap(A[i],A[z]);swap(B[i],B[z]);double x=A[i][i];if(fabs(x)<eps)return 0;for(int j=1;j<=n;j++)A[i][j]/=x,B[i][j]/=x; for(int j=1;j<=n;j++){if(i==j)continue;double rate=-A[j][i];for(int k=1;k<=n;k++)A[j][k]+=rate*A[i][k],B[j][k]+=rate*B[i][k];}}for(int i=n;i>=1;i--)for(int j=1;j<i;j++){double rate=-A[j][i];for(int k=1;k<=n;k++)A[j][k]+=rate*A[i][k],B[j][k]+=rate*B[i][k];	}return 1;
}
bool dfs(int x){for(int i=1;i<=n;i++)if(!v[i]&&fabs(B[x][i])>=eps){v[i]=1;if(!link[i]||dfs(link[i])){link[i]=x;return 1;}}return 0;
}
int calc(int x,int top){for(int i=1;i<=n;i++)if(!v[i]&&fabs(B[x][i])>=eps){v[i]=1;if(link[i]==top||(link[i]>top&&calc(link[i],top))){link[i]=x;return i;}}return 0;
}
int main()
{scanf("%d",&n);for(int i=1;i<=n;i++)for(int j=1;j<=n;j++)scanf("%lf",&A[j][i]);for(int i=1;i<=n;i++)for(int j=1;j<=n;j++)scanf("%lf",&B[j][i]);if(!GetInv())return puts("NIE")&0;for(int i=1;i<=n;i++){memset(v,0,sizeof(v));if(!dfs(i))return puts("NIE")&0;}puts("TAK");for(int i=1;i<=n;i++){memset(v,0,sizeof(v));printf("%d\n",calc(i,i));}return 0;
}

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/318317.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

数数字

今天中午1小时&#xff0c;定时两道&#xff0c;我就知道自己的太阳降落了。。。orz 文章目录题目题解代码实现题目 PB 带来了若干只蒟蒻。 众所周知&#xff0c;NTF 是数论学会的会长&#xff0c;于是 PB 准备用数字击败 NTF&#xff0c;以证明 PB 比 NTF 更强。 于是 PB 准…

YBTOJ:前缀数组(KMP)

文章目录题目描述解析代码题目描述 解析 题面脸上写着5个大字&#xff1a;我是KMP 但是本题没有自己做出来。。。我一开始的思路其实很接近题解了&#xff0c;只是被我舍弃了qwq。 后来卡在暴力nL2的瓶颈上&#xff0c;用了个倍增的诡异操作搞到了nLlogL&#xff0c;但是n5&am…

Strange Partition CodeForces - 1471A

Strange Partition CodeForces - 1471A 题意&#xff1a; 对于数组b可以两两合并其中的相连元素&#xff0c;beauty值就是对于数组b中的b[i]除以k向上取整的和&#xff0c;求出beauty值的范围。 题解&#xff1a; 题目给的是向上取整&#xff0c;也就是越合并值有可能越低&…

西安活动 | 2019年1月13号 拥抱开源, 又见.NET 线下交流活动报名进行中

随着.NET Core的发布和开源&#xff0c;.NET又重新回到人们的视野。.NET Core的下个3.0即将release&#xff0c;加入非常多的新功能&#xff0c;越来越拥抱变化&#xff0c;DevOps和Microservice的最佳实践已经在.NET Core落地&#xff0c;比如 Ocelot网关、GrpcConsul 服务注册…

C++STL(set……)

set 底层实现是用红黑树。 set 建立 set<int> s; // 不可重&#xff0c;默认升序 set<int,less> s; // 不可重&#xff0c;升序 set<int,greater> s; // 不可重&#xff0c;降序 multiset<int> s; // 可重集 set 也可以重载&#xff0c;利用结构体实现…

[ZOJ 3203] 灯泡

今天这个**又带着题走来了&#xff01;&#xff01; 文章目录[TOC] 题目题解代码实现题目 相比 wildleopard 的家&#xff0c;他的弟弟 mildleopard 比较穷。他的房子是狭窄的而且在他的房间里面仅有一个灯泡。每天晚上&#xff0c;他徘徊在自己狭小的房子里&#xff0c;思考如…

P5591-小猪佩奇学数学【单位根反演】

正题 题目链接:https://www.luogu.com.cn/problem/P5591 题目大意 给出n,p,kn,p,kn,p,k求 (∑i0n(ni)pi⌊ik⌋)mod998244353\left(\sum_{i0}^n\binom{n}{i}p^i\left\lfloor\frac{i}{k}\right\rfloor \right)\mod 998244353(i0∑n​(in​)pi⌊ki​⌋)mod998244353 1≤n,p<9…

YBTOJ:字符串题(KMP)

文章目录题目描述解析1.p[i]>0:2.p[i]0&#xff1a;代码题目描述 解析 我个人做起来很费劲的一道题&#xff0c;用vector&#xff0c;并查集等等东西搞来搞去过掉了&#xff08;竟然只WA了一次 &#xff09; 看题解思路就一下子清晰了&#xff0c;还是对KMP的理解不到位。 …

Strange List CodeForces - 1471B

题意&#xff1a; 对于长度为n的数组a&#xff0c;从第一位开始如果可以整除x&#xff0c;就将x个a/x的结果加到数组最后&#xff0c;然后对下一位进行一样的操作&#xff0c;直到第x位不可以整除x&#xff0c;到此结束&#xff0c;然后计算此时数组的总和 题解&#xff1a; …

Asp.Net Core 轻松学-正确使用分布式缓存

前言本来昨天应该更新的&#xff0c;但是由于各种原因&#xff0c;抱歉&#xff0c;让追这个系列的朋友久等了。上一篇文章 在.Net Core 使用缓存和配置依赖策略 讲的是如何使用本地缓存&#xff0c;那么本篇文章就来了解一下如何使用分布式缓存&#xff0c;通过本章&#xff0…

计算几何初步

两点之间距离 欧氏距离 即欧几里得距离。 平面内两点的距离为 \[\sqrt{(x_1-x_2)^2(y_1-y_2)^2} \]立体空间内两点的距离为 \[\sqrt{(x_1-x_2)^2(y_1-y_2)^2(z_1-z_2)^2} \]\(\dots\) \(n\) 维空间内两点的距离为 \[\sqrt{\sum_{i1}^{n}{(x_1-x_2)^2}} \]曼哈顿距离 二维空间内…

树链剖分概念及模板 + 例题 [POJ3237 tree + 软件包管理器]

文章目录概念模板例题1&#xff1a;软件包管理器题目题解代码实现例题2&#xff1a;POJ3237 tree题目题解代码实现概念 树链剖分主要是用于解决以下这两个问题。 1、更改树上点x到点y的最短路径上的所有结点的值 2、查询树上点x到点y的最短路径上的所有结点的和。 在讲树链剖分…

P3706-[SDOI2017]硬币游戏【高斯消元,字符串hash】

正题 题目链接:https://www.luogu.com.cn/problem/P3706 题目大意 给出 nnn 个长度为 mmm 的 H/TH/TH/T 串。 开始一个空序列&#xff0c;每次随机在后面加一个 H/TH/TH/T &#xff0c;求每个串第一次出现的概率。 1≤n,m≤3001\leq n,m\leq 3001≤n,m≤300 解题思路 数据范…

YBTOJ:字符串匹配(KMP)

文章目录题目描述解析代码题目描述 解析 看了题解。。。 这题的关键在于可以变换匹配的一个充要条件&#xff1a; 每个字符与前一个相同字符的距离相同 这个搞出来之后就可以以它为关键字进行KMP了 注意&#xff01; 当与前一个字符的距离超过匹配长度时&#xff0c;是没有意义…

Strange Birthday Party CodeForces - 1471C

Strange Birthday Party CodeForces - 1471C 题意&#xff1a; 我有n个朋友&#xff0c;商店有m种商品&#xff0c;这m种商品按序号价格从小到大排列&#xff0c;对于每一个朋友我给出一个序号k&#xff0c;我可以直接给朋友序号k的商品价格的金钱或给朋友买一个序号小于k的商…

TransactionScope事务处理方法介绍及.NET Core中的注意事项

作者&#xff1a;依乐祝原文链接&#xff1a;https://www.cnblogs.com/yilezhu/p/10170712.html今天在写CzarCms的UnitOfWork的使用使用到了这个TransactionScope事务&#xff0c;因此对它进行了相关资料的查阅并记录如下&#xff0c;希望对大伙在.NET Core中使用有所帮助。写在…

P4548-[CTSC2006]歌唱王国【概率生成函数,KMP】

正题 题目链接:https://www.luogu.com.cn/problem/P4548 题目大意 ttt次询问&#xff0c;给出一个长度为mmm的串SSS和一个空串TTT&#xff0c;每次在TTT后面随机加入1∼n1\sim n1∼n的字符&#xff0c;直到TTT中出现SSS为止&#xff0c;求期望次数。 1≤n≤105,t≤50,1≤m≤10…

YBTOJ:字符匹配(KMP)

文章目录题目描述解析代码题目描述 解析 显然应该要尝试套kmp的板子 关键是如何套 也就是那个判断匹配的条件是什么的问题 本题的关键是当kmp匹配时&#xff0c;匹配位之前的所有位大小关系的顺序都是匹配的&#xff0c;所以我们只需要看当前位即可 考虑对b预处理出3个数组&a…

历史上的今天(history)+ 勇者斗恶龙(dragon)

朋友们我来了&#xff0c;好久都没有更新了&#xff0c;手实在痒痒&#xff0c;不撸两道&#xff0c;内心过意不去 文章目录A&#xff1a;历史上的今天&#xff08;history&#xff09;题目题解代码实现B&#xff1a; 勇者斗恶龙&#xff08;dragon&#xff09;题目描述题解代…

Codeforces Round #694 Div. 2

Codeforces Round #694 Div. 2 CodeForces 1471A Strange Partition CodeForces 1471B Strange List CodeForces 1471C Strange Birthday Party 签到题 稍有难度 CodeForces 1471D Strange Definition CodeForces 1471E Strange Shuffle CodeForces 1471F Strange Housing