POJ1845 首先把A写成唯一分解定理的形式 分解时让A对所有质数从小到大取模就好了 然后就有:A = p1^k1 * p2^k2 * p3^k3 *...* pn^kn 然后有: A^B = p1^(k1*B) * p2^(k2*B) *...* pn^(kn*B); 约数和公式: 对于已经分解的整数A=(p1^k1)*(p2^k2)*(p3^k3)*....*(pn^kn) 有A的所有因子之和为 S = (1+p1+p1^2+p1^3+...p1^k1) * (1+p2+p2^2+p2^3…
+; int prime[maxn]; void marktable(int n){ memset(prime,,sizeof(prime)); ;i<=n;i++){ ]]=i; ;j<=prime[]&&prime[j]<=n/i;j++){ prime[prime[j]*i]=; ) break; } } } ][]; int fatCnt; int getFactors(long long x){ fatCnt=; long long tmp=x; ;prime[…
import java.util.ArrayList; import java.util.List; //一个数如果恰好等于它的因子之和,这个数就称为"完数".例如6=1+2+3.第二个完全数是28, //它有约数1.2.4.7.14.28,除去它本身28外,其余5个数相加, //编程找出1000以内的所有完数. public class Test { public static void main(String[] args) { System.out.print("完数有…
day12 --------------------------------------------------------------- 实例019:完数 题目 一个数如果恰好等于它的因子之和,这个数就称为"完数".例如6=1+2+3.编程找出1000以内的所有完数. 分析:如果能被a整除,那么a是因子,放入一个列表然后再求和,与原被除数相等,就是完数 1 for i in range(1,int(a)): 2 list = [] 3 if i >1: 4 for j in r…
Sumdiv Time Limit:1000MS     Memory Limit:30000KB     64bit IO Format:%I64d & %I64u Submit Status Practice POJ 1845 Appoint description:   System Crawler  (2015-05-27) Description Consider two natural numbers A and B. Let S be the sum of all natural…
阶乘之和取模 (25 分) 输入正整数n, 计算S = 1!+2!+...+n!的末6位(不含前导0). 这里1<=n<=10​9​​. 输入样例: 例如输入: 20 输出样例: 输出: 820313 题解:我看道这个题的第一想法就是用暴力求解,但是最后一组数据时间超限了,然后我就打了一个表,发现当n>=24时,后面六位数都是一样的.所以就有了下面这个代码. #include<iostream> #include<cstdio> #include<algor…
新冠疫情深刻和全面地影响着社会和生活,已经成为数学建模竞赛的背景帝. 本文收集了与新冠疫情相关的的数学建模竞赛赛题,供大家参考,欢迎收藏关注. 『Python小白的数学建模课 @ Youcans』带你从数模小白成为国赛达人. 0. 前言:新冠疫情成了数模竞赛的背景帝 新冠疫情爆发以来,不仅严重影响到全球的政治和经济,也深刻和全面地影响着社会和生活的方方面面,甚至已经成为数学建模竞赛的背景帝. 传染病模型本来就是数学建模课程中的常见问题和模型.随着疫情的影响越来越严重.广泛和持久,不仅疫情传播.疫…
题目 poj的1845 分解a的质因数a=p1^t1*p2^t1........ 每个质因数对sum的贡献: 当除去质因数p1时的因数和为sum,当计入p1时,因子和变成sum*p1^0+sum*p1^1+sum*p1^2......+sum*p1^t1 也就是所有的sum=[1+p1+p1^2+p1^3+...+p1^t1]*[p2.....][p3...] 然后由于是a^b,所以最后是 sum=sum=[1+p1+p1^2+p1^3+...+p1^(t1*b)]*[p2.....][p3..…
ISSN Abbreviated Journal Title Full Title Category Subcategory Country total Cites IF        2013-2014 IF 2012-2013 IF 2011-2012 IF 2010-2011 IF 2009-2010 IF 2008-2009 IF 2007-2008 5-Year Impact Factor Immediacy Index Articles Cited Half-Life Eigenfa…
a.b只要数字a能被数字b整除,不论b是不是质数,都算是a的因子.比如:8的质因子是 2, 2, 2,但8的因子就包括 1,2,4. import math   for i in range(2, 1000):     factors = []  #因子列表,i 每次循环都清空     for j in range(1, math.floor(i/2)+1):         if i%j == 0:             factors.append(j)     if sum(factor…
package a; public class Wanshu { public static void main(String[] args) { for (int i = 1; i <= 1000; i++) { int sum=0; for (int j = 1; j < i; j++) { if(i%j==0) { sum+=j; } } if(i==sum) { System.out.println(i); } } } }…
理论依据: 代码: /* 显然,数据够大的时候,数组要用 __int64 */ #include<iostream> #include<map> #include<cstdio> #include<cstdlib> #include<cstring> using namespace std; ]; ];//用来%j ];//保存和. void make_ini() //全部扫一遍.筛选一下. { int i,j,k,t; ;i<=;i++)…
大整数的因子 总时间限制:  1000ms 内存限制:  65536kB 描述 已知正整数k满足2<=k<=9,现给出长度最大为30位的十进制非负整数c,求所有能整除c的k. 输入 一个非负整数c,c的位数<=30. 输出 若存在满足 c%k == 0 的k,从小到大输出所有这样的k,相邻两个数之间用单个空格隔开:若没有这样的k,则输出"none". 样例输入 30 样例输出 2 3 5 6 代码 高精模运算 #include<iostream> #inc…
题意:给出一大数K(4 <= K <= 10^100)与一整数L(2 <= L <= 106),K为两个素数的乘积(The cryptographic keys are created from the product of two primes) 问构成K的最小素数是否绝对小于L,若是,则输出BAD p,p为最小素数,否则输出GOOD; 分析:从小到大枚举1~10^6内的素数p,while(p<L)时,判断K是否能被p整除,若能则证明构成K的最小素数绝对小于L,反之则大于L…
n对mod求模,它的值在0到mod-1之间,如果要求模整除的时候转化成mod可以用下面的式子: n = (n - 1 % mod + mod) % mod +1 这里先减一,模上mod再加一,这样如果是整除mod的话先减一模上之后就变成mod-1,最后+1就变成mod了 这里模mod的方法是先模mod再加mod再模mod,这是常用的做法,目的是为了使得模的结过为正数,而不是负数.…
先放知识点: 莫比乌斯反演 卢卡斯定理求组合数 乘法逆元 快速幂取模 GCD of Sequence Alice is playing a game with Bob. Alice shows N integers a 1, a 2, -, a N, and M, K. She says each integers 1 ≤ a i ≤ M. And now Alice wants to ask for each d = 1 to M, how many different sequences b…
The Children's Day has passed for some days .Has you remembered something happened at your childhood? I remembered I often played a game called hide handkerchief with my friends. Now I introduce the game to you. Suppose there are N people played the…
一.知识目录 字符串处理 ................................................................. 3 1.KMP 算法 ............................................................ 3 2.扩展 KMP ............................................................ 6 3.Manacher 最长回文子串 .......…
\(>Topcoder \space Srm \space 673 \space Div2 \space 1000 \space BearPermutations2<\) 题目大意 : 对于一个长度为 \(n\) 的排列,定义其的贡献为对其建笛卡尔树,树上有两个儿子的节点其左右儿子在原排列中的距离之和,给出 \(n, Mod\),求所有长度为 \(n\) 的排列的贡献之和对 \(Mod\) 取模的值 \(1 \leq n \leq 100\) 解题思路 : 考虑一个最暴力的 \(dp\) ,设…
BUPT2017 wintertraining(16) #5 C HDU - 1021 题意 1, 12, 123, 1234, ..., 12345678910, ... 问第a到第b个数(inclusive)里有几个可以被3整除. 题解 前几项的各位数之和对3取模,可以找到规律--余数是1,0,0,1,0,0,1,0,0,1....也就是每三个数有两个可被3整除. 代码 #include <cstdio> #include <cstring> #include <algo…
目录 语法 c++ java 动态规划 多重背包 最长不下降子序列 计算几何 向量(结构体) 平面集合基本操作 二维凸包 旋转卡壳 最大空矩形 | 扫描法 平面最近点对 | 分治 最小圆覆盖 | 随机增量法 三维向量(结构体) 三维凸包 几何杂项 数据结构 ST表 单调队列 树状数组 线段树 并查集 左偏树 珂朵莉树,老司机树 莫队 二叉搜索树 一些建议 图论 图论的一些概念 图论基础 最短路径 最小生成树 树论的一些概念 最近公共祖先 联通性相关 图上的NP-hard问题 弦图+区间图 | 最…
题目链接 题意:求 A^B的所有约数之和对9901取模后的结果. 分析: 看了小优的博客写的. 分析来自 http://blog.csdn.net/lyy289065406/article/details/6648539 (1)   整数的唯一分解定理: 任意正整数都有且只有一种方式写出其素因子的乘积表达式. A=(p1^k1)*(p2^k2)*(p3^k3)*....*(pn^kn)   其中pi均为素数 (2)   约数和公式: 对于已经分解的整数A=(p1^k1)*(p2^k2)*(p3^…
2691. Sumdiv ★★★   输入文件:sumdiv.in   输出文件:sumdiv.out   简单对比时间限制:1 s   内存限制:12 MB [题目描述] 考虑两个自然数A和B.定义S是A ^ B的所有自然因数的总和.确定S模9901的值. [输入格式] 唯一的行包含由空格分隔的两个自然数A和B(0 <= A,B <= 50000000). [输出格式] 输出一行,即S模9901. [样例输入] 2 3 [样例输出] 15 样例解释:2 ^ 3 = 8. 8的自然因数是:1,…
MySQL基础知识 ----------- MySQL数据库简介 内存:一断电数据就没有了---------数据库:可以实现数据的持久化存储,有完整的管理系统管理,方便查询--------- DB:数据库,保存一组有组织的数据的容器DBMS:数据库管理系统,也叫数据库软件或产品,用于管理DB中的数据SQL:结构化查询语言,用于和DBMS通信的语言 DB:database 数据库:存储数据的大仓库,保存的是有组织,有规范的数据,DBMS:DatabseManagement System 数据库管理…
DIV2 有部分div1的题会写 div1的大佬真的太强了 向他们学习 (好像和zqc大佬说过话了hhh,zqc大佬真的是一个超有意思的人啊,羡慕有妹子队友的zqc大佬) A: 你有一棵树,你想把它画在平面上,使得没有边相交. 题解:dfs(u,dep)表示我在第dep层,第cnt[dep]列来放置第u个节点 代码: #include<bits/stdc++.h> using namespace std; ; int n,m; struct EDGE{ int v,nxt; }edge[max…
题目背景 xht喜欢研究数学函数,他特别喜欢反比例函数. 题目描述 我们知道,反比例函数xy=a的图象是双曲线. xht于是想:把它推广到三维是什么样的呢? 定义曲面C(k)为方程xyz=k所确定的曲面. 又定义曲面的美观程度P(k)为曲面C(k)上所有整点(x,y,z坐标均为整数)到原点的曼哈顿距离的平方之和. (点(x,y,z)到原点的曼哈顿距离为|x|+|y|+|z|). 现在,xht把一些曲面{C(a),C(a+1)...C(b)}排成一列,你要求出它们美观程度之和对10007取模的结果…
龟速补题.目前基本弃坑.已暂时放弃 D.I 两题. 下面不再写题意了直接说解法注意事项之类,直接放contest链接. https://vjudge.net/contest/151537 A.The Perfect Stall 很明显的二分图匹配,可以跑最大流 我是直接写的匈牙利算法,hungary,不是hungry (虽然读音差不多??? 事实证明浪了一个寒假我连hungary算法都不会写了 #include <cstdio> #include <vector> #include…
题目大意 从 \((0,0)\) 开始,每次只可走上下左右一个单位长度,可走重复路,求第 \(T\) 步正好走到 \((n,m)\) 的方案数. 答案要求对 \(MOD\) 取模,\(MOD\) 保证是几个不同质数的乘积. 数据范围 \(1\leq T\leq 100000;-T\leq n,m\leq T;1\leq MOD\leq 10^9+7\) 样例输入 4 10 2 2 样例输出 6 思路 枚举向左走了 \(l\) 步,则向右走了 \(r\) 步,向上走了 \(u\) 步,向下走了 \…
计算几何 A Farmer Greedy 题意:n个点选3个组成三角形,问m个点在三角形内的数字是奇数的这样的三角形个数. 分析:暴力O(N^3*M)竟然能过!我写的搓,加了优化才过掉.正解是先处理出每条线段正下方点的个数,然后枚举每个三角形O(1)计算,cnt[i][j] + cnt[j][k] - cnt[i][k](i,j,k已经对应的点按照x坐标排序),复杂度(n^2*m+n^3). 贴上别人的代码: #include <cstdio> #include <cstring>…
转自:http://blog.csdn.net/hellolwl/article/details/6164449 目录(?)[-] 协议介绍 模块协议介绍 1            命令包格式 2            数据包格式 下串口通信系统的组成 协议的实现 1            协议实现时的考虑 2            的实现方案 3            复用协议特殊情况的处理 摘要: 介绍 3GPP提出的一种终端设备和移动台串行通信的复用协议及嵌入式 Linux系统下串行通信模…