Mysterious Bacteria

17 篇文章 0 订阅

这里写图片描述这里写图片描述

题目大意:

T组样例,给你一个x,求出满足x=b^p的最大幂p。

解题思路:

使用暴力用对数函数(log(n)/log(m)即以m为底n的对数)枚举必然超时,正确的做法是用唯一素数分解定理分解质因数然后求出所有次幂的最大公约数。如果n是负数需要特殊处理一下次幂,奇数次幂不需要处理,偶数次幂需要不断降幂(即除2)得到最终的奇数次幂。
代码如下:

#include<iostream>
#include<cstdio>
#include<fstream>
#include<set>
#include<cmath>
#include<cstring>
#include<string>
#include<map>
#include<vector>
#include<iomanip>
#include<cstdlib>
#include<list>
#include<queue>
#include<stack>
#include<algorithm>
#define inf 0x3f3f3f3f
#define MOD 1000000007
#define mem0(a) memset(a,0,sizeof(a))
#define mem1(a) memset(a,-1,sizeof(a))
#define meminf(a) memset(a,inf,sizeof(a))
#define MAXN 1e5+5
//vector ::iterator it;
typedef long long ll;
typedef unsigned long long ull;
using namespace std;
ll prime[1000010];
bool book[1000010];
int GCD(int a,int b)//求解这些幂的最大公约数
{
  if(b==0)return a;
  else GCD(b,a%b);
}
int getprime()
{
  int cnt=0;
  for(int i=2;i<=MAXN;i++)
  {
    if(!book[i])prime[cnt++]=i;
    for(int j=0;j<cnt,i*prime[j]<=MAXN;j++)
    {
      book[i*prime[j]]=1;
      if(i%prime[j]==0)break;
    }
  }
  return cnt;
}
ll solve(ll n,int len)//唯一素数分解定理
{
  int cur=0,g;
  bool flag=true;
  while(cur<len&&n>1)
  {
    if(n%prime[cur]==0)
    {
      int mi=0;
      while(n%prime[cur]==0)
      {
        mi++;
        n/=prime[cur];
      }
      if(flag)
      {
        flag=false;
        g=mi;
      }
      g=GCD(g,mi);
    }
    cur++;
  }
  if(n!=1)g=1;
  return g;
}
int main()
{
  std::ios::sync_with_stdio(false);
  cin.tie(0);
  //freopen("test.txt","r",stdin);
//  freopen("output0.txt","w",stdout);
  ll T,cas=1;
  int len=getprime();
  cin>>T;
  while(T--)
  {
    ll n,ma=0,result;
    bool flag=true;
    cin>>n;
    if(n<0)
    {
      flag=false;
      n=-n;
    }
    result=solve(n,len);
    if(!flag&&result%2==0)//负数且幂是偶数的要特殊处理一下
    {
      while(result%2==0)//将幂分解为奇数
      {
        result/=2;
      }
    }
    if(n==1||n==-1)result=1;
    cout<<"Case "<<cas++<<": "<<result<<endl;
  }
  return 0;
}
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值