数学专题1

最新推荐文章于 2023-11-09 13:57:37 发布
最新推荐文章于 2023-11-09 13:57:37 发布
阅读量318 收藏
点赞数
分类专栏： 数学 文章标签： 数论 矩阵快速幂 素数
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/leohujx/article/details/49025475
版权

A:相当于已知 f[i] = 2*(f[i-1]+f[i-2])，然后求f(n),n大到2^60，所以我们无法用递推，此时我们就用矩阵快速幂。

#include<stdio.h>
#include<iostream>
#include<string>
#include<string.h>
#include<algorithm>
#include<iomanip>
#include<vector>
#include<time.h>
#include<queue>
#include<stack>
#include<iterator>
#include<math.h>
#include<stdlib.h>
#include<limits.h>
#include<map>
#include<set>
#include<bitset>
//#define ONLINE_JUDGE
#define eps 1e-5
#define INF 0x7fffffff
#define FOR(i,a) for((i)=0;i<(a);(i)++)
#define MEM(a) (memset((a),0,sizeof(a)))
#define sfs(a) scanf("%s",a)
#define sf(a) scanf("%d",&a)
#define sfI(a) scanf("%I64d",&a)
#define pf(a) printf("%d\n",a)
#define pfI(a) printf("%I64d\n",a)
#define pfs(a) printf("%s\n",a)
#define sfd(a,b) scanf("%d%d",&a,&b)
#define sft(a,b,num) scanf("%d%d%d",&a,&b,&num)
#define for1(i,a,b) for(int i=(a);i<b;i++)
#define for2(i,a,b) for(int i=(a);i<=b;i++)
#define for3(i,a,b)for(int i=(b);i>=a;i--)
#define MEM1(a) memset(a,0,sizeof(a))
#define MEM2(a) memset(a,-1,sizeof(a))
#define ll long long
const double PI=acos(-1.0);
template<class T> T gcd(T a,T b){return b?gcd(b,a%b):a;}
template<class T> T lcm(T a,T b){return a/gcd(a,b)*b;}
template<class T> inline T Min(T a,T b){return a<b?a:b;}
template<class T> inline T Max(T a,T b){return a>b?a:b;}
using namespace std;
//#pragma comment(linker,"/STACK:1024000000,1024000000")
ll n,m,L;
#define N 210
#define M 1000100
#define Mod 1000000007
#define p(x,y) make_pair(x,y)
int a[100010];
struct Matrix{
    ll a[2][2];
    Matrix(){
        memset(a,0,sizeof a);
    }
};
Matrix Multi(Matrix c, Matrix d) {
    Matrix s;
    for (int i = 0; i < 2; i++)
        for (int j = 0; j < 2; j++)
            for (int k = 0; k < 2; k++) {
                    s.a[i][j] = (s.a[i][j] + c.a[i][k] * d.a[k][j]) % Mod;
            }
    return s;
}
Matrix Multi1(Matrix c, Matrix d) {
    Matrix s;
    for (int i = 0; i < 2; i++)
        for (int j = 0; j < 1; j++)
            for (int k = 0; k < 2; k++) {
                    s.a[i][j] = (s.a[i][j] + c.a[i][k] * d.a[k][j]) % Mod;
            }
    return s;
}
Matrix matrix;
void Init(){
    matrix.a[0][0]=2,matrix.a[0][1]=2;
    matrix.a[1][0]=1,matrix.a[1][1]=0;
}
ll quick_pow(ll n){
    Matrix unit;
    unit.a[0][0] = 8;
    unit.a[1][0] = 3;
    while(n){
        if(n&1)
            unit=Multi1(matrix,unit);
            matrix=Multi(matrix,matrix);
            n>>=1;
    }
    ll res=unit.a[0][0]*1%Mod;
    return res;
}
int main() {
    while(sfI(n)!=EOF){
        Init();
        if(n == 1)
            printf("3\n");
        else if(n == 2)
            printf("8\n");
        else{
            ll ans = quick_pow(n-2);
            printf("%I64d\n",ans);
        }
    }
    return 0;
}

B:求一个区间内的所有回文素数。区间右边界可达到1e8，所以我们不可能求出所有素数。我们看，既然是回文的，那么我们可以构造啊，构造出一个回文数，然后判断下是否为素数就可以了。我们还知道，除了11，是没有偶数位的回文素数了，因为偶数位的回文数都可以被11整除。我们只需要构造奇数位的回文数就行了。这里我们拿5位回文数来看：XYZYX，所以我们只需要3个不同的数，

for(int i=1;i<10;i+=2)
  for(int j=0;j<10;j++)
      for(int k=0;k<10;k++)
         tmp = i*10001+j*1010+k*100;

这样子我们就可以构造出一个5位的回文数，然后用素数的定义判断下它是否为素数即可。

#include<stdio.h>
#include<iostream>
#include<string>
#include<string.h>
#include<algorithm>
#include<iomanip>
#include<vector>
#include<time.h>
#include<queue>
#include<stack>
#include<iterator>
#include<math.h>
#include<stdlib.h>
#include<limits.h>
#include<map>
#include<set>
#include<bitset>
//#define ONLINE_JUDGE
#define eps 1e-5
#define INF 0x7fffffff
#define FOR(i,a) for((i)=0;i<(a);(i)++)
#define MEM(a) (memset((a),0,sizeof(a)))
#define sfs(a) scanf("%s",a)
#define sf(a) scanf("%d",&a)
#define sfI(a) scanf("%I64d",&a)
#define pf(a) printf("%d\n",a)
#define pfI(a) printf("%I64d\n",a)
#define pfs(a) printf("%s\n",a)
#define sfd(a,b) scanf("%d%d",&a,&b)
#define sft(a,b,num) scanf("%d%d%d",&a,&b,&num)
#define for1(i,a,b) for(int i=(a);i<b;i++)
#define for2(i,a,b) for(int i=(a);i<=b;i++)
#define for3(i,a,b)for(int i=(b);i>=a;i--)
#define MEM1(a) memset(a,0,sizeof(a))
#define MEM2(a) memset(a,-1,sizeof(a))
#define ll long long
const double PI=acos(-1.0);
template<class T> T gcd(T a,T b){return b?gcd(b,a%b):a;}
template<class T> T lcm(T a,T b){return a/gcd(a,b)*b;}
template<class T> inline T Min(T a,T b){return a<b?a:b;}
template<class T> inline T Max(T a,T b){return a>b?a:b;}
using namespace std;
//#pragma comment(linker,"/STACK:1024000000,1024000000")
int n,m,L;
#define N 210
#define M 1000100
#define Mod 1000000007
#define p(x,y) make_pair(x,y)
vector<int>v;
bool check(int x){
    for(int i=2;i*i<=x;i++){
        if(x%i == 0) return false;
    }
    return true;
}
void Init(){
    v.push_back(2);
    v.push_back(3);
    v.push_back(5);
    v.push_back(7);
    v.push_back(11);
    for(int i=1;i<10;i+=2){  //构造三位回文素数
        for(int j=0;j<10;j++){
            int tmp = i*101+10*j;
            if(check(tmp)) v.push_back(tmp);
        }
    }
    for(int i=1;i<10;i+=2){      //构造五位回文素数
        for(int j=0;j<10;j++){
            for(int k=0;k<10;k++){
                int tmp = i*10001+j*1010+100*k;
                if(check(tmp)) v.push_back(tmp);
            }
        }
    }

    for(int i=1;i<10;i+=2){  //构造七位回文素数
        for(int j=0;j<10;j++){
            for(int k=0;k<10;k++){
                for(int s=0;s<10;s++){
                    int tmp = i*1000001+j*100010+k*10100+s*1000;
                    if(check(tmp)) v.push_back(tmp);
                }
            }
        }
    }
}
int main() {
  Init();
//  printf("%d\n",(int)v.size());
    while(sfd(n,m)!=EOF){
        int size = (int)v.size();
        for(int i=0;i<size;i++){
            if(v[i]>=n&&v[i]<=m)
                printf("%d\n",v[i]);
        }
    }
    return 0;
}

C:给出n,k，求1-n中有多少个数的最大素因子小于等于k，n的范围是1<=N<=10^10,2<=k<=100,我们可以看出,k很小，所以我们可以从k下手。
我们利用<=k之内的所有素数，来通过dfs构造所有小于等于N且最最大素因子<=k的数。

#include<stdio.h>
#include<iostream>
#include<string>
#include<string.h>
#include<algorithm>
#include<iomanip>
#include<vector>
#include<time.h>
#include<queue>
#include<stack>
#include<iterator>
#include<math.h>
#include<stdlib.h>
#include<limits.h>
#include<map>
#include<set>
#include<bitset>
//#define ONLINE_JUDGE
#define eps 1e-5
#define INF 0x7fffffff
#define FOR(i,a) for((i)=0;i<(a);(i)++)
#define MEM(a) (memset((a),0,sizeof(a)))
#define sfs(a) scanf("%s",a)
#define sf(a) scanf("%d",&a)
#define sfI(a) scanf("%I64d",&a)
#define pf(a) printf("%d\n",a)
#define pfI(a) printf("%I64d\n",a)
#define pfs(a) printf("%s\n",a)
#define sfd(a,b) scanf("%d%d",&a,&b)
#define sft(a,b,num) scanf("%d%d%d",&a,&b,&num)
#define for1(i,a,b) for(int i=(a);i<b;i++)
#define for2(i,a,b) for(int i=(a);i<=b;i++)
#define for3(i,a,b)for(int i=(b);i>=a;i--)
#define MEM1(a) memset(a,0,sizeof(a))
#define MEM2(a) memset(a,-1,sizeof(a))
#define ll long long
const double PI=acos(-1.0);
template<class T> T gcd(T a,T b){return b?gcd(b,a%b):a;}
template<class T> T lcm(T a,T b){return a/gcd(a,b)*b;}
template<class T> inline T Min(T a,T b){return a<b?a:b;}
template<class T> inline T Max(T a,T b){return a>b?a:b;}
using namespace std;
//#pragma comment(linker,"/STACK:1024000000,1024000000")
ll n,m,L;
#define N 1000005
#define M 1000100
#define Mod 1000000007
#define p(x,y) make_pair(x,y)
int prime[30] = {2,3,5,7,11,13,17,19,23,29,31,37,41,43,47,53,59,61,67,71,73,79,83,89,97,101};
ll ans;
void dfs(ll x,int k){
    ans++;
    for(int i=k;prime[i]<=m && x*prime[i]<=n;i++)
        dfs(x*prime[i],i);
}
int main() {
    while(scanf("%I64d%I64d",&n,&m)!=EOF){
        ans=0;
        dfs(1ll,0);
        printf("%I64d\n",ans);
    }
    return 0;
}

D.求出一个N，使得N!%(x^y)=0，1<=X<=1,000,000,000, 1<=Y<=30,000。这种题目一般肯定是素数分解，因为任意一个数都可以被分解为素数相乘的形式。
我们先打个素数表，范围到sqrt(X)就行了，然后开始分解x,用两个数组，一个记录x能被哪几个素数整除，另外一个数组记录该素数出现了几次，次数其实就是单个x中出现的次数*y就行了。然后开始计算N，因为N越大，它所含的素数就越多，所以这是个单调的过程，我们可以用二分，得到一个N，我们就计算它N!中所含的X^Y中出现过的素数的个数，如果小于X^Y中出现的次数，那么说明N还不够大，左边界变大，否则右边界变小。

#include<stdio.h>
#include<iostream>
#include<string>
#include<string.h>
#include<algorithm>
#include<iomanip>
#include<vector>
#include<time.h>
#include<queue>
#include<stack>
#include<iterator>
#include<math.h>
#include<stdlib.h>
#include<limits.h>
#include<map>
#include<set>
#include<bitset>
//#define ONLINE_JUDGE
#define eps 1e-5
#define INF 0x7fffffff
#define FOR(i,a) for((i)=0;i<(a);(i)++)
#define MEM(a) (memset((a),0,sizeof(a)))
#define sfs(a) scanf("%s",a)
#define sf(a) scanf("%d",&a)
#define sfI(a) scanf("%I64d",&a)
#define pf(a) printf("%d\n",a)
#define pfI(a) printf("%I64d\n",a)
#define pfs(a) printf("%s\n",a)
#define sfd(a,b) scanf("%d%d",&a,&b)
#define sft(a,b,num) scanf("%d%d%d",&a,&b,&num)
#define for1(i,a,b) for(int i=(a);i<b;i++)
#define for2(i,a,b) for(int i=(a);i<=b;i++)
#define for3(i,a,b)for(int i=(b);i>=a;i--)
#define MEM1(a) memset(a,0,sizeof(a))
#define MEM2(a) memset(a,-1,sizeof(a))
#define ll long long
const double PI=acos(-1.0);
template<class T> T gcd(T a,T b){return b?gcd(b,a%b):a;}
template<class T> T lcm(T a,T b){return a/gcd(a,b)*b;}
template<class T> inline T Min(T a,T b){return a<b?a:b;}
template<class T> inline T Max(T a,T b){return a>b?a:b;}
using namespace std;
//#pragma comment(linker,"/STACK:1024000000,1024000000")
ll n,m,L;
#define N 50005
#define M 1000100
#define Mod 1000000007
#define p(x,y) make_pair(x,y)
ll f[N],prime[N];
ll fac[N],num[N];
int tot;
int len;
void Init(){            //素数打表(到sqrt(n)即可)
    memset(f,0,sizeof f);
    tot=0;
    for(int i=2;i<N;i++){
        if(!f[i])
            f[i] = prime[tot++] = i;
        for(int j=0;j<tot && prime[j]<=f[i] && i*prime[j]<N;j++)
            f[i*prime[j]] = prime[j];
    }
    return;
}
void div1(ll x,ll y){       //质因数分解
    len = 0;
    for(int i=0;i<tot && prime[i]*prime[i]<=x;i++){
        if(x%prime[i] == 0){
            ll s = 0;
            while(x%prime[i] == 0){
                x/=prime[i];
                s++;
            }
            fac[len] = prime[i];    //记录素数
            num[len++] = s*y;       //记录这个素数在x^y里面出现了几次
        }
    }
    if(x > 1){
        fac[len] = x;
        num[len++] = y;
    }
    return;
}
bool check(ll x){
    for(int i=0;i<len;i++){
        ll tmp = x,p=0;
        while(tmp){     //计算tmp!里面中几个fac[i]因子
            p += (tmp/fac[i]);
            tmp /= fac[i];
        }
        if(p<num[i]) return false;  //如果tmp!里面含有的fac[i]因子小于x^y里面fac[i]的因子数，说明x不够大
    }
    return true;
}
int main() {
     Init();
     ll X,Y;
    while(scanf("%I64d%I64d",&X,&Y)!=EOF){
        div1(X,Y);
        ll l=0,r=(1ll<<60);
        while(l<=r){
            ll mid = (l+r)>>1;
            if(check(mid))
                r = mid-1;
            else
                l = mid+1;
        }
        printf("%I64d\n",l);
    }
    return 0;
}
leohujx
关注
  • 0
    点赞
  • 踩
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
数学专题
WMXNLFD的博客
06-11 288
LeetCode 268. Missing Number 给定一个包含 0, 1, 2, …, n 中 n 个数的序列，找出 0 … n 中没有出现在序列中的那个数。 示例 1: 输入: [3,0,1] 输出: 2 示例 2: 输入: [9,6,4,2,3,5,7,0,1] 输出: 8 说明: 你的算法应具有线性时间复杂度。你能否仅使用额外常数空间来实现? class Solution { pub...
ZZULIOJ 1096: 水仙花数（函数专题） (C/C++)
无
05-10 836
在不知道m到n之间有多少个水仙花数的情况下，如何控制两个数之间有空格，而最后一个水仙花数之后没有空格？对于每个测试实例，要求输出所有在给定范围内的水仙花数，就是说，输出的水仙花数必须大于等于m,并且小于等于n，如果有多个，则要求从小到大排列在一行内输出，之间用一个空格隔开;要求程序定义一个narcissus()函数和一个main()函数，narcissus()函数判断一个整数n是否为水鲜花数，其余功能在main()函数中实现。输入数据有多组，每组占一行，包括两个整数m和n（100<=m<=n<=999）。
（小白看了都会！）zzulioj 1096: 水仙花数（函数专题）
Yolo_lxt的博客
06-14 1642
哥们儿，姐们儿，代买部分可复制粘贴，但是希望你能看看注释理解之后，自己敲出来，如果能给我一个小小的赞就更好啦~谢谢你~ 代码及注释部分： #include <stdio
ZZULIOJ-1096，水仙花数（函数专题）（Python）
manianuts的博客
10-04 303
【代码】ZZULIOJ-1096，水仙花数（函数专题）（Python）
python 1096: 水仙花数（函数专题）
m0_52141241的博客
12-20 1078
while True: m,n=map(int,input().split()) e=[] for i in range(m,n+1): b=i//100 c=(i%100)//10 d=i%10 if ib3+c3+d**3: e.append(str(i)) if e[]: print(‘no’) else: print(’ '.join(e))
Oj1096 水仙花数（函数专题,多实例）
2303_79080543的博客
11-09 259
int n, m;= EOF)int t = 0;++t;++m;if(t == 0)return 0;x = i/100;y = i%10;z = i/10;n = z%10;if(t == 0)elsereturn 1;elsereturn 0;
小升初数学专题1：数与代数（2）数的整除、因数、倍数、合数、质数、奇数、偶数.pdf
07-02
小升初数学专题1：数与代数（2）数的整除、因数、倍数、合数、质数、奇数、偶数.pdf
小升初数学专题1：数与代数（1）整数与小数.pdf
07-02
小升初数学专题1：数与代数（1）整数与小数.pdf
小升初数学专题1：数与代数（5）简易方程.pdf
07-02
小升初数学专题1：数与代数（5）简易方程.pdf
水仙花函数
03-03
水仙花函数 功能：判断输入 nValue 是否为水仙花数 输入： nValue为正整数 输出：无 返回：如果输入为水仙花数 返回1 否则返回0 unsigned int IsDaffodilNum unsigned int nValue
小升初数学专题1：数与代数（7）比与比例.pdf
07-02
小升初数学专题1：数与代数（7）比与比例.pdf
小升初数学专题1：数与代数（4）运算与规律.pdf
07-02
小升初数学专题1：数与代数（4）运算与规律.pdf
zzulioj 1096: 水仙花数（函数专题）
11-09 7676
题目链接：http://acm.zzuli.edu.cn/problem.php?id=1096 1096: 水仙花数（函数专题） 时间限制:1Sec内存限制:128 MB 提交:20358解决:7802 [状态] [讨论版] [提交] [命题人:admin] 题目描述 春天是鲜花的季节，水仙花就是其中最迷人的代表，数学上有个水仙花数，他是这样定义的： “水仙花数”...
1096: 水仙花数(函数专题)_C语言 | 输出水仙花数
weixin_29009971的博客
01-29 839
例54：C语言编程输出100-1000之间所有的“水仙花数”，所谓的“水仙花数”是指一个3位数，其各位数字立方和等于该数本身。解析：153是一个水仙花数，因为153=1^3+5^3+3^3，逐个位数判断即可，写代码的时候要考虑到每一位。源代码演示：#include<stdio.h>//头文件 int main()//主函数 { int i,j,k,narcissistic_Nu...
1239水仙花数
u013087645的专栏
12-14 2358
题目描述 春天是鲜花的季节，水仙花就是其中最迷人的代表，数学上有个水仙花数，是这样定义的：  “水仙花数”是指一个三位数，它的各位数字的立方和等于其本身，比如：153=13+53+33。 现在要求输出所有在m和n范围内的水仙花数。 输入 输入数据有多组，每组占一行，包括两个整数m和n（100 输出 对于每个测试实例，要求输出所有在给定范围内的水仙花数，就是说，输出
并查集的几道题目
Triose的专栏
08-11 2726
先从最简单的开始： 题意： 一共n个强盗，m条线索，输出有多少个独立团伙（同伙的同伙是同伙，老鼠的儿子会打洞） 样例输入： 10 9 1 2 3 4 5 2 4 6 2 6 8 7 9 7 1 6 2 4 样例输出： 3 代码： #include #include #include #include #include #include #include #inc
1096: 水仙花数（函数专题）
热门推荐
青梅煮茶的博客
09-21 1万+
1096: 水仙花数（函数专题） Time Limit: 1 Sec  Memory Limit: 128 MB Submit: 12194  Solved: 4322 SubmitStatusWeb Board Description 春天是鲜花的季节，水仙花就是其中最迷人的代表，数学上有个水仙花数，他是这样定义的：  “水仙花数”是指一个三位数，它的各位数字的立方和等于
html css js网页设计
最新发布
07-12
HTML、CSS和JavaScript是构建网页和网站的基本技术，它们共同工作来创建和设计用户界面。下面是关于这三种技术的详细介绍： ### HTML (HyperText Markup Language) - **定义**：HTML是构建网页内容的标准标记语言。 - **作用**：用于创建网页的结构和内容，如段落、链接、图片、表格等。 - **语法**：使用标签（如 `<p>`, `<div>`, `<a>`, `<img>` 等）来定义网页元素。 ### CSS (Cascading Style Sheets) - **定义**：CSS是一种样式表语言，用于描述HTML文档的呈现方式。 - **作用**：用于设置网页的布局、颜色、字体和其他视觉元素。 - **语法**：通过选择器（如 `p`, `.myclass`, `#myid` 等）应用样式规则。 ### JavaScript - **定义**：JavaScript是一种脚本语言，通常用于网页上实现交互功能。 - **作用**：允许网页与用户进行交互，如响应用户操作、动态更新内容、动画效果等。 - **语法**：Java
2020年高考数学 考点分析与突破性讲练 专题22 简单线性规划 理.doc
10-12
2020年高考数学考点分析与突破性讲练专题22简单线性规划，是数学考试中的一个重要内容，涉及了二元一次不等式组、二元一次不等式的几何意义以及简单的二元线性规划问题的解决。考生需要掌握确定二元一次不等式组表示...

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值