BZOJ4128: Matrix 矩阵BSGS

最新推荐文章于 2018-04-11 20:18:38 发布

Drenight

最新推荐文章于 2018-04-11 20:18:38 发布

阅读量203

点赞数

分类专栏： BSGS

本文链接：https://blog.csdn.net/qq_34927456/article/details/77918827

版权

博客介绍了如何利用矩阵快速幂(BSGS算法)解决指数模运算问题，避免求逆元的过程，从而提高计算效率。

摘要由CSDN通过智能技术生成

A^x=B(mod C)

令x=im-j

A^(im)=BA^j(mod C)

这就不用求逆元了

#include<bits/stdc++.h>
#include<stdio.h>
#include<algorithm>
#include<queue>
#include<string.h>
#include<iostream>
#include<math.h>
#include<set>
#include<map>
#include<vector>
#include<iomanip>
using namespace std;
#define ll long long
#define ull unsigned long long
#define pb push_back
#define FOR(a) for(int i=1;i<=a;i++)
const int inf=0x3f3f3f3f;
const ll Linf=9e18;
const int maxn=1e6+7; 
const ll mod=1e9+7;
const double eps=1e-6;

#define base 231
ll n,p;
struct MATRIX{
	ll x[80][80];
	void init(){memset(x,0,sizeof x);}
	MATRIX operator * (MATRIX a){
		MATRIX ret;
		ret.init();
		FOR(n)for(int j=1;j<=n;j++){
			for(int k=1;k<=n;k++)ret.x[i][j]=ret.x[i][j]+x[i][k]*a.x[k][j];
			ret.x[i][j