题意:
若一个大于1的整数M的质因数分解有k项,其最大的质因子为Ak,并且满足Ak^K<=N,Ak<128,我们就称整数M为N-伪光滑数。现在给出N,求所有整数中,第K大的N-伪光滑数。
题解:
好题,只能%:beginend
话说下次遇到这种第k大是不是也可以这么搞啊
code:
#include<queue>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<iostream>
#include<algorithm>
#define LL long long
using namespace std;
int prime[50],pr=0;
bool v[130];
struct trnode{
int lc,rc,d;
LL c,t;
}tr[17000005];int tot=0;
int g[50][110],f[50][110];
LL n;int K;
void pre()
{
memset(v,true,sizeof(v));
for(int i=2;i<128;i++)
if(v[i])
{
prime[++pr]=i;
for (int j=i*2;j<128;j+=i) v[j]=false;
}
}
int new_node(int x,LL t)
{
tot++;tr[tot]=tr[x];
tr[tot].c*=t;tr[tot].t*=t;
return tot;
}
void pushdown(int x)
{
if(tr[x].t==1) return;
LL t=tr[x].t;tr[x].t=1;
if(tr[x].lc) tr[x].lc=new_node(tr[x].lc,t);
if(tr[x].rc) tr[x].rc=new_node(tr[x].rc,t);
}
int merge(int x,int y)
{
if(!x||!y) return x+y;
pushdown(x);pushdown(y);
if(tr[x].c<tr[y].c) swap(x,y);
int v=new_node(x,1);tr[v].rc=merge(tr[v].rc,y);
if(tr[tr[v].lc].d<tr[tr[v].rc].d) swap(tr[v].lc,tr[v].rc);
tr[v].d=tr[tr[v].rc].d+1;
return v;
}
priority_queue<pair<LL,pair<int,int> > > q;
void make_heap()
{
pre();
g[0][0]=++tot;tr[tot].c=tr[tot].t=1;
for(int i=1;i<=pr;i++)
{
g[i][0]=1;
for(LL j=1,p=prime[i];p>0&&p<=n;j++,p*=prime[i])
{
for(LL k=1,pri=prime[i];k<=j;k++,pri*=prime[i]) f[i][j]=merge(f[i][j],new_node(g[i-1][j-k],pri));
g[i][j]=merge(g[i-1][j],f[i][j]);
q.push(make_pair(tr[f[i][j]].c,make_pair(i,j)));
}
}
}
void solve()
{
LL ans;
while(K--)
{
ans=q.top().first;int i=q.top().second.first,j=q.top().second.second;q.pop();
pushdown(f[i][j]);f[i][j]=merge(tr[f[i][j]].lc,tr[f[i][j]].rc);
q.push(make_pair(tr[f[i][j]].c,make_pair(i,j)));
}
printf("%lld",ans);
}
int main()
{
scanf("%lld %d",&n,&K);
make_heap();solve();
}