HDU4542 小明系列故事——未知剩余系

大赞的数论题;

大致思路:

对于TYPE=1的情况,认为    X 中有  X-K个约数,求最小的X,X-K>0 那么化为B+K的约数为B,

我们知道(B+K)的约数<=2*SQRT(B+K);这个应该知道,

再化简一下:B*B<=4(B+K),可以直接枚举B,有人担心会TLE,我们再证明一下:

化简得 :B*B-4*B+4<=4*K+4;然后K<=47777,所以可以枚举了。。

对于TYPE=0的情况, 我们可以预先求出有X个数的最小数,因为根据求约数个数的方程N=a[1]^x1+a[2]^x2.....;

num=(1+x1)*(1+x2)*。。。。;

这里我们可以DFS搜索出来。具体看看代码应该能明白

#include<cstdio>
#include<algorithm>
#include<cmath>
#include<string.h>
typedef long long ll;
using namespace std;
#define maxn 200000
const long long inf=(1ll<<62)+1;
int prime[maxn];
ll a[maxn*10];
int b[maxn*2],t; void pri()
{
for (int i=2;i*i<maxn;i++)
if (!b[i])
for (int j=i*2;j<maxn;j+=i)
b[j]=1;
for (int i=2;i<maxn;i++)
if (!b[i]) prime[++t]=i;
} int get(ll x){
int ans=1;
ll tmp=x;
for (int i=1;prime[i]<=tmp/prime[i];i++)
{
ll t=1;
while (tmp%prime[i]==0)
{
t++;
tmp/=prime[i];
}
if (t>1) ans*=t;
}
if (tmp!=1) ans*=2;
return ans;
} void dfs(int i,ll x,int n)//DFS部分
{
if (n>47777) return;
if (x<inf&&(a[n]==0||a[n]>x)) a[n]=x;//类似DP的思想
for (int j=1;j<=62;j++)
{
if (inf/prime[i]<x) break;//防止溢出,
x*=prime[i];
if (x>=inf) break;
dfs(i+1,x,n*(j+1));
}
} void solve1(int x)
{
if (a[x]!=0) printf("%I64d\n",a[x]);
else printf("INF\n");
}
void solve(ll x)
{
ll num=1;
while (num*num<=(num+x)*4)
{
if (num==get(num+x)) {printf("%I64d\n",num+x);return; }
num++;
}
printf("Illegal\n");
} int main()
{
int T;
scanf("%d",&T);
pri();
dfs(1,1,1);
for (int o=1;o<=T;o++)
{
printf("Case %d: ",o);
int type,k;
scanf("%d%d",&type,&k);
if (type) solve(k);
else solve1(k);
}
return 0;
}

  

上一篇:springmvc的几点见解


下一篇:java 调用微信截图工具