HDU5666 Segment (俄罗斯乘法(快速积) / Java大数取模 / 大数乘法取模转换)

最新推荐文章于 2023-10-26 23:21:00 发布

原创最新推荐文章于 2023-10-26 23:21:00 发布 · 435 阅读

0 ·

本内容遵循CC 4.0 BY-SA版权协议

数学问题同时被 3 个专栏收录

3 篇文章

订阅专栏

3.位运算

2 篇文章

订阅专栏

2.计算技巧

2 篇文章

订阅专栏

本博客详细介绍了如何解决HDU5666问题，主要涉及俄罗斯乘法算法，Java大数取模以及大数乘法取模的转换技巧。针对给定的直线在第一象限内形成的三角形整数点数量，通过优化算法在大数运算下进行高效取模计算，包括O(log2q)复杂度的俄罗斯乘法和O(1)复杂度的Java大数及作差求余数方法。

题目链接：http://acm.hdu.edu.cn/showproblem.php?pid=5666

题目大意：给定 $T$ 组数据，每组数据给定 $q$ 和 $P$ ，问直线 $x+y=q$ 在第一象限形成的三角形内部有多少个整数点，结果对 $P$ 取模;

数据范围： $1\leq{T}\leq10$ ， $1\leq{q}\leq10^{18}$ ， $1\leq{P}\leq10^{18}$

题目解答：俄罗斯乘法(复杂度： $O(\log_2q)$ ) / Java大数取模(复杂度： $O(1)$ ) / 大数乘法取模转换为作差求余数(复杂度： $O(1)$ )
原题实际上为求 $\sum\limits_{i = 1}^{q-2} i\bmod P= (q-2)*(q-1)/2\bmod P$ ，由于 $q$ 的范围很大，直接乘爆long long，采用俄罗斯乘法，或者用Java大数
另外一种思路，由于 $a*b\bmod p = a*b-\lfloor a*b/P\rfloor*P$ ， $\lfloor \rfloor$ 表示向下取整，可以把大数乘法取模转换为作差求余数

题目代码：
1.俄罗斯乘法

#include <cstdio>
#include <cstring>
#include <iostream>
#include <bitset>
using namespace std;

//俄罗斯乘法（快速积）
long long mul(long long a, long long b, long long p)
{
    long long ans = 0;
    for(; b; b >>= 1)
    {
        if(b & 1) ans = (ans + a) % p;
        a = a * 2 % p;
    }
    return ans;
}

int main()
{
    int t;
    cin >> t;
    while(t--)
    {
        long long q, p, ans;
        cin >> q >> p;
        if(q % 2)
            ans = mul((q-1)/2, q-2, p) % p;
        else
            ans = mul(q-1, (q-2)/2, p) % p;
        cout << ans << endl;
    }
}

2.Java大数

import java.util.*;
import java.math.*;
public class Main {
    public static void main(String[] args) {
        Scanner cin=new Scanner(System.in);
        int x;
        x=cin.nextInt();
        while(x!=0)
        {
            x--;
            BigInteger c=new BigInteger("2");
            BigInteger e=new BigInteger("1");
            BigInteger a=cin.nextBigInteger();
            BigInteger d=a.subtract(c);
            BigInteger f=a.subtract(e);
            BigInteger b=cin.nextBigInteger();
            BigInteger ans=f.multiply(d);
            ans=ans.divide(c);
            System.out.println(ans.remainder(b));
        }
    }
}

3.大数乘法取模转换

#include <cstdio>
#include <cstring>
#include <iostream>
#include <bitset>
using namespace std;

long long mul(long long a, long long b, long long p)
{
    a %= p, b %= p;
    long long c = (long double)a * b / p;
    long long ans = a * b - c * p;
    if(ans < 0) ans += p;
    else if (ans >= p) ans -= p;
    return ans;
}

int main()
{
    int t;
    cin >> t;
    while(t--)
    {
        long long q, p, ans;
        cin >> q >> p;
        if(q % 2)
            ans = mul((q-1)/2, q-2, p) % p;
        else
            ans = mul(q-1, (q-2)/2, p) % p;
        cout << ans << endl;
    }
}