线性代数(矩阵、高斯、线性基……)

矩阵

矩阵加法:

相同位置相加。

矩阵乘法:

满足分配率、结合律,不满足交换律(矩阵与逆矩阵之间除外)

矩阵转置:

记矩阵为 \(A\) ,则 \(A\) 的转置记为 \(A^T\)

性质:

  • \[{(A^T)}^T=A \]
  • \[{(A+B)}^T=A^T+B^T \]
  • \[{(k\times A)}^T=k\times A^T \]
  • \[{(AB)}^T=A^TB^T \]

矩阵求逆:

P4783 【模板】矩阵求逆

\[\begin{bmatrix}2&-1&0&1&0&0\\-1&2&-1&0&1&0\\0&-1&2&0&0&1\end{bmatrix} \]

对左半边的矩阵做高斯消元,同时更新右半边的部分,(交换时也一起交换,但最终不用再换回来了)。而做完之后的右半边部分就是求得的逆矩阵。

矩阵快速幂:

  • 对于不含常数项的递推式:(比较正常的矩阵快速幂)

  • 对于含有常数项的递推式:加上一维,在转移矩阵中不更改。

  • 对于含有关于 \(n^1\) 的递推式:加上两维,每次后一位给前一位加一。

  • 对于含有关于 \(n^k\) 的递推式:加上 \(k+1\) 维,例:

\[f[n]=2\times f[n-1]+n^k \]
\[\Downarrow \]
\[\begin{bmatrix}s[n]&n^k&n^{k-1}&\cdots&n^0\end{bmatrix}\times\begin{bmatrix}2&0&0&0&\cdots&0\\C_k^0&C_k^0&0&0&\cdots&0\\C_k^1&C_k^1&C_{k-1}^0&0&\cdots&0\\C_k^2&C_k^2&C_{k-1}^1&C_{k-2}^0&\cdots&0\\\vdots&\vdots&\vdots&\vdots&\ddots&\vdots\\C_k^k&C_k^k&C_{k-1}^{k-1}&C_{k-2}^{k-2}&\cdots&C_0^0\end{bmatrix}=\begin{bmatrix}s[n+1]&{(n+1)}^k&{(n+1)}^{k-1}&\cdots&{(n+1)}^0\end{bmatrix} \]

注意:!!

由于转移矩阵与答案举证的大小不同,应该在 struct 的矩阵中记录这个矩阵的大小,防止将 \(O(n^2)\) 变为 \(O(n^3)\) !!!

高斯消元

复杂度(朴素): \(O(n^3)\)

主要代码:

scanf("%d",&n);
for(int i=1;i<=n;i++) for(int j=1;j<=n+1;j++) scanf("%lf",&a[i][j]);
for(int i=1,Max=1;i<=n;Max=++i)
{for(int s=i+1;s<=n;s++) if(fabs(a[s][i])>fabs(a[Max][i])) Max=s; // 找出绝对值最大的 for(int j=1;j<=n+1;j++) swap(a[i][j],a[Max][j]);if(a[i][i]<10e-8 && a[i][i]>-10e-8) { p=false; break; } // 记得 double 的精度问题 for(int s=1;s<=n;s++) if(s!=i) // 这样省去了第二步处理的麻烦 {double tmp=0-(a[s][i]/a[i][i]);a[s][i]=0;for(int j=i+1;j<=n+1;j++) a[s][j]+=tmp*a[i][j];}
}
if(p) for(int i=1;i<=n;i++) printf("%.2lf\n",a[i][n+1]/a[i][i]);
else printf("No Solution\n");

线性基

线性基为一个数集构造出来的新数集,满足以下性质:

  • 线性基的元素能相互异或得到原集合的元素的所有相互异或得到的值

  • 线性基是满足性质 \(1\)最小的集合

  • 线性基没有异或和为 \(0\) 的子集。

  • 线性基中不同的异或组合异或出的数都是不一样的。

  • 线性基中每个元素的二进制最高位互不相同

用处:

  • 快速查询一个数是否可以被一堆数异或出来

  • 快速查询一堆数可以异或出来的最大 \(/\) 最小值

  • 快速查询一堆数可以异或出来的第 \(k\) 大值

处理线性基:

void Insert(ll x)
{for(int i=62;i>=0;i--){if(!(x & (1ll<<(ll)i))) continue; // 防止对高位影响 if(!p[i]) { p[i]=x; break; }x^=p[i]; // 更新 [0,i-1] 位的更优答案 }if(!x) zero=1ll; // 特判 0 
}

查询一个元素是否可以被异或出来:

bool ask(ll x)
{for(int i=62;i>=0;i--) if(x&(1ll<<(ll)i)) x^=p[i];return x==0;
}

查询异或最大值:

ll query_max()
{ll ret=0;for(int i=62;i>=0;i--) if((ans^p[i])>ans) ans^=p[i];return ans;
}

查询异或最小值:

ll query_min()
{for(int i=0;i<=62;i++) if(p[i]) return p[i];return 0;
}

查询异或第 \(k\) 小:

void rebuild()
{// 重建 d 数组,求出哪些位可以被异或为 1// d[i] 只有第 i 个二进制位为 1 for(int i=62;i>=1;i--) // 从高到低防止后效性 for(int j=i-1;j>=0;j--)if(p[i] & (1ll<<(ll)j)) p[i]^=p[j];for(int i=0;i<=62;i++) if(p[i]) d[cnt++]=p[i];
}
ll kth(ll k)
{if(!k) return 0ll; // 特判 0 if(k>=(1ll<<(ll)cnt)) return -1ll; // k 大于可以表示出的数的个数 ll ret=0;for(int i=62;i>=0;i--) if(k & (1ll<<(ll)i)) ret^=d[i];return ret;
}

变式:

O(logn) 求区间异或最大值:

P3292 [SCOI2016]幸运数字

题意:给定一棵树,求 \(x\)\(y\) 路径的异或最大值。

\(p[x][]\) 表示点 \(x\) 到根之间所有点的线性基,同时维护 \(pos[x][]\) 表示这一线性基由哪一个点转移而来。

\(\text{Dfs}\) 加入一个新的点时,贪心将贡献相同或更高,且深度更大的点代替原来线性基中的值。

这样查询的时候就能在 \(O(\log_{2} n)\) 复杂度内求出在点 \(x\)\(Lca\) 路径中的点的最大异或和(更深的点已经加入线性基)。

完整代码

CF1100F Ivan and Burgers

题意:求出序列中 \(l\)\(r\) 的区间最大异或和。

和上一题差不多,将维护更深的点转化为维护更靠后的点即可。

完整代码

线性基还可以推广至非二进制的情况。

P3265 [JLOI2015]装备购买

这里需要维护一个 \(k\) 进制的线性基

主要代码:

int p[Maxn]; // 由于浮点数存储不方便,这里 p 记录的是值的下标 
struct Data
{ld val[Maxk];int Cost;
}a[Maxn];
void Insert(Data k) // 这里传入的也是下标 
{for(int i=1;i<=m;i++){if(fabs(a[k].val[i])<10e-8) continue;if(!p[i]) { p[i]=k,cnt+=1,sum+=a[k].Cost; break; }ld tmp=a[k].val[i]/a[p[i]].val[i];for(int j=i;j<=m;j++) a[k].val[j]-=a[p[i]].val[j]*tmp;}
}

较复杂的情况

P4151 [WC2011]最大XOR和路径

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

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

相关文章

Tavan

Tavan–简单的k进制题 呵呵真的简单 也就是考试没做出来罢了&#xff0c;不慌不慌&#xff0c;奶我一口&#xff0c;还能活 【题目摘要】 题目描述 小 Zeljko 一直在阁楼里读他奶奶的旧信&#xff0c;并且发现了一个长度为 N 的单词。 不幸的是&#xff0c;由于溢出的墨水&…

.NET Core 2.x中使用Named Options处理多个强类型配置实例

来源&#xff1a; Using multiple instances of strongly-typed settings with named options in .NET Core 2.x作者&#xff1a; Andrew Lock译者&#xff1a; Lamond Lu.NET Core从1.0版本开始&#xff0c;就已经开始使用Options模式绑定强类型配置对象。从那时起到现在&…

CF235D-Graph Game【LCA,数学期望】

正题 题目链接:https://www.luogu.com.cn/problem/CF235D 题目大意 给出一棵基环树&#xff0c;每次随机选择一个点让权值加上这个点的连通块大小然后删掉这个点。 求删光所有点时期望权值。 1≤n≤30001\leq n\leq 30001≤n≤3000 解题思路 先找到环&#xff0c;然后考虑暴…

YBTOJ:魔法数字(数位dp)

文章目录题目描述解析题目描述 解析 迷惑。。。 首先&#xff0c;比较容易想到用二进制状态压缩记录1-9是否在十进制中出现过 然后就是整除的问题 如果记录余数&#xff0c;它的模数又有9个 开九维余数直接爆炸。。。 怎么办嘞&#xff1f; 有一个结论&#xff1a; 若&#xf…

斯特林数

第一类斯特林数 咕咕咕 第二类斯特林数 定义&#xff1a;把 \(n\) 个不同的球放入 \(r\) 个相同的盒子的方案数(盒子不能为空&#xff0c;记为&#xff1a;\(S(n,r)\) 或 \(\begin{Bmatrix}n\\r\end{Bmatrix}\) 。 递推式&#xff1a; \[\begin{Bmatrix}n\\r\end{Bmatrix}\begi…

Prosjecni(构造)

Prosjecni 【题目摘要】 描述 Slavko很无聊&#xff0c;所以他把正整数填到N*N的方阵中。 如果他填出来的方阵满足以下条件&#xff0c;他会特别高兴&#xff1a; ●每行中的数字的平均值是一个位于同一行的整数。 ●每列中的数字的平均值是一个位于同一列的整数。 ●表中的所…

撒花!中文翻译仓库链接已加入 ML.NET 官方示例网站首页

从2018年12月02日决定开始做ML.NET 示例中文版https://github.com/feiyun0112/machinelearning-samples.zh-cn&#xff0c;然后以每天一篇的速度进行翻译&#xff0c;总共耗时15天&#xff0c;将现有的官方实例全部翻译成了中文&#xff0c;并提交了添加中文链接PR&#xff0c;…

P3803 【模板】多项式乘法(FFT)

P3803 【模板】多项式乘法&#xff08;FFT&#xff09; 题目描述 给定一个 n 次多项式 F(x)&#xff0c;和一个 m 次多项式 G(x)。 请求出 F(x)和 G(x)的卷积。 从低到高输出F(x)*G(x)的系数 另一种问法&#xff1a; 如果有两个无限序列a和b&#xff0c;那么它们卷积的结果是…

CF666E-Forensic Examination【广义SAM,线段树合并】

正题 题目链接:https://www.luogu.com.cn/problem/CF666E 解题思路 给出一个串SSS和nnn个串TiT_iTi​。mmm次询问Sa∼bS_{a\sim b}Sa∼b​在Tl∼rT_{l\sim r}Tl∼r​中出现的最多次数并且输出这个串的编号。 1≤∣s∣≤5105,∑Ti≤5104,1≤m≤51051\leq |s|\leq 5\times 10^5…

简简单单组合数学

简简单单组合数学 P3158 [CQOI2011]放棋子 \(\uparrow\) 假组合数学&#xff0c;真 \(\text{DP}\) 。 \(f[i][j][k]\) : 用了 \(i\) 行 \(j\) 列&#xff0c;涂了前 \(k\) 种颜色的方案数。 \(g[i][j][k]\) : 用了 \(i\) 行 \(j\) 列&#xff0c;涂了第 \(k\) 种颜色的方案数(用…

YBTOJ:单词频率(AC自动机)

解析 我对力量一无所知 通过这题&#xff0c;可以看出我对AC自动机还是完全没有理解 qwq 首先容易想到&#xff1a; 建一课trie树&#xff0c;然后建树时记录每个串s的终点&#xff0c;这个点后面每被经过一次&#xff0c;就相当于出现一次该单词s 但是&#xff0c;这种“出现”…

一元一次方程

一元一次方程–逆波兰栈 【题目摘要】 题目描述 SLON是一个调皮的学生&#xff0c;为了让他静下心来&#xff0c;老师给他出了一道数学题&#xff1a; 给定表达式A&#xff0c;A中含有变量x和,-,*,(,)这些符号&#xff0c;括号成对出现&#xff0c;一个算术运算符均对应两个操…

.net core上 K8S(七).netcore程序的服务发现

正文上一章我们分享了k8s的网络代理模式&#xff0c;今天我们来分享一下k8s中的服务发现。1.环境变量模式的服务发现k8s默认为我们提供了通过环境变量来实现服务发现的功能&#xff0c;前提是1.需要service在pod之前创建2.适用于同一命名空间1.1创建servicekubectl create -f n…

质数与合数

质数与合数 题意&#xff1a; FFF和GGG正在玩一个质数与合数的游戏 一开始有N个石头 FFF和GGG轮流对这堆石头进行操作&#xff0c;FFF每次选择1到K之间的一个数x&#xff0c;并拿走x个石头&#xff0c;拿走之后剩下的石头数量必须是质数 接着GGG进行同样的操作&#xff0c;但…

P6628-[省选联考 2020 B 卷] 丁香之路【欧拉回路,最小生成树】

正题 题目链接:https://www.luogu.com.cn/problem/P6628 题目大意 给出nnn个点的一张完全无向图&#xff0c;i∼ji\sim ji∼j的边权是∣i−j∣|i-j|∣i−j∣。 然后给出mmm条必经边&#xff0c;和起点sss。 求对于每个终点经过所有必经边的最短路径。 1≤n≤2500,0≤m≤n(n…

ASP.NET Core 实战:使用ASP.NET Core Web API 和 Vue.js 搭建前后端分离项目

一、前言这几年前端的发展速度就像坐上了火箭&#xff0c;各种的框架一个接一个的出现&#xff0c;需要学习的东西越来越多&#xff0c;分工也越来越细&#xff0c;作为一个 .NET Web 程序猿&#xff0c;多了解了解行业的发展&#xff0c;让自己扩展出新的技能树&#xff0c;对…

【CF1199 D,E, F】Welfare State // Matching vs Independent Set // Rectangle Painting 1

2019-08-15下午三道练习题CF1199 思路有点难想 but很好实现 这是原网站链接&#xff1a;传送门 这里只完成D, E, F三题 文章目录D&#xff1a;Welfare State题目大意正解瞅瞅代码E:Matching vs Independent Set题目大意正解代码实现F:Rectangle Painting 1题目大意正解代码实现…

YBTOJ:前缀匹配(AC自动机)

文章目录题目描述解析代码题目描述 解析 做的不错 把trie树真的当成一棵树递归即可 注意一个标记时的问题&#xff1a; void AC(){int lstrlen(s01),pl1;for(int i1;i<l;i){int aask(s0[i]);pltr[pl][a];int kpl;while(k>1){if(ok[k]) break;//注意&#xff01;ok[k]1;…

[HAOI2006]均分数据

[HAOI2006]均分数据 题解&#xff1a; 题目稍微解释一下&#xff1a; 把n个数以分为m组&#xff0c;计算每一组的和&#xff0c;求得到的这m个数的方差。由于分法是任意的&#xff0c;我们要求这些方差中的最小值 我们先用STL中的函数random_shuffle()用来对一个元素序列进行…

开机指南

有些博客上了密码&#xff0c;密码是一个 \(100\) 以内的数字(与ACwisher有关) 腾讯QQ PC版 Notepad Div-c 编译选项&#xff1a; -Wl,--stack1024000000 -Wall -lm -O2 -stdc14 // Author:A weak man named EricQian #include<bits/stdc.h> using namespace std; #defin…