[Hdu] 数学基础 day2

Sun, 25 May 2025 00:00:00 +0000

快速幂

求 A^B 的最后三位数表示的整数

输入格式

输入数据包含多个测试实例
每个实例占一行，由两个正整数A和B组成（1 ≤ A,B ≤ 10000）
如果A=0且B=0，则表示输入结束，不做处理

输出格式

对于每个测试实例，输出A^B的最后三位表示的整数
每个输出占一行

输入输出示例

输入样例	输出样例
2 3	8
12 6	984
6789 10000	1
0 0	(无输出)

算法讲解 :

对于 123^(234). 我们可以优化成 (123^2)^117, 再进一步优化成 ((123^2) * 123 )^116
同理往下, 我们可以发现, 我们的计算量由次幂来控制, 同时我们可以根据就性质来操作次幂的优化 3. 对于奇数我们只需要拆成 ans = ans * a 的形式 4. 对于偶数我们只需要拆成 a = a * a 的形式

#include <bits/stdc++.h>
using namespace std;
int qmi(int a, int b) {
    int ans = 1 ; 
    while(b != 0) {
        if(b&1) ans = ans * a % 1000 ; 
        a = a  * a % 1000 ; 
        b = b >> 1; 
    }
    return ans % 1000 ; 
}
int main() { 
    int a, b; 
    while(1) {
        cin>>a>>b;
        if (a  == 0 && b == 0 ) {
            break ; 
        }
        cout<<qmi(a,b)<<endl;
    }
    return 0;
}

vj题目链接

todo

[Hdu] 数学基础

Sun, 04 May 2025 00:00:00 +0000

引导问题

整数求和

Hey, welcome to HDOJ(Hangzhou Dianzi University Online Judge).

In this problem, your task is to calculate SUM(n) = 1 + 2 + 3 + ... + n.

Input
The input will consist of a series of integers n, one integer per line.

Output
For each case, output SUM(n) in one line, followed by a blank line. You may assume the result will be in the range of 32-bit signed integer.

常规做法是 for-each 后对 sum 进行累加，然后输出
我们可以根据等差数列高斯公式求得 (a + b) * n / 2
不过这题需要注意的是 4. 我们计算乘法的时候可能暴 int , 因为对于两个 int32 的数相乘必然会爆 int32 所以我们需要开辟成 int64

题目链接如下, 由于不支持 GO 所以不做

https://vjudge.net/problem/HDU-1001

例题

最小公倍数

输入两个整数 a 和 b，请你编写一个函数，int lcm(int a, int b)，计算并输出 a 和 b 的最小公倍数。

输入格式
共一行，包含两个整数 a 和 b。

输出格式
共一行，包含一个整数，表示 a 和 b 的最小公倍数。

数据范围
1≤a,b≤1000

输入样例：
6 8
输出样例：
24

朴素做法

我们从最大的数开始枚举，如 8,9,10… 每次都判断是否能整除最小数
简单优化，我们可以直接枚举最大数的倍数。可是最大数的倍数, 我们有很大概率爆int 的风险

正解

lcm(A,B) = A * B / gcd(A,B)
优化一下 = A / gcd(A,B) * B
问题引导为 如何求最大公约数

Acwing-最小公倍数

最大公约数

我们现在要求 10, 14 的最大公约数。
我们假设 X 为这两个数的最大公约数，可以知道，对于 (14%10) 的余数，也应该是 X 的倍数
从而我们可以依次类推 (10,14) , (10,4) , (2,4), (2,0) 则 2 就是最大公约数

todo : 很抓马的一件事，我用 go 写这个程序 TLE 了

package main

import(
    "fmt"
)

func gcd(a,b int) int {
    if a > b {
        a,b = b,a
    }    
    for ; a != 0 ;  {
        a, b = b%a,a
    }
    return b 
}

func main() {
    var n int 
    fmt.Scan(&n)
    for i := 0 ; i < n ; i ++ {
        var a,b int 
        fmt.Scan(&a,&b)
        fmt.Println(gcd(a,b))
    }
}

#include <iostream> 
using namespace std ;
int gcd(int a,int b)  {
    int temp ; 
    if(a > b) {
        temp = a;
        a = b ;
        b = temp;
    } 
    for( ; a != 0 ; ) {
        temp = a;
        a = b%a;
        b = temp ;
    }
    return b ; 
}

int main() {
    int a,b ;
    int n;
    cin>>n;
    for(int i  = 1; i <= n ; i ++ ) {
        cin>>a>>b;
        cout<<gcd(a,b)<<endl;
    }
}

循环节 - N * N 的个位数

Given a positive integer N, you should output the most right digit of N^N.

Input

The input contains several test cases. The first line of the input is a single integer T which is the number of test cases.

T test cases follow.

Each test case contains a single positive integer N(1<=N<=1,000,000,000).

Output

For each test case, you should output the rightmost digit of N^N.

时间范围 : T * N 必然做法是 T * O(K) 的

2,4,8,16,32,64,128 ... 3,9,27,81,243 ...

显然，奇数 * 奇数 = 奇数，偶数 * 偶数 = 偶数。个位上的偶数也就 (0,2,4,6,8) , 奇数 (1,3,5,7,9) 所以肯定是一个循环,
因此我们只需要枚举出前5个,找出不重复的循环节即可解出这道题

#include <iostream> 
#include <math.h>
using namespace std;
int main() {
    int t;cin>>t;
    while(t -- ) {
        int n ; cin >> n; 
        int st[10] = {0};
        
        int last = n % 10 ;
        int ans[6] = {0} ;
        int idx = 0 ;
        
        for(int i = 0 ; i < 5 ; i++ ) {
            int k = int(pow(last, i + 1)) % 10; 
            if(st[k]){
                break ; 
            }
            ans[++idx] = k ; 
            st[k] = 1 ; 
        }
        ans[0] = ans[idx] ; 
        cout<<ans[(n%(idx)) ] <<endl; 
        idx = 0 ;
    }
    return 0; 
}

题目链接

Hdu-Icpc on GzzhM