1.两数不能组成的最大的数
引理:给定a,ba,b,若d=gcd(a,b)>1d=gcd(a,b)>1,则一定不能凑出最大数
结论:如果 a,ba,b 均是正整数且互质,那么由 ax+by,x≥0,y≥0不能凑出的最大数是 (a−1)(b−1)−1
a,b不能构成的最大的数为a*b-a-b;
2,连续区间的判断
对于1到N的不确定顺序的数列,对【a,b】区间,若max【a,b】-min【a,b】==b-a,则该区间数递增排序后是连续序列
3.求约数
一个数可以写成若干p^c相乘
例:100=2的平方乘5的平方,然后100的约数可由2的次方和5的次方组合
#include<bits/stdc++.h>
#define ll long long
#define PI 3.141592653589793
#define E 2.718281828459045
#define HalF (l + r)>>1
#define lsn rt<<1
#define rsn rt<<1|1
#define Lson lsn, l, mid
#define Rson rsn, mid+1, r
#define QL Lson, ql, qr
#define QR Rson, ql, qr
#define myself rt, l, r
#define mem(a,b) memset(a,b,sizeof(a))
#define FOR( i , a , b ) for ( int i = a ; i <= b ; ++ i )
#define FO( i , n ) for ( int i = 0 ; i < n ; ++ i )
#define lowbit(a) ((a)&-(a))
#define PII pair<ll ,ll >
#define ft first
#define sd second
typedef unsigned long long ull;
const ll mod=10007;
const ll INF=0x3f3f3f3f;
const ll Max=1e4+10;
using namespace std;
ll t,n,m,l;
ll ans;
ll a[Max];
/*queue<ll> q;
stack<ll> s;
//升序队列
priority_queue <int,vector<int>,greater<int> > q;
//降序队列
priority_queue <int,vector<int>,less<int> >q;*/
int main()
{
ios::sync_with_stdio(false);
while(cin>>n,n)
{
m=1;
for(ll i=2;n!=1;i++)
{
ll w=1;
while(n%i==0)
{
w++;
n/=i;
}
m*=w;
}
cout<<m<<endl;
}
return 0;
}
4.最大公约数(GCD)
利用的是辗转相除法,具体的证明过程太麻烦了,有时间的可以去找找,下面是代码
int gcd(int x,int y)
{
if(y==0)
return x;
return gcd(y,x%y);
}
5.最大公倍数(LCM)
a,b的最大公倍数等于a,b相乘除以a,b的最大公约数
int LCM(int x,int y)
{
return (x*y)/gcd(x,y);
}
6.拓展欧几里得算法
解决ax+by=1的问题
定理一:
定理二:
算法:
.
7.计算斜率和截距
double a=(double)(y2-y1)/(double)(x1-x2);//斜率
double b=(double)(x2*y1-x1*y2)/(double)(x2-x1);//截距
8.循环小数化为分数
例子;
0.12121212
则可以化为:12/99,
a/b=(循环体)/(10^(循环体长度)-1),a为分子,b为分母