Lightoj1220三种解法(唯一分解,二分答案,pow开根)

Lightoj1220
题意:给你一个数n,要你求这个数可以是一个几次方数。
比如:n=64, n = 2 6 n=2^{6} n=26,所以答案为6。
大多数人给的题解都是下面这种唯一分解的。把我坑惨了QAQ
为什么呢?我在比赛中碰到一个一样的题,但数据范围不一样n<=1e18,t<=100000
我比赛的时候就不停的往这个唯一分解这个方向想答案,一直想不出,因为你素数筛打表开不出1e9*int的空间,就算有这么大的空间,预处理也超时了,如果不预处理,那么也超时。。。。想了半天没做出来,比赛结果很遗憾TAT
如果比赛有重现赛的话我就把题目链接发出来。

  1. 唯一分解
//12ms过
//当数据量大,时这种做法就做不了了,时间空间都过不了。所以推荐用二分或者pow开根做
#include <iostream>
#include <cstdio>
#include <algorithm>
#include <cstring>
#include <cmath>
using namespace std;
typedef long long ll;

const int N=1e6+5;
int p[N],k=0;
bool vis[N];

void prime(int n){
	memset(vis,0,sizeof(vis));
	int m=(int)sqrt(n+0.5);
	for(int i=2;i<=m;i++)if(!vis[i])
		for(int j=i*i;j<=n;j+=i)vis[j]=1;	
	for(int i=2;i<=n;i++)if(!vis[i])p[k++]=i;
}

ll cal(ll x){//算x的因子个数                                                                                                                                                                                                                                                                                                                                                                                                                                       
	ll cnt=0;
	bool sign=0;
	if(x<0)x=-x,sign=1;
	for(ll i=0;p[i]*p[i]<=x&&i<k;i++){
		if(x%p[i]==0){
			ll t=0;
			while(x%p[i]==0){x/=p[i],t++;}
			if(!cnt)cnt=t;
			else cnt=__gcd(cnt,t);
		} 
	}
	if(x>1)cnt=1;//别写return cnt 
	if(sign)while(cnt%2==0)cnt/=2;
	return cnt;
} 

int main(){
	prime(N);
	int t,kas=0;
	ll n;
	scanf("%d",&t);
	while(t--){
		scanf("%lld",&n);
		printf("Case %d: %lld\n",++kas,cal(n));	
	} 
}
  1. 二分解法:
#include <iostream>
#include <cstdio>
#include <algorithm>
#include <cstring>
#include <cmath>
using namespace std;
typedef long long ll;
// 1ms过 
ll n;
int judge(ll x,int y){
	ll temp=1;
	for(int i=0;i<y;i++){
		temp*=x;
		if(temp>n)return 1;
	}
	if(temp<n)return -1;
	return 0;//相等 
} 

int main(){
	int t,kas=0;
	scanf("%d",&t);
	while(t--){
		int ans=1;
		scanf("%lld",&n);
		bool f=1,flag=0;
		if(n<0)n=-n,flag=1;
		for(int i=33;i>=2&&f;i--){//从高次幂开始枚举,找到就退出 
			ll l=1,r=n;
			while(l<=r){
				ll mid=r+l>>1;	
				int temp=judge(mid,i);
				if(temp==1)r=mid-1;	
				else if(temp==-1)l=mid+1;
				else {
					ans=i,f=0;
					break;
				}	
			}	
		}
		if(flag)while(ans%2==0)ans/=2;
		printf("Case %d: %lld\n",++kas,ans);	
	} 
}

  1. pow开根
    可以用pow(n,1/x)去开方然后判断,开方结果是否是整数,但是注意结果会产生浮点误差会导致你结果错误,不过误差不会很大,所以我们需要将结果在附近找,然后回代验算一下就好了。
#include<bits/stdc++.h>
using namespace std;
typedef long long ll;

int main(){
	ios::sync_with_stdio(false);
	int T,k=1;
	cin>>T;
	while(T--){
		ll n;
		cin>>n;
		cout<<"Case "<<k++<<": ";
		bool f=0;
		if(n<0)f=1,n=-n;	
		int ans=0;
		for(int i=32;i>=1;i--){ 
			if(!(i&1)&&f)continue;
			ll temp=(ll)pow(n,1.0/(double)i);
			for(int j=-3;j<=3;j++)
				if(n==(ll)pow(temp+j,i))ans=i;
			if(ans){
				cout<<ans<<endl;
				break;
			}
		}	
	}	
} 
  • 1
    点赞
  • 2
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
Sigma函数是指一个数字的所有因子之和。给定一个数字n,需要求出有多少个数字的Sigma函数是偶数。\[2\] 为了解决这个问题,可以先筛选出n范围内的素数(范围在10^6即可),然后对n进行素因子分解。对于每个因子,如果它的Sigma函数中连乘的每一项都是偶数,那么整个Sigma函数就是偶数。具体实现中,可以判断每个因子的平方根是否为偶数,如果是偶数,则减去(平方根+1)/2。\[1\] 另外,还可以使用O(1)的做法来解决这个问题。根据观察,所有的完全平方数及其两倍的值都会导致Sigma函数为偶数。因此,可以直接计算n的平方根,然后减去(平方根+1)/2即可得到结果。\[3\] #### 引用[.reference_title] - *1* [Sigma Function](https://blog.csdn.net/PNAN222/article/details/50938232)[target="_blank" data-report-click={"spm":"1018.2226.3001.9630","extra":{"utm_source":"vip_chatgpt_common_search_pc_result","utm_medium":"distribute.pc_search_result.none-task-cask-2~all~insert_cask~default-1-null.142^v91^control_2,239^v3^insert_chatgpt"}} ] [.reference_item] - *2* *3* [【LightOJ1336】Sigma Function(数论)](https://blog.csdn.net/qq_30974369/article/details/79009498)[target="_blank" data-report-click={"spm":"1018.2226.3001.9630","extra":{"utm_source":"vip_chatgpt_common_search_pc_result","utm_medium":"distribute.pc_search_result.none-task-cask-2~all~insert_cask~default-1-null.142^v91^control_2,239^v3^insert_chatgpt"}} ] [.reference_item] [ .reference_list ]

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

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

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值