首页 > 其他分享 >整除相关

整除相关

时间:2022-09-03 12:47:26浏览次数:52  
标签:phi gcd int ans 相关 整除 lcm operatorname

gcd 相关

首先要熟悉几个转化:
\(\gcd (a,b) = \gcd(a, a + qb)\)
\(\gcd (a, b) = \gcd(a, a - b)\)
\(\gcd (a, b) = \gcd(a, a \% b)\)

这几个转化很常见。
还有:
\(\operatorname{lcm}(a, \gcd(b, c)) = \gcd(\operatorname{lcm}(a, b), \operatorname{lcm}(a, c))\)
\(\gcd(a, \operatorname{lcm}(b, c)) = \operatorname{lcm}(\gcd(a, b), \gcd(a, c))\)

CF1717E

求 \(\sum \limits_{a+b+c=n} \operatorname{lcm}(c, \gcd(a, b))\)。

暴力枚举 \(c\),然后 \(a+b=n-c\)。利用第一个转化发现原式 \(= \sum \limits _{a \in [1, n - c - 1]} \operatorname{lcm}(c, \gcd(a, n - c))\)。

这时候要用到一个经典的例题(本题重中之重):计算 \(a \in [1, n-c-1], \gcd(a, n-c)\) 为各个数的 \(a\) 的方案数。

首先,\(\gcd(a, n-c)\) 肯定是 \(n-c\) 的因子。那么我们枚举 \(\gcd\),假设为 \(g\),那么如果 \(\gcd(a, n-c) = g\),\(a\) 要满足:\(a | g, (a / g) ⊥ ((n-c) / g)\)。

刚好,我们把 \(1 \sim n-c-1\) 中所有整除 \(g\) 的数除了 \(g\) 之后的取值映射到 \(1 \sim (n-c) / g\) 中,并且需要与 \((n-c) / g\) 互质。方案数显然是 \(\phi((n-c)/g)\)。

需要注意,当 \(g = n-c\) 时算出来的 \(\phi\) 为 \(1\),但这个 \(1\) 其实是 \(a = n - c\) 贡献的,不符合要求,不能计算。

#include <bits/stdc++.h>
using namespace std;
#define int long long
#define f(i, a, b) for (int i = (a); i <= (b); i++)
#define cl(i, n) i.clear(), i.resize(n);
#define endl '\n'
typedef long long ll;
typedef unsigned long long ull;
typedef pair<int, int> pii;
const int inf = 1e9;
const int mod = 1e9 + 7;
int qpow(int x, int k)
{
    int ans = 1;
    while (k)
    {
        if (k & 1)
            ans = ans * x % mod;
        x = x * x % mod;
        k >>= 1;
    }
    return ans;
}
int lcm(int x, int y)
{
    return x * y / __gcd(x, y);
}
int n;
int phi[100010];
int calc(int x)
{
    int ans = x;
    for (int i = 2; i * i <= x; i++)
    {
        if (x % i == 0)
        {
            ans = ans * (i - 1) / i;
            while (x % i == 0)
            {
                x /= i;
            }
        }
    }
    if (x > 1)
        ans = ans * (x - 1) / x;
    return ans;
}
signed main()
{
    ios::sync_with_stdio(0);
    cin.tie(NULL);
    cout.tie(NULL);
    time_t start = clock();
    // think twice,code once.
    // think once,debug forever.
    int ans = 0;
    cin >> n;
    f(i, 1, n) phi[i] = calc(i);
    f(c, 1, n - 2)
    {
        int k = n - c;
        for (int i = 1; i * i <= k; i++)
        {
            if (k % i == 0)
            {
                ans += phi[k / i] * lcm(c, i);
                ans %= mod;

                if (i * i != k && i != 1)
                {
                    ans += phi[i] * lcm(c, k / i);
                    ans %= mod;
                }
            }
        }
    }
    cout << ans << endl;
    time_t finish = clock();
    // cout << "time used:" << (finish-start) * 1.0 / CLOCKS_PER_SEC <<"s"<< endl;
    return 0;
}

标签:phi,gcd,int,ans,相关,整除,lcm,operatorname
From: https://www.cnblogs.com/Zeardoe/p/16652355.html

相关文章

  • http相关协议2
    协议怎么诞生的呢在tcp协议簇的基础之上1991年英国人tom研究了一套改变世界数据交互的创新 1991年tom发明第五个创新hhtp协议超文本传输......
  • 会话与认证介绍以及相关漏洞案例分析
    很开心能在这边跟大家分享有关于会话跟认证相关的内容以及跟认证相关的漏洞案例分析,在这个过程中过大家有什么疑惑或者有什么我讲得不对的地方,欢迎大家及时提出,我们一起交......
  • 【基础整理】Mapping representation 机器人所用地图种类及相关介绍
    参考与前言本文主要介绍建图Mapping方面的一些基础知识介绍与相关下游任务使用涉及知识较为基础,SLAM大佬们可以提前退出了主要针对应用为移动机器人与物流无人驾驶......
  • osi相关学习
    1.分层的优点2.七层模型3.TCP/IP参考模型4.理想的网络设计5.数据封装的方法 1.分层的优点层次划分的方法网络的每层应当具有相对独立的功能(便于排错)这个功能用......
  • 关于『数论』:整除理论
    「一本书上每多一个公式,就会死掉减少一半读者。」——霍金「这上面每多一个公式,我就会失明一次。」——JQ序言  数论真有趣(弥天大雾。  (投放Myblogscsdn......
  • 05_Linux基础-NGINX编译安装^判断是否启动^修改端口^启动停止重启^相关路径^中文乱码-
    05_Linux基础-NGINX编译安装判断是否启动修改端口启动停止重启相关路径^中文乱码-Windows、Linux文件传输博客......
  • 实用工具相关
    mindmanager破解与使用快捷键大全一、格式快捷键①Ctrl+B将字体加粗②Ctrl+U对选定文本加下划线③Ctrl+I使选定文本变成斜体④Ctrl+Shift+.增加字体大小⑤Ctrl+Sh......
  • 认识Node.js以及相关说明
    首先我们先要认识一下,Node是用JavaScript来编写的。问:JavaScript时怎么组成的?答:JavaScript主要是由两块核心组成,第一个核心是js的核心语法,该语法也被称为ECMAScript,语法......
  • 论文笔记 — CheXGAT:基于胸部 X 光图像诊断胸部疾病的疾病相关感知网络
    论文笔记—车行:用于从胸部X射线图像诊断胸部疾病的疾病相关感知网络Photoby阿辛·K·苏雷什on不飞溅介绍ChestX-Ray(CXR)成像是临床诊断中最常见的诊断成像......
  • pytorch相关知识点
    1.super的作用super()函数可以隐式的将子类里的method和父类里的method进行关联,这样就不需要再子类里显式的重新创建父类method里的属性说人话就是继承父类需要对用到的......