D2. Reverse Card (Hard Version)

D2. Reverse Card (Hard Version)

The two versions are different problems. You may want to read both versions. You can make hacks only if both versions are solved.

You are given two positive integers n, m.

Calculate the number of ordered pairs (a,b) satisfying the following conditions:

  • 1an, 1bm;
  • bgcd(a,b) is a multiple of a+b.

Input

Each test contains multiple test cases. The first line contains the number of test cases t (1t104). The description of the test cases follows.

The first line of each test case contains two integers n, m (1n,m2106).

It is guaranteed that neither the sum of n nor the sum of m over all test cases exceeds 2106.

Output

For each test case, print a single integer: the number of valid pairs.

Example

input

6
1 1
2 3
3 5
10 8
100 1233
1000000 1145141

output

0
1
1
6
423
5933961

Note

In the first test case, no pair satisfies the conditions.

In the fourth test case, (2,2),(3,6),(4,4),(6,3),(6,6),(8,8) satisfy the conditions.

 

解题思路

  令 d=gcd(a,b),则有 a=pdb=qd,且 gcd(p,q)=1。对于 (a+b)bd,等价于 (p+q)qd。又因为 gcd(q,p+q)=gcd(q,p)=1,因此有 (p+q)d

  因此 (a+b)bgcd(a,b) 等价于 (p+q)qd

  此时考虑能否枚举数对 (p,q)。事实上 p<d=npp2<n,同理 q2<m,因此这样的数对最多有 O(nm)。所以我们直接暴力枚举找到满足 gcd(p,q)=1pq,由于此时 d 最大能取到 min{np,mq},因此满足 (p+q)qd 对应的 (a,b) 的数量就是 min{np,mq}p+q

  AC 代码如下,时间复杂度为 O(nmlogmax{n,m})

#include <bits/stdc++.h>
using namespace std;

typedef long long LL;

void solve() {
    int n, m;
    scanf("%d %d", &n, &m);
    LL ret = 0;
    for (int i = 1; i * i < n; i++) {
        for (int j = 1; j * j < m; j++) {
            if (__gcd(i, j) == 1) ret += min(n / i, m / j) / (i + j);
        }
    }
    printf("%lld\n", ret);
}

int main() {
    int t;
    scanf("%d", &t);
    while (t--) {
        solve();
    }
    
    return 0;
}

 

参考资料

  Codeforces Round 942 (Div. 1, Div. 2) Editorial:https://codeforces.com/blog/entry/129027

posted @   onlyblues  阅读(127)  评论(0编辑  收藏  举报
相关博文:
阅读排行:
· 单线程的Redis速度为什么快?
· 展开说说关于C#中ORM框架的用法!
· Pantheons:用 TypeScript 打造主流大模型对话的一站式集成库
· SQL Server 2025 AI相关能力初探
· 为什么 退出登录 或 修改密码 无法使 token 失效
Web Analytics
点击右上角即可分享
微信分享提示