2020牛客国庆集训派对day2 F-SUM OF SUB RECTANGLE AREAS

本文介绍了一种计算所有子矩阵权值之和的方法,适用于N*N的矩阵,其中每个元素都为1。通过数学推导得出了简洁的公式,并使用Java大数BigInteger实现了高效算法。

2020牛客国庆集训派对day2 F-SUM OF SUB RECTANGLE AREAS 数学推导+大数相乘


传送门: https://ac.nowcoder.com/acm/contest/7818/F

题意

给一个N*N的矩阵,每一位上都是1,求所有子矩阵的权值之和。

思路

设 n = 3 ; 设n = 3; n=3
考虑子矩阵大小为 i ∗ j i * j ij的个数 x x x,即该大小的所有子矩阵的权值为 x ∗ i ∗ j x* i * j xij
1 ∗ 2 1 * 2 12子矩阵:每行个数为 2 2 2个,有 3 3 3列,即总权值为 2 ∗ 3 ∗ ( 1 ∗ 2 ) − − 2 ∗ 3 2 * 3 * (1 * 2)--2 * 3 231223为该矩阵在n*n矩阵中的个数, 1 ∗ 2 1 * 2 12为该子矩阵里的权值。
2 ∗ 3 2 * 3 23子矩阵:每两行有 1 1 1个,一共 2 2 2个两行,权值为 2 ∗ 3 2*3 23,即总权值为 1 ∗ 2 ∗ 2 ∗ 3 1 * 2 * 2 * 3 1223

根据上面推导出:在 n ∗ n n*n nn的矩阵中,有 i ∗ j i*j ij子矩阵的个数为 ( n − i + 1 ) ∗ ( n − j + 1 ) (n - i + 1) * (n - j + 1) (ni+1)(nj+1),即总权值为 ( n − i + 1 ) ∗ ( n − j + 1 ) ∗ i ∗ j (n - i + 1) * (n - j + 1) * i * j (ni+1)(nj+1)ij
即总答案为:
a n s = ∑ i = 1 n ∑ j = 1 n ( n − i + 1 ) ∗ ( n − j + 1 ) ∗ i ∗ j ans=\sum_{i=1}^n\sum_{j=1}^n(n - i + 1) * (n - j + 1) * i * j ans=i=1nj=1n(ni+1)(nj+1)ij

a n s = [ ∑ i = 1 n ( n − i + 1 ) i ] 2 ans=[\sum_{i=1}^n(n-i+1)i]^2 ans=[i=1n(ni+1)i]2

a n s = [ n 2 ( n + 1 ) 2 − n ( n + 1 ) ( 2 n + 1 ) 6 + n ( n + 1 ) 2 ] 2 ans=[\frac{n^2(n+1)}{2}-\frac{n(n+1)(2n+1)}{6}+\frac{n(n+1)}{2}]^2 ans=[2n2(n+1)6n(n+1)(2n+1)+2n(n+1)]2

a n s = [ n ∗ ( n + 1 ) ∗ ( n + 2 ) 6 ] 2 ans=[\frac{n * (n + 1)*(n+ 2)}{6}]^2 ans=[6n(n+1)(n+2)]2

因为题目说会爆longlong,所以我这里用的是JAVA大数BigInteger,python也可(不过我不会)。

Code(71MS)


import java.util.*;
import java.math.*;

public class Main
{
	public static void main(String[] args){
		Scanner in = new Scanner(System.in);
        int T;
        T = in.nextInt();
        for(int i = 1;i <= T; i++) {
		    long n;
		    n = in.nextLong();
		    BigInteger ans = new BigInteger("1");
            ans = ans.multiply(BigInteger.valueOf(n));
            ans = ans.multiply(BigInteger.valueOf(n + 1));
            ans = ans.multiply(BigInteger.valueOf(n + 2));
		    ans = ans.divide(BigInteger.valueOf(6));
		    ans = ans.multiply(ans);
		    System.out.println(ans);
	    }
    }
}
### 关于2020寒假算法基础集训营中的欧几里得算法 在2020年的寒假算法基础集训营中,确实存在涉及欧几里得算法的相关题目。具体来说,在第四场竞赛的第一题即为“A. 欧几里得”,该题目的核心在于利用扩展欧几里得定理来解决问题[^5]。 #### 扩展欧几里得算法简介 扩展欧几里得算法主要用于求解形如 ax + by = gcd(a, b) 的线性不定方程的一组特解(x,y),其中gcd表示最大公约数。此方法不仅能够计算两个整数的最大公因数,还能找到满足上述条件的具体系数x和y。 对于给定的数据范围较小的情况可以直接通过递归来实现;而对于较大数据则需考虑效率优化问题。下面给出了一段基于C++语言编写的用于解决此类问题的模板代码: ```cpp #include<bits/stdc++.h> #define int long long using namespace std; // 定义全局变量存储结果 int x, y; void ex_gcd(int a, int b){ if(b == 0){ x = 1; y = 0; return ; } ex_gcd(b, a % b); int tmp = x; x = y; y = tmp - (a / b) * y; } ``` 这段程序实现了经典的扩展欧几里得算法逻辑,并且可以作为处理类似问题的基础工具函数调用。 #### 实际应用案例分析 回到原题本身,“A. 欧几里得”的解答思路就是先预处理斐波那契数列前若干项数值存入数组`a[]`内以便快速查询,之后针对每一次询问直接输出对应位置处两相邻元素之和即可得出最终答案。这实际上巧妙运用到了广为人知的裴蜀定理——任意一对互质正整数都可由它们自身的倍数组合而成,而这里正是借助了这一性质简化了解决方案的设计过程。
评论 2
成就一亿技术人!
拼手气红包6.0元
还能输入1000个字符
 
红包 添加红包
表情包 插入表情
 条评论被折叠 查看
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值