D. Soldier and Number Game/牛牛的“质因数”(欧拉筛处理最小/最大质因数)

 https://ac.nowcoder.com/acm/contest/9982/I

https://codeforces.com/problemset/problem/546/D

思路:总的思路就是要o(n)拿出每个数的最大/最小质因数。然后就可以logn处理。

欧拉筛本身的处理就是可以把v[i]存成i的最小质因数来筛,详见F. 1.小W 的质数(prime)[欧拉筛再理解]

如果我要处理出每个数的最大质因数呢?类似一个dp的思路

Maxp[x]:表示x的最大质因子

枚举的i是x的最大因子

Maxp[i*primes[j]]=max(primes[j],Maxp[i]);不断更新能整除的小的质因子得出最大的/

	n=read();
	for(int i=2;i<=10000000;++i){
		if(!vis[i]){p[++cnt]=i;Maxp[i]=i;}
		for(int j=1;p[j]*i<=10000000;++j){
			vis[p[j]*i]=1;//线性筛
			Maxp[p[j]*i]=max(p[j],Maxp[i]);//求最大的质因子
			if(i%p[j]==0)break;
		}
	}
	for(int i=1;i<=n;++i){
		a[i]=read();
		int tmp=a[i];
		while(tmp!=1){
			int t=Maxp[tmp];
			while(tmp%t==0)tmp/=t,printf("%d ",t);
		}
		puts("");
	}

#include<iostream>
#include<vector>
#include<queue>
#include<cstring>
#include<cmath>
#include<map>
#include<set>
#include<cstdio>
#include<algorithm>
#define debug(a) cout<<#a<<"="<<a<<endl;
using namespace std;
const int maxn=4e6+100;
typedef int LL;
const int mod=1e9+7;
inline LL read(){LL x=0,f=1;char ch=getchar();	while (!isdigit(ch)){if (ch=='-') f=-1;ch=getchar();}while (isdigit(ch)){x=x*10+ch-48;ch=getchar();}
return x*f;}
LL primes[maxn];
LL v[maxn];
LL m=0;
LL ksm(LL a,LL b){
    LL ans=1;
    while(b){
        if(b&1)ans=ans*a%mod;
        a=a*a%mod;
        b/=2;
    }
    return ans%mod;
}
void getprimes(int n)
{
	for(LL i=2;i<=n;i++)
	{
		if(v[i]==0) {v[i]=i;primes[++m]=i;}
		for(LL j=1;j<=m;j++)
		{
			if(primes[j]>v[i]||primes[j]>n/i)  break;
			v[i*primes[j]]=primes[j];
		}
	}
}
LL get(LL x){
    LL sum=0;
    while(x>0) x/=10,sum++;
    return  sum;
}
LL solve(LL x){
    long long sum=0;
    if(x==1) return 0;
    LL tmp=x;
    while(tmp!=1){
        LL t=v[tmp];
        while(tmp%t==0){
            tmp/=t;
            sum=(sum*ksm(10,get(t))%mod+t%mod)%mod;
        }
    }
    if(tmp>1) sum=(sum*ksm(10,get(tmp))%mod+tmp%mod)%mod;
    return sum%mod;
}
int main(void)
{
  getprimes(4000000);

  LL n;n=read();
  long long ans=0;
  for(LL i=1;i<=n;i++){
    ans=(ans%mod+solve(i)%mod)%mod;
  }
  printf("%lld\n",ans);
return 0;
}

 

  • 0
    点赞
  • 3
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值