6.10 省选模拟赛 小C的利是 高斯消元 矩阵行列式
LINK:小C的利是
想起来把这道题的题解写了 。一个常识:利是在广东那边叫做红包。
关于行列式的题目 不过我不太会23333..口胡还是可以的。
容易想到10分的状压.不过没什么意思。
仔细观察要求的东西 在每一行中选择一个数字 选择的位置还是相应的排列不过这个是排列之和.
容易联想到行列式的那个定义式。
此时容易发现如果把每个位置上的数字变成 \(x^{a_{i,j}}\)
那么就把乘法变成了加法 也就是最后求出来的行列式是一个nk多项式.
直接利用拉格朗日插值法 那么就得到了一个\(n^4k\)的做法。
不过这个有很大的几率不对 因为存在一个\((-1)^k\)的东西 可能抵消掉了
所以为了不让其抵消掉 那么此时可以rand一个随机的值 然后再做 这样几乎是卡不掉的。
考虑进行优化 期望得到一个\(n^3k\)的做法或者更优做法。
满分做法:
考虑一个质数P 其中P=sk+1.可以证明必然存在且不是很大。
然后找到其原根 此时就可以得到K次单位根了。
将这k次单位根带入x 然后将求出的行列式相加 然后就能得到 \(x^0\)位置的系数的K倍。
证明:
设一个排列p 最后的和为S 那么最后的贡献为 \(1+g^S+g^{2S}+g^{3S}+...\) (系数省略 也显然系数是相同的。
考虑 S%k==0时 上面的每一项都为1 系数得到累加。
当S%k!=0时 等比数列求和 得到\(\frac{1-g^{(p-1)S}}{1-g}\) 得到0.
那么其实最后我们只需要判断是不是0即可 可以忽略k倍这个东西。
//#include<bits\stdc++.h>
#include<iostream>
#include<iomanip>
#include<cstdio>
#include<cstring>
#include<string>
#include<ctime>
#include<cmath>
#include<cctype>
#include<cstdlib>
#include<queue>
#include<deque>
#include<stack>
#include<vector>
#include<algorithm>
#include<utility>
#include<bitset>
#include<set>
#include<map>
#define ll long long
#define db double
#define INF 1000000000000000000ll
#define ldb long double
#define pb push_back
#define put_(x) printf("%d ",x);
#define get(x) x=read()
#define gt(x) scanf("%d",&x)
#define gi(x) scanf("%lf",&x)
#define put(x) printf("%d\n",x)
#define putl(x) printf("%lld\n",x)
#define gc(a) scanf("%s",a+1)
#define rep(p,n,i) for(RE int i=p;i<=n;++i)
#define go(x) for(int i=lin[x],tn=ver[i];i;tn=ver[i=nex[i]])
#define fep(n,p,i) for(RE int i=n;i>=p;--i)
#define vep(p,n,i) for(RE int i=p;i<n;++i)
#define pii pair<int,int>
#define mk make_pair
#define RE register
#define P 1000000007
#define gf(x) scanf("%lf",&x)
#define pf(x) ((x)*(x))
#define uint unsigned long long
#define ui unsigned
#define EPS 1e-8
#define sq sqrt
#define S second
#define F first
#define Set(a,v) memset(a,v,sizeof(a))
using namespace std;
char buf[1<<15],*fs,*ft;
inline char getc()
{
return (fs==ft&&(ft=(fs=buf)+fread(buf,1,1<<15,stdin),fs==ft))?0:*fs++;
}
inline int read()
{
RE int x=0,f=1;RE char ch=getc();
while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getc();}
while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getc();}
return x*f;
}
const int MAXN=110,maxn=55;
int n,k,mod;
int a[MAXN][MAXN];
int b[MAXN][MAXN];
int c[MAXN][MAXN];
vector<int>d;
inline int ksm(int b,int p)
{
if(p<0)return 0;
int cnt=1;
while(p)
{
if(p&1)cnt=(ll)cnt*b%mod;
b=(ll)b*b%mod;p=p>>1;
}
return cnt;
}
inline int check_p(int x)
{
for(int i=2;i*i<=x;++i)if(x%i==0)return 0;
return 1;
}
inline int check_g(int x)
{
vep(0,d.size(),i)if(ksm(x,(mod-1)/d[i])==1)return 0;
return 1;
}
inline int det()
{
int ans=1;
rep(1,n,i)
{
int p=i;
rep(i+1,n,j)if(abs(c[j][i])>abs(c[p][i]))p=j;
if(p!=i)
{
rep(1,n,j)swap(c[i][j],c[p][j]);
ans=(mod-ans)%mod;
}
ans=((ll)ans*c[i][i]+mod)%mod;
int ww=ksm(c[i][i],mod-2);
rep(i+1,n,j)
{
int cc=(ll)ww*c[j][i]%mod;
rep(1,n,k)
{
c[j][k]=(c[j][k]-(ll)cc*c[i][k]+mod)%mod;
}
}
}
return ans;
}
signed main()
{
freopen("luckymoney.in","r",stdin);
freopen("luckymoney.out","w",stdout);
srand(time(0));
get(n);get(k);
rep(1,n,i)rep(1,n,j)get(a[i][j]),b[i][j]=rand();
mod=k+1;while(!check_p(mod))mod+=k;
int ww=mod-1;
for(int i=2;i*i<=ww;++i)
if(ww%i==0)
{
d.pb(i);
while(ww%i==0)ww/=i;
}
if(ww>1)d.pb(ww);
int g=2;
while(!check_g(g))++g;
int wn=ksm(g,(mod-1)/k);//mod的k次单位根.
int D=1,ans=0;
vep(0,k,T)
{
rep(1,n,i)
rep(1,n,j)
{
c[i][j]=(ll)b[i][j]*ksm(D,a[i][j])%mod;
}
ans=(ans+det())%mod;
D=(ll)D*wn%mod;
}
if(ans)puts("Yes");
else puts("No");
return 0;
}
6.10 省选模拟赛 小C的利是 高斯消元 矩阵行列式的更多相关文章
- 2019-ACM-ICPC-沈阳区网络赛-K. Guanguan's Happy water-高斯消元+矩阵快速幂
2019-ACM-ICPC-沈阳区网络赛-K. Guanguan's Happy water-高斯消元+矩阵快速幂 [Problem Description] 已知前\(2k\)个\(f(i)\),且 ...
- 【BZOJ-4031】小z的房间 Matrix-Tree定理 + 高斯消元解行列式
4031: [HEOI2015]小Z的房间 Time Limit: 10 Sec Memory Limit: 256 MBSubmit: 937 Solved: 456[Submit][Statu ...
- 5.19 省选模拟赛 小B的图 最小生成树 LCT
LINK:小B的图 这道题就比较容易了. 容易想到将询问离线 然后 从小到大排序 那么显然是优先放正图(x+k)的边. 考虑随着x的增大 那么负图上的边会逐渐加进来 一条边被加进来当且仅当 其权值小于 ...
- 5.19 省选模拟赛 小B的夏令营 概率 dp 前缀和优化dp
LINK:小B的夏令营 这道题是以前从没见过的优化dp的方法 不过也在情理之中. 注意读题 千万不要像我这个sb一样 考完连题意都不知道是啥. 一个长方形 要求从上到下联通的概率. 容易发现 K天只是 ...
- 5.10 省选模拟赛 拍卖 博弈 dp
LINK:拍卖 比赛的时候 前面时间浪费的有点多 写这道题的时候 没剩多少时间了. 随便设了一个状态 就开始做了. 果然需要认真的思考.其实 从我的状态的状态转移中可以看出所有的结论. 这里 就不再赘 ...
- 5.10 省选模拟赛 tree 树形dp 逆元
LINK:tree 整场比赛看起来最不可做 确是最简单的题目. 感觉很难写 不过单独考虑某个点 容易想到树形dp的状态. 设f[x]表示以x为根的子树内有黑边的方案数. 白边方案只有一种所以不用记录. ...
- 【洛谷比赛】[LnOI2019]长脖子鹿省选模拟赛 T1 题解
今天是[LnOI2019]长脖子鹿省选模拟赛的时间,小编表示考的不怎么样,改了半天也只会改第一题,那也先呈上题解吧. T1:P5248 [LnOI2019SP]快速多项式变换(FPT) 一看这题就很手 ...
- 10.16 NOIP模拟赛
目录 2018.10.16 NOIP模拟赛 A 购物shop B 期望exp(DP 期望 按位计算) C 魔法迷宫maze(状压 暴力) 考试代码 C 2018.10.16 NOIP模拟赛 时间:2h ...
- 10.17 NOIP模拟赛
目录 2018.10.17 NOIP模拟赛 A 咒语curse B 神光light(二分 DP) C 迷宫maze(次短路) 考试代码 B 2018.10.17 NOIP模拟赛 时间:1h15min( ...
随机推荐
- 如何嵌套一个网页html到另一个html中
在常规网页开发中(单页应用除外哈),经常会遇到把一些通用内容的页面集中到一个页面中,需要使用这些页面只需要包含引入即可,这样有利于维护和修改,当通用页面修改时只需更改一个文件就可以了,不需要每个文件单 ...
- box-shadow,text-shadow
box-shadow:inset 30px 40px 20px #f00; 如上实例,总共五个参数,其中第一个代表阴影是向内阴影还是向外阴影,第二个参数代表向右(从左向右)的偏移量,第三个参数代表向下 ...
- 状压DP之吃奶酪
题目 传送们 思路 1≤n≤15,妥妥的状压,数据这么小, 这道题的状压思路还是很好想的,我们定义f[i][s]代表以i为起点,吃掉状态为s的奶酪所需要跑的最短距离,那么显然,我们先枚举状态s,然后枚 ...
- Flask-实现下载功能
1. 接口返回真实的文件 这种情况比较简单, flask里带有此类api, 可以用send_from_directory和send_file. 核心代码如下: from flask import se ...
- URL编码转换函数:escape()、encodeURI()、encodeURIComponent()讲解
转自:https://www.cnblogs.com/douJiangYouTiao888/p/6473874.html 函数出现时间: escape() ...
- LintCode笔记 - 145.大小写转换 - 极简之道 - 最短代码
这道题目一眼就能看出是送分题,当然在这里也不谈高难度的实现逻辑,肯定有同学会想直接用自带函数实现不就可以了吗? 对的,就是这么简单,然而今天的重点是如何把代码简写到最短. 本文章将带你把代码长度从 一 ...
- IIFE中的函数是函数表达式,而不是函数声明
下面的代码打印什么内容,为什么? var b = 10; (function b(){ b = 20; console.log(b); })(); 针对这题,在知乎上看到别人的回答说: 函数表达式与函 ...
- JVM 学习笔记(四)
回顾: 在之前的文章中,我们主要体现了当堆内存设置的比较小的情况下,比如:-Xmx20M -Xms20M,在项目运行的过程中,不断往内存中去添加对象, 这时候就会出现OOM,也就是内存溢出,本文章将展 ...
- python eval函数,将列表样式的字符串转化为列表
python eval函数,将列表样式的字符串转化为列表 >>> str_1 = '[1,2,3,4,5,6]'>>> type(str_1)<type 's ...
- 当我谈 HTTP 时,我谈些什么?
当我们打开网站时也许不会去留意网站前面的HTTP是怎么来的.但是它毫无疑问在网络中有着举足轻重的地位.本文从起源到发展,详说HTTP从1到3的演变. 说在前面 本文不致力于讲完 HTTP 的全部内容, ...