【BZOJ 4128】 matrix

题目其实很简单,BSGS可以一眼看出来,但这就需要矩阵求逆,然而我并不会。。。
于是发现了一种BSGS的非求逆方法,借此介绍一下。

aixy=b(modP)
推出
aix=bay(modP)

因此我们只需要对每个 bay hash一下就好了,然后穷举 i <script type="math/tex" id="MathJax-Element-157">i</script>,验证是否存在即可。code:

#include<cstdio>
#include<cstring>
#include<cmath>
#include<map>
using namespace std;
int n,P;
struct hp{
    int num[71][71];
    bool operator < (const hp& a)const
      {
        int i,j;
        for (i=1;i<=n;++i)
          for (j=1;j<=n;++j)
            {
              if (num[i][j]<a.num[i][j]) return true;
              if (num[i][j]>a.num[i][j]) return false;
            }
        return false;
      }
    bool operator > (const hp& a)const
      {
        int i,j;
        for (i=1;i<=n;++i)
          for (j=1;j<=n;++j)
            {
              if (num[i][j]<a.num[i][j]) return false;
              if (num[i][j]>a.num[i][j]) return true;
            }
        return false;
      }
    bool operator == (const hp& a)const
      {
        int i,j;
        for (i=1;i<=n;++i)
          for (j=1;j<=n;++j)
            {
              if (num[i][j]<a.num[i][j]) return false;
              if (num[i][j]>a.num[i][j]) return false;
            }
        return true;
      }
}a,b;
hp mult(hp a,hp b)
{
    hp c;int i,j,k;
    for (i=1;i<=n;++i)
      for (j=1;j<=n;++j)
        {
          c.num[i][j]=0;
          for (k=1;k<=n;++k)
            c.num[i][j]=(c.num[i][j]+a.num[i][k]*b.num[k][j])%P;
        }
    return c;
}
int bsgs()
{
    int i,j,m;
    hp e,mi,xi;
    map<hp,int>x; 
    for (i=1;i<=n;++i)
      for (j=1;j<=n;++j)
        e.num[i][j]=b.num[i][j];
    m=(int)(sqrt(P)+0.5); x[e]=0;
    for (i=1;i<=m;++i)
      {
        e=mult(e,a);
        x[e]=i;
      }
    for (i=1;i<=n;++i) 
      for (j=1;j<=n;++j)
        {
          mi.num[i][j]=0;
          if (i==j)
            mi.num[i][j]=1;
        }
    for (i=1;i<=m;++i)
      mi=mult(mi,a);
    xi=mi;
    for (i=1;i<=m;++i)
      {
        if (x.count(mi)) return i*m-x[mi];
        mi=mult(mi,xi);     
      }
}
int main()
{
    int i,j;
    scanf("%d%d",&n,&P);
    for (i=1;i<=n;++i)
      for (j=1;j<=n;++j)
        scanf("%d",&a.num[i][j]);
    for (i=1;i<=n;++i)
      for (j=1;j<=n;++j)
        scanf("%d",&b.num[i][j]);
    printf("%d\n",bsgs());
}
  • 2
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值