洛谷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-1\equiv k*9(mod m)$$

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

然后直接套上BSGS的板子

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

 //minamoto
 #include<bits/stdc++.h>
 using namespace std;
 #define int __int128
 #define getc() (p1==p2&&(p2=(p1=buf)+fread(buf,1,1<<21,stdin),p1==p2)?EOF:*p1++)
 char buf[<<],*p1=buf,*p2=buf;
 inline int read(){
     #define num ch-'0'
     char ch;bool flag=;int res;
     while(!isdigit(ch=getc()))
     (ch=='-')&&(flag=true);
     for(res=num;isdigit(ch=getc());res=res*+num);
     (flag)&&(res=-res);
     #undef num
     return res;
 }
 inline void print(int x) {
     int sta[],top=;
     while (x) sta[++top]=x%,x/=;
     while (top) putchar(sta[top--]+'');
 }

 inline int BSGS(int a,int b,int p){
     map<int,int> mp;mp.clear();
     b%=p;int t=sqrt((double)p)+,tot=;
     for(int j=;j<t;++j){
         int val=b*tot%p;
         mp[val]=j,tot=tot*a%p;
     }
     if(!tot) return b==?:-;
     a=tot,tot=;
     for(int i=;i<=t;++i){
         int j=mp.find(tot)==mp.end()?-:mp[tot];
         if(j>=&&i*t-j>=) return i*t-j;
         tot=tot*a%p;
     }
     return -;
 }
 signed main(){
 //    freopen("testdata.in","r",stdin);
     int k=read(),m=read(),ans=BSGS(,(*k+)%m,m);
     print(ans);
     return ;
 }