loj 572 Misaka Network 与求和 —— min_25筛

题目：https://loj.ac/problem/572

推式子：https://www.cnblogs.com/cjoieryl/p/10150718.html

又学习了一下杜教筛hh；

原来 unsigned int 的输出是 %u 啊；

注意各处还是要用 (ll)，不要不小心都写成 (uint) 了；

然而递归版很慢...

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
using namespace std;
typedef unsigned int uint;
typedef long long ll;
int const xn=1e6+;
int n,m,K,pri[xn],cnt,w[xn],sqr;
uint prk[xn],h[xn],G[xn];
bool vis[xn];
uint pw(uint a,int b){uint ret=; for(;b;b>>=,a=a*a)if(b&)ret=ret*a; return ret;}
void init(int mx)
{
  for(int i=;i<=mx;i++)
    {
      if(!vis[i])pri[++cnt]=i,prk[cnt]=pw(i,K);
      for(int j=;j<=cnt&&(ll)i*pri[j]<=mx;j++)
    {
      vis[i*pri[j]]=;
      if(i%pri[j]==)break;
    }
    }
}
int Id(int x)
{
  if(x>sqr)return n/x;
  return m-x+;
}
uint F(int x,int y)
{
  if(pri[y]>x)return ;
  uint ret=;
  for(int i=y;i<=cnt&&(ll)pri[i]*pri[i]<=x;i++)//ll
    for(ll p0=pri[i];p0*pri[i]<=x;p0*=pri[i])
      ret+=F(x/p0,i+)+(uint)prk[i]*(h[Id(x/p0)]-i+);
  return ret;
}
uint S(int x)
{
  if(G[Id(x)]!=-)return G[Id(x)];
  uint ret=F(x,)+h[Id(x)];
  for(int i=,j;i<=x;i=j+)j=x/(x/i),ret-=(j-i+)*S(x/i);
  return G[Id(x)]=ret;
}
int main()
{
  scanf("%d%d",&n,&K); sqr=sqrt(n); init(sqr);
  for(int i=,j;i<=n;i=j+)
    {w[++m]=n/i; j=n/w[m]; h[m]=w[m]-;}
  for(int j=;j<=cnt;j++)
    for(int i=;i<=m&&(ll)pri[j]*pri[j]<=w[i];i++)
      h[i]=h[i]-h[Id(w[i]/pri[j])]+j-;//w[i]
  memset(G,-,sizeof G);
  uint ans=;
  for(int T=,nxt;T<=n;T=nxt+)
    {
      nxt=n/(n/T);
      ans+=(uint)(n/T)*(n/T)*(S(nxt)-S(T-));
    }
  printf("%u\n",ans);//
  return ;
}

递归版

于是写成了循环版，真的变快了^_^

代码如下：

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
using namespace std;
typedef unsigned int uint;
typedef long long ll;
int const xn=1e6+;
int n,m,K,pri[xn],cnt,w[xn],sqr;
uint prk[xn],h[xn],G[xn],f[xn];
bool vis[xn];
uint pw(uint a,int b){uint ret=; for(;b;b>>=,a=a*a)if(b&)ret=ret*a; return ret;}
void init(int mx)
{
  for(int i=;i<=mx;i++)
    {
      if(!vis[i])pri[++cnt]=i,prk[cnt]=pw(i,K);
      for(int j=;j<=cnt&&(ll)i*pri[j]<=mx;j++)
    {
      vis[i*pri[j]]=;
      if(i%pri[j]==)break;
    }
    }
}
int Id(int x)
{
  if(x>sqr)return n/x;
  return m-x+;
}
void getf()
{
  for(int i=cnt;i;i--)
    for(int j=;j<=m&&(ll)pri[i]*pri[i]<=w[j];j++)
      for(ll p0=pri[i];p0*pri[i]<=w[j];p0*=pri[i])
    f[j]+=f[Id(w[j]/p0)]+(uint)prk[i]*(h[Id(w[j]/p0)]-i+);
}
uint S(int x)
{
  if(G[Id(x)]!=-)return G[Id(x)];
  uint ret=f[Id(x)]+h[Id(x)];
  for(int i=,j;i<=x;i=j+)j=x/(x/i),ret-=(j-i+)*S(x/i);
  return G[Id(x)]=ret;
}
int main()
{
  scanf("%d%d",&n,&K); sqr=sqrt(n); init(sqr);
  for(int i=,j;i<=n;i=j+)
    {w[++m]=n/i; j=n/w[m]; h[m]=w[m]-;}
  for(int j=;j<=cnt;j++)
    for(int i=;i<=m&&(ll)pri[j]*pri[j]<=w[i];i++)
      h[i]=h[i]-h[Id(w[i]/pri[j])]+j-;//w[i]
  memset(G,-,sizeof G);
  uint ans=; getf();
  for(int T=,nxt;T<=n;T=nxt+)
    {
      nxt=n/(n/T);
      ans+=(uint)(n/T)*(n/T)*(S(nxt)-S(T-));
    }
  printf("%u\n",ans);//
  return ;
}