Golang版 YY的GCD(莫比乌斯反演)

学会一种语言最好的办法就是拿它写一道算法题{

鉴于sort和stl实现不了,就选取了莫比乌斯反演的题来实现

golang语法还是挺臭屁的,初学者可能会引起强烈不适,但是大的上面难度还是很低的,这一番操作go也熟悉了不少

package main
 
import "fmt"
 
func Min(x, y int) int {
    if x < y {
        return x
    }
    return y
}
func Max(x, y int) int{
    if x > y {
        return x
    }
    return y
}
const limit int = 10000005;
var n int;
var num[10000005] int;
var phi[10000005] int;
var miu[10000005] int;
var sum[10000005] int64;
var prime[10000005] int;
var tot int;
func Memset(){
    for i := 0 ; i < 10000005; i++{
        num[i] = 1;
        phi[i] = 0;
        sum[i] = 0;
        miu[i] = 0;
    }
}//全清空
func Ssss(){
    Memset();
    tot = 0;
    num[1] = 0;
    num[0] = 0;
    miu[1] = 1;
    for i := 2 ;i < 10000005 ; i++{
        if num[i] != 0{
            tot++;
            prime[tot] =i;
            miu[i] = -1;
        }
        for j := 1; j <= tot && prime[j] * i < 10000005 ;j++{
            num[prime[j] * i] = 0;
            if i % prime[j] == 0{
                miu[i * prime[j]] = 0;
                break;
            }else{
                miu[i * prime[j]] = -1 * miu[i];//莫比乌斯函数
            }
        }
    }
    for i := 1 ; i <= tot; i++{
        for j := prime[i] ; j < 10000005 ; j += prime[i]{
            sum[j] += int64(miu[j / prime[i]]);//几个平方因子
        }
    }
    for i := 1; i < 10000005 ;i++{
        sum[i] += sum[i-1];//前缀和
    }
}
var m int;
func Calc1()int64{
    var ans int64 = 0;
    var r int;
    for l := 1 ; l <= Min(n,m);l = r + 1 {
        //值域分块
        var t = n;
        var s = m;
        r = Min(t / (t / l), s / (s / l));
        ans = ans + (int64(sum[r] - sum[l - 1]) * int64(int64(t / l) * int64(s / l)));
        //fmt.Printf("%d\n",(sum[r] - sum[l - 1]) * (int64(t / l) * int64(s / l)));
    }
    return ans;
}
func main() {
    Ssss();
    var kase int;
    fmt.Scanf("%d",&kase);
    for i := 1; i<=kase;i++{
        fmt.Scanf("%d%d",&n,&m);
        fmt.Println(Calc1());
    }
}
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

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

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

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值