【题目概括】

计算\(C_n^k\)的因子个数。

【思路要点】

  • 首先考虑将组合数展开,展开后就是\(\frac {n!}{k!\times (n-k)!}\)。
  • 这样就是计算出这些质因子的个数,然后将插值相乘就是答案了。
  • 但是数据较大,所以我们考虑一开始先将所有答案预处理来。
  • 如果是暴力求解,那么就会超时。
  • 我们考虑用归纳法来计算,如果计算出\(C_{n}^{k}\),那么如何计算出\(C_n^{k+1}\)。
  • 那么将式子展开,就相差了\(k+1\)和\(n-k+1\),就\(O(n)\)计算一次。
  • 时间复杂度 \(\mathcal O(n^3)\)
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <cmath>

using namespace std;

typedef long long ll;

const int N = 505;

int x, y, primeCnt;
bool ntPrime[N];
ll tot[N];
int prime[N];
ll ans[N][N];

void getPrime(int lim) {
  for (int i = 2; i <= lim; i++) {
    if (!ntPrime[i])
      prime[++primeCnt] = i;
    for (int j = 1; j <= primeCnt && i * prime[j] <= lim; j++) {
      ntPrime[i * prime[j]] = 1;
      if (i % prime[j])
        break;
    }
  }
}

void divide(int x, int tp) {
  for (int i = 1, cnt; i <= primeCnt; i++) {
    if (x % prime[i] == 0) {
      cnt = 0;
      while (x % prime[i] == 0) {
        x /= prime[i];
        cnt++;
      }
      tot[i] += tp * cnt;
    }
  }
}

ll solve(int x, int y) {
  for (int i = x; i >= (x - y + 1); i--)
    divide(i, 1);
  for (int i = 1; i <= y; i++)
    divide(i, -1);
  ll ans = 1;
  for (int i = 1; i <= primeCnt; i++)
    ans *= tot[i] + 1;
  return ans;
}

void init() {
  for (int i = 0; i <= 432; i++)
    for (int j = 0; j <= 432; j++)
      ans[i][j] = 1ll;
  for (int i = 1; i <= 432; i++) {
    memset(tot, 0, sizeof tot);
    ans[i][1] = solve(i, 1);
    for (int j = 2; j <= (i >> 1); j++) {
      divide(j, -1);
      divide(i - j + 1, 1);
      for (int k = 1; k <= primeCnt; k++)
        ans[i][j] *= tot[k] + 1;
    }
  }
}

int main() {
  getPrime(450);
  init();
  while (scanf("%d %d", &x, &y) != EOF) {
    y = min(y, x - y);
    printf("%lld\n", ans[x][y]);
  }
  return 0;
}

【POJ2992】Divisors的更多相关文章

  1. 【BZOJ】【4146】 【AMPPZ2014】Divisors

    暴力 由于值的范围很小($ \leq 2*10^6$),所以用一个cnt数组统计每个值有多少个数,然后从小到大,统计每个数的倍数即可. 根据调和数?的神奇性质= =这样是$O(nlogn)$的…… / ...

  2. 【BZOJ4146】[AMPPZ2014]Divisors

    [BZOJ4146][AMPPZ2014]Divisors Description 给定一个序列a[1],a[2],...,a[n].求满足i!=j且a[i]|a[j]的二元组(i,j)的个数. In ...

  3. Sigma Function (LightOJ - 1336)【简单数论】【算术基本定理】【思维】

    Sigma Function (LightOJ - 1336)[简单数论][算术基本定理][思维] 标签: 入门讲座题解 数论 题目描述 Sigma function is an interestin ...

  4. 【LeetCode】数学(共106题)

    [2]Add Two Numbers (2018年12月23日,review) 链表的高精度加法. 题解:链表专题:https://www.cnblogs.com/zhangwanying/p/979 ...

  5. Python高手之路【六】python基础之字符串格式化

    Python的字符串格式化有两种方式: 百分号方式.format方式 百分号的方式相对来说比较老,而format方式则是比较先进的方式,企图替换古老的方式,目前两者并存.[PEP-3101] This ...

  6. 【原】谈谈对Objective-C中代理模式的误解

    [原]谈谈对Objective-C中代理模式的误解 本文转载请注明出处 —— polobymulberry-博客园 1. 前言 这篇文章主要是对代理模式和委托模式进行了对比,个人认为Objective ...

  7. 【原】FMDB源码阅读(三)

    [原]FMDB源码阅读(三) 本文转载请注明出处 —— polobymulberry-博客园 1. 前言 FMDB比较优秀的地方就在于对多线程的处理.所以这一篇主要是研究FMDB的多线程处理的实现.而 ...

  8. 【原】Android热更新开源项目Tinker源码解析系列之一:Dex热更新

    [原]Android热更新开源项目Tinker源码解析系列之一:Dex热更新 Tinker是微信的第一个开源项目,主要用于安卓应用bug的热修复和功能的迭代. Tinker github地址:http ...

  9. 【调侃】IOC前世今生

    前些天,参与了公司内部小组的一次技术交流,主要是针对<IOC与AOP>,本着学而时习之的态度及积极分享的精神,我就结合一个小故事来初浅地剖析一下我眼中的“IOC前世今生”,以方便初学者能更 ...

随机推荐

  1. v-if 和v-show的区别

    在切换 v-if 块时,Vue.js 有一个局部编译/卸载过程,因为 v-if 之中的模板也可能包括数据绑定或子组件.v-if 是真实的条件渲染,因为它会确保条件块在切换当中合适地销毁与重建条件块内的 ...

  2. C语言I博客作业12-学期总结

    一.我学到的内容 二.我的收获 1. https://www.cnblogs.com/1076022899-lj/p/11576442.html 收获:第一次接触到C语言和PTA,第一次学习了博客园和 ...

  3. [Python3] 032 常用模块 random

    目录 random 1. random.random() 2. random.choice() 3. random.shuffle() 4. random.randint() 5. random.ra ...

  4. yarn以及mapreduce部署

    修改hadoop的配置文件yarn-site.xml: 复制该配置文件到其他服务器 scp yarn-site.xml ubuntu-01:$PWD yarn启动命令: start-yarn.sh M ...

  5. tree_cuttting(树形dp求解树的重心)

    Tree Cutting After Farmer John realized that Bessie had installed a "tree-shaped" network ...

  6. 开启linux服务器防火墙

    启用ufwsudo ufw enablesudo ufw default deny运行以上两条命令后,开启了防火墙,并在系统启动时自动开启.关闭所有外部对本机的访问,但本机访问外部正常. 开启和禁用s ...

  7. C#模态对话框和非模态对话框

    模态对话框弹出窗口阻止调用窗口的所有消息响应.只有在弹出窗口结束后调用窗口才能继续.在模态窗口“关闭”后,可以读取模态窗口中信息,包括窗口的返回状态,窗口子控件的值. 非模态对话框可以在弹出窗口和调用 ...

  8. LOJ 103 字串查找 题解

    题面 这道题是KMP的模板. KMP需要注意的细节有很多,所以把这篇文章发上来供参考: #include <bits/stdc++.h> using namespace std; char ...

  9. [Vue] vue的一些面试题4

    1.你知道 nextTick 的原理吗? 用法:在下次 DOM 更新循环结束之后执行延迟回调.在修改数据之后立即使用这个方法,获取更新后的 DOM. 异步更新队列提到 DOM 的更新是异步执行的,只要 ...

  10. QtSpim使用Tips

    QtSpim使用记录 垃圾QtSpim,输入中文会死机 MIPS的中文资料奇缺,如果有问题建议google参考英文资料,许多外国大学的网站上有对MIPS各方面的详细介绍 QtSpim是MIPS处理器的 ...