LightOJ1007 Mathematically Hard

LightOJ1007
题目的意思很容易理解,求一个区间[a,b]内的欧拉函数平方之和
但是题目的数据限很大 a,b都是1e12级。按照一般思路写好程序后,实验以下从2到10000发现long long也会超范围,所以用unsigned long long。
此外还需要预处理前缀和,不然就会超时。
预处理前缀和时不能够再开数组,不然就会爆内存,WA点多多.
AC代码:

#include <bits/stdc++.h>
typedef unsigned long long ull;
const int maxn=5000000+10;
using namespace std;
ull euler[maxn],a,b;
void geteuler_table(){          //欧拉函数打表
    memset(euler,0,sizeof(euler));
    euler[1]=0;
    for(int i=2;i<=maxn;i++){
            if(euler[i]==0)
            for(int j=i;j<=maxn;j+=i){
                if(euler[j]==0)
                    euler[j]=j;
                euler[j]=euler[j]/i*(i-1);
            }
    }
    //cout<<euler[2]<<endl;
    for(int i=2;i<=maxn;i++)        //预处理前缀和
        euler[i]=euler[i-1]+euler[i]*euler[i];
}
int main(){
    //freopen("input.txt","r",stdin);
    int T,cnt=0;
    scanf("%d",&T);
    geteuler_table();
    while(T--){
        ull ans;
        scanf("%lld%lld",&a,&b);
        printf("Case %d: ",++cnt);
        ans=euler[b]-euler[a-1];
        printf("%llu\n",ans);
    }
}
  • 0
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 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、付费专栏及课程。

余额充值