zoukankan      html  css  js  c++  java
  • 洛谷P4884 多少个1?(BSGS)

    传送门

    模数好大……__int128好麻烦……而且BSGS第一次写有点写蒙了……

    $11...1(N个1)equiv k(mod m)$很难算,那么考虑转化一下

    先把$11...1(N个1)$写成$frac{10^n-1}{9}$

    则$$frac{10^n-1}{9}equiv k(mod m)$$

    $$10^n-1equiv k*9(mod m)$$

    $$10^nequiv k*9+1(mod m)$$

    然后直接套上BSGS的板子

    然后因为模数是long long级别的,所以要么用龟速乘,要么像我一样懒得只会用__int128了(记得手打输入输出)

     1 //minamoto
     2 #include<bits/stdc++.h>
     3 using namespace std;
     4 #define int __int128
     5 #define getc() (p1==p2&&(p2=(p1=buf)+fread(buf,1,1<<21,stdin),p1==p2)?EOF:*p1++)
     6 char buf[1<<21],*p1=buf,*p2=buf;
     7 inline int read(){
     8     #define num ch-'0'
     9     char ch;bool flag=0;int res;
    10     while(!isdigit(ch=getc()))
    11     (ch=='-')&&(flag=true);
    12     for(res=num;isdigit(ch=getc());res=res*10+num);
    13     (flag)&&(res=-res);
    14     #undef num
    15     return res;
    16 }
    17 inline void print(int x) {
    18     int sta[30],top=0;
    19     while (x) sta[++top]=x%10,x/=10;
    20     while (top) putchar(sta[top--]+'0');
    21 }
    22 
    23 
    24 inline int BSGS(int a,int b,int p){
    25     map<int,int> mp;mp.clear();
    26     b%=p;int t=sqrt((double)p)+1,tot=1;
    27     for(int j=0;j<t;++j){
    28         int val=b*tot%p;
    29         mp[val]=j,tot=tot*a%p;
    30     }
    31     if(!tot) return b==0?1:-1;
    32     a=tot,tot=1;
    33     for(int i=0;i<=t;++i){
    34         int j=mp.find(tot)==mp.end()?-1:mp[tot];
    35         if(j>=0&&i*t-j>=0) return i*t-j;
    36         tot=tot*a%p;
    37     }
    38     return -1;
    39 }
    40 signed main(){
    41 //    freopen("testdata.in","r",stdin);
    42     int k=read(),m=read(),ans=BSGS(10,(9*k+1)%m,m);
    43     print(ans);
    44     return 0;
    45 }
  • 相关阅读:
    C# 冒泡排序
    C# IO流 File.Exists,Directory.Exists, File.Create,Directory.CreateDirectory
    UGUI 哪些显示在前方的问题
    UGUI Image血条或者进度条效果
    Unity 切换场景的时候让某个游戏对象不消失
    C# string型的转换成float型的
    Json 解析Json
    logging模块

    模块和包
  • 原文地址:https://www.cnblogs.com/bztMinamoto/p/9740439.html
Copyright © 2011-2022 走看看