hihoCoder-1839 榶榶榶 数学

题面

题意:给你一个500000长度的数字,然后环形的让每位做头,例如123,就有123,231,312三个,然后问这n个数字的和S,S的最小非1因子是多少

题解:每个数字在每个位置都会有一次,如果说所有数字之和为A(1234的A就是10),那么原始长度为N的数,S就一定长成AAAAAAAAAAAAAA

也就是  S=A*10^0+A*10^1+A*10^2+...+A*10^(n-1)

10S=             A*10^1+A*10^2+...+A*10^(n-1)+A*10^n

9S=A*10^n-A

那其实S就是A*1111111(n个1),也就是找A和111111(n个1)因子中最小的那个 这个数当然也一定是质数

对于A,我们知道500000*9=4500000,也就是A最大这么多,可以暴力sqrt(A),找到他最小的因子

但是对于11111(n个1),他的最小因子就没有规律,通过打表可以发现,除了每3项有一个可以被3整除,还常有11,41,并时不时的出现他自身(也就是说这个数是质数)

但是题目说了保证答案是小于5*10^6 所以我们可以预处理这范围内的所有素数,然后挨个判断能否被1111(n个1)整除

那这样其实我们就可以直接对S进行操作,不需要拆分了,S=(A*10^n-A)/9,枚举质因子p 意思就是要满足(A*10^n-A)%(9p)==0就行了

 #include<bits/stdc++.h>
 using namespace std;
 char c;
 int x,tot,ans;
 int f[];
 int pow(int a,int b,int mod)
 {
     if (b==) return a;
     long long t=pow(a,b/,mod);
     t=t*t%mod;
     if (b%) t=t*a%mod;
     return t;
 }
 int check(int n,int s,int m)
 {
     int res=(pow(,n,*m)+*m-)%(*m)/;
     return 1LL*res*s%m==;
 }
 int main()
 {
     for(int i=;i<;i++)
         if (!f[i])
             for(int j=i+i;j<;j+=i) f[j]=;
     while (scanf("%c",&c)!=EOF)
     {
         if (c>='' && c<='')
         {
             tot++;
             x+=c-'';
         }else break;
     }
     for(int i=;i<=;i++)
         if (!f[i] && check(tot,x,i))
         {
             printf("%d",i);
             return ;
         }
 }