给一个n*m的格点图,问其中有多少个矩形?

$ \sum_{x=1}^{nm} \sum_{ab=x} [a + b \leq n](n - a - b + 1)\sum_{cd=x} [c + d \leq m](m - c - bd + 1)$

#include <bits/stdc++.h>
using namespace std;
typedef long long LL;
typedef long long ll;
#define rep(i, j, k) for (int i = int(j); i <= int(k); ++ i)
const int N = ;
LL f[N * N], g[N * N]; int main() {
LL n, m;
cin >> n >> m;
n --; m --;
rep(i, , n) rep(j, , n - i)
if (i * j < N * N) f[i * j] += (LL)(n - i - j + );
rep(i, , m) rep(j, , m - i)
if (i * j < N * N) g[i * j] += (LL)(m - i - j + );
LL ans = ;
rep(i, , N * N - ) {
// if (f[i] * g[i]) cout << i << '\n';
ans += f[i] * g[i];
}
ans += (n * (n + ) / ) * (m * (m + ) / );
cout << ans << '\n';
}

最新文章

  1. 查看mysql语句运行时间
  2. ExpandableListView 里面嵌套GridView实现高度自适应
  3. 从oracle数据表中读取表结构
  4. Class类的理解
  5. Spring的注入问题
  6. bzoj1019 [SHOI2008]汉诺塔
  7. C#调用存储过程详解
  8. nodejs hello world (调用文件执行没反应)
  9. angularjs uigrid 中celltemplate的写浮动框
  10. java关于jdbc的配置与使用步骤
  11. Java动态代理实现及实际应用
  12. 重新使用linux的一些事情
  13. 区间dp(入门题)
  14. shell读取mysql数据库
  15. Jquery 图片延迟加载技术
  16. HTML for android 移动小球
  17. Exception的ToString()方法究竟返回的是什么
  18. Kafka:ZK+Kafka+Spark Streaming集群环境搭建(二十)ES6.2.2 Client API
  19. Python - WebDriver 识别登录验证码
  20. Java IO流-File类

热门文章

  1. java JFR
  2. Ubuntu虚拟机识别不了SD卡的解决办法
  3. 20190404 Oracle忘记登陆密码
  4. [js]浏览器同源策略(same-origin policy)
  5. Delphi中类的运行期TypeInfo信息结构说明(转载)
  6. 158A Next Round
  7. 基于ROS和beaglebone的串口通信方式,使用键盘控制移动机器人
  8. php 数组数字 补零
  9. CSS 边框样式
  10. Python socket粘包解决