求 \(n^k\) 的因子和, \(n \leq 2^{16}, k \leq 20\)

Solution

\[\prod_i \frac{p_i^{q_ik+1}-1}{p_i-1}
\]

#include <bits/stdc++.h>
using namespace std; const int maxlen = 105; class HP {
public:
int len, s[maxlen];
HP() { (*this) = 0; }
HP(int inte) { (*this) = inte; }
HP(const char *str) { (*this) = str; }
friend ostream &operator<<(ostream &cout, const HP &x);
HP operator=(int inte);
HP operator=(const char *str);
HP operator=(const HP &b);
HP operator*(const HP &b);
HP operator+(const HP &b);
HP operator-(const HP &b);
HP operator/(const HP &b);
HP operator%(const HP &b);
bool operator<(const HP &b);
bool operator>(const HP &b);
int Compare(const HP &b);
}; ostream &operator<<(ostream &cout, const HP &x) {
for (int i = x.len; i >= 1; i--) cout << x.s[i];
return cout;
}
HP HP::operator=(const char *str) {
len = (int)strlen(str);
for (int i = 1; i <= len; i++) s[i] = str[len - i] - '0';
return *this;
}
HP HP::operator=(int inte) {
if (inte == 0) {
len = 1;
s[1] = 0;
return (*this);
}
for (len = 0; inte > 0;) {
s[++len] = inte % 10;
inte /= 10;
}
return *this;
}
HP HP::operator=(const HP &b) {
len = b.len;
for (int i = 1; i <= len; i++) s[i] = b.s[i];
return *this;
}
HP HP::operator*(const HP &b) {
int i, j;
HP c;
c.len = len + b.len;
for (i = 1; i <= c.len; i++) c.s[i] = 0;
for (i = 1; i <= len; i++)
for (j = 1; j <= b.len; j++) c.s[i + j - 1] += s[i] * b.s[j];
for (i = 1; i < c.len; i++) {
c.s[i + 1] += c.s[i] / 10;
c.s[i] %= 10;
}
while (c.s[i]) {
c.s[i + 1] = c.s[i] / 10;
c.s[i] %= 10;
i++;
}
while (i > 1 && !c.s[i]) i--;
c.len = i;
return c;
}
HP HP::operator+(const HP &b) {
int i;
HP c;
c.s[1] = 0;
for (i = 1; i <= len || i <= b.len || c.s[i]; i++) {
if (i <= len)
c.s[i] += s[i];
if (i <= b.len)
c.s[i] += b.s[i];
c.s[i + 1] = c.s[i] / 10;
c.s[i] %= 10;
}
c.len = i - 1;
if (c.len == 0)
c.len = 1;
return c;
}
HP HP::operator-(const HP &b) {
int i, j;
HP c;
for (i = 1, j = 0; i <= len; i++) {
c.s[i] = s[i] - j;
if (i <= b.len)
c.s[i] -= b.s[i];
if (c.s[i] < 0) {
j = 1;
c.s[i] += 10;
} else
j = 0;
}
c.len = len;
while (c.len > 1 && !c.s[c.len]) c.len--;
return c;
}
int HP::Compare(const HP &y) {
if (len > y.len)
return 1;
if (len < y.len)
return -1;
int i = len;
while ((i > 1) && (s[i] == y.s[i])) i--;
return s[i] - y.s[i];
}
bool HP::operator<(const HP &b) {
if (len < b.len)
return 1;
if (len > b.len)
return 0;
int i = len;
while ((i > 1) && (s[i] == b.s[i])) i--;
return s[i] < b.s[i];
}
bool HP::operator>(const HP &b) {
if (len > b.len)
return 1;
if (len < b.len)
return 0;
int i = len;
while ((i > 1) && (s[i] == b.s[i])) i--;
return s[i] > b.s[i];
}
HP HP::operator/(const HP &b) {
int i, j;
HP d(0), c;
for (i = len; i > 0; i--) {
if (!(d.len == 1 && d.s[1] == 0)) {
for (j = d.len; j > 0; j--) d.s[j + 1] = d.s[j];
++d.len;
}
d.s[1] = s[i];
c.s[i] = 0;
while ((j = d.Compare(b)) >= 0) {
d = d - b;
c.s[i]++;
if (j == 0)
break;
}
}
c.len = len;
while ((c.len > 1) && (c.s[c.len] == 0)) c.len--;
return c;
}
HP HP::operator%(const HP &b) {
int i, j;
HP d(0);
for (i = len; i > 0; i--) {
if (!(d.len == 1 && d.s[1] == 0)) {
for (j = d.len; j > 0; j--) d.s[j + 1] = d.s[j];
++d.len;
}
d.s[1] = s[i];
while ((j = d.Compare(b)) >= 0) {
d = d - b;
if (j == 0)
break;
}
}
return d;
} int n,k,a[99],b[99],top;
HP ans;
int main() {
cin>>n>>k;
ans=1;
int t=n;
for(int i=2;i<=n;i++) {
if(t%i==0) {
a[++top]=i;
}
while(t%i==0) {
b[top]++;
t/=i;
}
}
for(int i=1;i<=top;i++) {
HP tmp=1,xx;
for(int j=1;j<=b[i]*k+1;j++) xx=a[i],tmp=tmp*xx;
xx=1;
tmp=tmp-xx;
xx=a[i]-1;
tmp=tmp/xx;
ans=ans*tmp;
}
cout<<ans;
}

[AHOI2002] 芝麻开门 - 数论的更多相关文章

  1. 2021record

    2021-10-14 P2577 [ZJOI2004]午餐 2021-10-13 CF815C Karen and Supermarket(小小紫题,可笑可笑) P6748 『MdOI R3』Fall ...

  2. Codeforces Round #382 Div. 2【数论】

    C. Tennis Championship(递推,斐波那契) 题意:n个人比赛,淘汰制,要求进行比赛双方的胜场数之差小于等于1.问冠军最多能打多少场比赛.题解:因为n太大,感觉是个构造.写写小数据, ...

  3. NOIP2014 uoj20解方程 数论(同余)

    又是数论题 Q&A Q:你TM做数论上瘾了吗 A:没办法我数论太差了,得多练(shui)啊 题意 题目描述 已知多项式方程: a0+a1x+a2x^2+..+anx^n=0 求这个方程在[1, ...

  4. 数论学习笔记之解线性方程 a*x + b*y = gcd(a,b)

    ~>>_<<~ 咳咳!!!今天写此笔记,以防他日老年痴呆后不会解方程了!!! Begin ! ~1~, 首先呢,就看到了一个 gcd(a,b),这是什么鬼玩意呢?什么鬼玩意并不 ...

  5. hdu 1299 Diophantus of Alexandria (数论)

    Diophantus of Alexandria Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 65536/32768 K (Java ...

  6. 【BZOJ-4522】密钥破解 数论 + 模拟 ( Pollard_Rho分解 + Exgcd求逆元 + 快速幂 + 快速乘)

    4522: [Cqoi2016]密钥破解 Time Limit: 10 Sec  Memory Limit: 512 MBSubmit: 290  Solved: 148[Submit][Status ...

  7. bzoj2219: 数论之神

    #include <iostream> #include <cstdio> #include <cstring> #include <cmath> #i ...

  8. hdu5072 Coprime (2014鞍山区域赛C题)(数论)

    http://acm.hdu.edu.cn/showproblem.php?pid=5072 题意:给出N个数,求有多少个三元组,满足三个数全部两两互质或全部两两不互质. 题解: http://dty ...

  9. ACM: POJ 1061 青蛙的约会 -数论专题-扩展欧几里德

    POJ 1061 青蛙的约会 Time Limit:1000MS     Memory Limit:10000KB     64bit IO Format:%lld & %llu  Descr ...

随机推荐

  1. Redis 为什么这么快?

    1. 纯内存操作,肯定快 数据存储在内存中,读取的时候不需要进行磁盘的 IO 2. 单线程,无锁竞争损耗 单线程保证了系统没有线程的上下文切换 使用单线程,可以避免不必要的上下文切换和竞争条件,没有多 ...

  2. Android布局管理器-使用TableLayout表格布局管理器实现简单的用户登录页面

    场景 Android布局管理器-使用FrameLayout帧布局管理器显示层叠的正方形以及前景照片: https://blog.csdn.net/BADAO_LIUMANG_QIZHI/article ...

  3. ACP知识总结

    由于ACP是一个敏捷开发的系统性知识,下面只针对我自身学习的知识总结,若需要完整的考试学习资料,可评论区或私聊我拿.   敏捷估计与规划.png   ACP知识点锦集.png   敏捷项目软件总结.p ...

  4. [PHP] 使用PHP在mongodb中进行嵌套查询

    作为文档数据库,数据库中存储的数据是类似json的结构,比如{“modelInfo”:{"status":1,audited:"1"}},想要查询status是 ...

  5. List集合去重各种方式汇总

    package com.sb.test; import java.util.*; import java.util.concurrent.ConcurrentHashMap; import java. ...

  6. Linux学习Day1:开班第一天

    其实这篇博客应该昨天就要写完的,算是补作业吧. 昨天(2020年2月14日)是参加Linux线上培训的第一天,当天培训结束后,老师要求学员每天写一篇博客来记录自己学到的知识,于是就有了这篇博客的诞生. ...

  7. PTA Is Topological Order

    Write a program to test if a give sequence Seq is a topological order of a given graph Graph. Format ...

  8. [TJOI2007] 足彩投注

    足彩投注 题目概述 题目背景 了解足球彩票的人可能知道,足球彩票中有一种游戏叫做"胜负彩",意为猜比赛的胜负.下面是一些与胜负彩有关的术语 注 :每一组有效组合数据. 投 注:彩民 ...

  9. QQ常用表情

    以下表情均为QQ官方表情原图,版权归QQ所有,禁止用于商业用途. ![3nEdY9.png](https://s2.ax1x.com/2020/02/21/3nEdY9.png) ![3nEaFJ.p ...

  10. linux 下查看json 文件 使用jq工具

    安装 文档 yum 安装 yum search jq yum -y install jq.x86_64 apt-get install jq jq支持查看 jq . json 文件 查看json文件 ...