題目的意思大概就是問是否存在一串全是8的數字是L的倍數

直接想沒有什么想法，要想到用簡潔的形式將這個數字表示出來，對于每一位都是8的數字我們可以用
X=8*(10^k-1)/9的形式表示出來，那么題目的意思就是求X使L|X，我們先處理一下8和L，即除去他們的最大公約數，然后就是L’|(10^k-1)/9，即就是10^k-1|9L’我們用L’'表示9L’

問題就轉化成了要求10^k-1%L’’==0，10^k=1（mod L’’）(其實找的是在模L’'剩余系中10的階)

如果10和L’'不互質那么10沒有階

由歐拉定理我們得知，10^f(L’’)=1（mod L’’），因此階一定是f(L’’)的因子，其中f(n)代表的是n的歐拉函數，因此我們從小到大查找歐拉函數值的所有因子直到找到階

在這種思路下有如下代碼

#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<algorithm>
#include<iostream>
#include<cmath>
#include<ctime>
#include<climits>
#include<queue>
#include<vector>
#include<set>
#include<map>
using namespace std;typedef long long ll;
const int INF=0x3f3f3f3f;
const int MAXN=1e5+5;ll L;ll gcd(ll a,ll b)
{return b==0?a:gcd(b,a%b);
}ll phi(ll x)
{ll ret=1;for(ll i=2;i*i<=x;i++){if(x%i==0){ret*=(i-1);x/=i;while(x%i==0){ret*=i;x/=i;}}if(x==1) break;}if(x>1) ret*=(x-1);return ret;
}ll mult(ll a,ll b,ll t)
{a%=t; b%=t;ll ret=0;while(b){if(b&1){ret+=a; if(ret>t) ret-=t;}a<<=1; b>>=1;if(a>t) a-=t;}return ret;
}ll quick_pow(ll a,ll b,ll t)
{ll ret=1;a%=t;while(b){if(b&1) ret=mult(ret,a,t);a=mult(a,a,t);b>>=1;}return ret;
}int main()
{int Case=0;while(~scanf("%lld",&L) && L){++Case;ll t=gcd(8,L);L=L/t*9;ll ans;if(L%2==0 || L%5==0){ans=0;}else{t=phi(L);ans=-1;//printf("t=%lld\n",t);ll tt=sqrt(t)+1;for(ll i=1;i<=tt;i++){if(t%i==0 && quick_pow(10,i,L)==1){ans=i;break;}}if(ans<0)for(ll i=tt;i>0;i--){if(t%i==0 && quick_pow(10,t/i,L)==1){ans=t/i;break;}}}printf("Case %d: %lld\n",Case,ans);}return 0;
}

這里求歐拉函數值的用了p為素數，f(p^k)=(k-1)*p^k-1。然后在遍歷歐拉函數因子的時候用到一個技巧：任何一個數的因子都可以和另一個因子相乘得到這個數，這兩個因子中一個大一個小，小的一定小于等于sqrt(n)，因此我們在查找因子的時候不要遍歷1—n，而是先遍歷1—sqrt(n)查找較小的因子，如果沒有找到， 再從sqrt(n)—1,查找相對應的較大的因子，這樣就將原來O(n)的復雜度降低為O(logn)的復雜度,十分巧妙

這里放一個優化版本的

#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<algorithm>
#include<iostream>
#include<cmath>
#include<ctime>
#include<climits>
#include<queue>
#include<vector>
#include<set>
#include<map>
using namespace std;typedef long long ll;
const int INF=0x3f3f3f3f;
const int MAXN=1e5+5;ll L;ll gcd(ll a,ll b)
{return b==0?a:gcd(b,a%b);
}ll phi(ll x)
{ll ret=x;for(ll i=2;i*i<=x;i++){if(x%i==0){ret=ret/i*(i-1);x/=i;while(x%i==0){x/=i;}}if(x==1) break;}if(x>1) ret=ret/x*(x-1);return ret;
}ll mult(ll x, ll y, ll p) 
{long double d=1;d=d*x/p*y;return ((x*y-((ll)d)*p)%p+p)%p;
}ll quick_pow(ll a,ll b,ll t)
{ll ret=1;a%=t;while(b){if(b&1) ret=mult(ret,a,t);a=mult(a,a,t);b>>=1;}return ret;
}int main()
{int Case=0;while(~scanf("%lld",&L) && L){++Case;ll t=gcd(8,L);L=L/t*9;ll ans;if(L%2==0 || L%5==0){ans=0;}else{t=phi(L);ans=-1;//printf("t=%lld\n",t);ll tt=sqrt(t)+1;for(ll i=1;i<=tt;i++){if(t%i==0 && quick_pow(10,i,L)==1){ans=i;break;}}if(ans<0)for(ll i=tt;i>0;i--){if(t%i==0 && quick_pow(10,t/i,L)==1){ans=t/i;break;}}}printf("Case %d: %lld\n",Case,ans);}return 0;
}

主要優化了快速冪中long long 乘法部分和歐拉函數值的求值.

這里歐拉函數值的求法用到歐拉函數容斥原理求法.即
$f(n)=n?(p1?1)/p1?(p2?1)/p2?.....?(pk?1)/pk$這里pi為n的素因子