-
A -- Eular质数筛法
- Description
https://vjudge.net/problem/HihoCoder-1295
- Input
第1行:1个正整数n,表示数字的个数,2≤n≤1,000,000
- Output
第1行:1个整数,表示从1到n中质数的个数
- Sample Input
9
- Sample Output
4
- 题目理解
核心思想是一个非质数进行因数分解一定能够分解成为许多的质数因子,打一个数组标记质数和非质数,然后在外面进行一次
循环的遍历,当访问到一个质数的时候就将它的倍数都标记为非质数,通过这个方法就可以将所有的非质数都找出来剩下的就是质数了,然后求一遍前缀和就可以得到答案
当时我本来想一边筛质数一边求前缀和,但是如果外层对的
循环不进行开方处理,内层对
的
循环就会在赋初值的时候
溢出导致程序崩溃出错
#include<cstdio>
#include<cstring>
#include<cmath>
const int maxn=1000005;
int vis[maxn],sum[maxn];
void init()
{
int i;
memset(vis,0,sizeof(vis));
memset(sum,0,sizeof(sum));
vis[0]=vis[1]=1;//非质数
for(i=2;i<=sqrt(maxn);++i){
if(!vis[i]){//质数
sum[i]=sum[i-1]+1;//增加
for(int j=i*i;j<maxn;j=j+i){
vis[j]=1;//非质数
}
}else sum[i]=sum[i-1];
}
while(i<maxn){
if(!vis[i])
sum[i]=sum[i-1]+1;
else
sum[i]=sum[i-1];
i++;
}
return ;
}
int main()
{
int n;
init();
while(~scanf("%d",&n)){
printf("%d\n",sum[n]);
}
return 0;
}
-
B -- 高斯消元
- Description
https://vjudge.net/problem/HihoCoder-1195
- Input
第1行:2个正整数,N,M。表示商品的数量N,组合的数量M。1≤N≤500, N≤M≤2*N
第2..M+1行:N+1个非负整数,第i+1行第j列表示在第i个组合中,商品j的数量a[i][j]。第i+1行第N+1个数表示该组合的售价c[i]。0≤a[i][j]≤10, 0≤c[i]≤10^9
- Output
若没有办法计算出每个商品单独的价格,输出"No solutions"
若可能存在多个不同的结果,输出"Many solutions"
若存在唯一可能的结果,输出N行,每行一个非负整数,第i行表示第i个商品单独的售价。数据保证如果存在唯一解,那么解一定恰好是非负整数解
- Sample Input
2 2
2 1 5
1 2 4
- Sample Output
2
1
- 题目理解
题目已经给出伪代码也没有理解就是照着伪代码直接敲出来的
#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
double a[1050][550],b[1050],ans[550];//a为每个商品的数量,b为总价,ans为最后的单价
const double fbs=1e-6;
int n,m;
int solve()
{
for(int i=1;i<=n;i++){
int flag=0;
for(int j=i;j<=m;j++){ //从第i行开始,找到第i列不等于0的行j
if(a[j][i]){
for(int k=1;k<=n;k++)
swap(a[i][k],a[j][k]); //交换第i行和第j行
swap(b[i],b[j]);
flag=1;
}
if(flag) break;
}
if(!flag){ //若无法找到,则存在多个解
printf("Many solutions\n");
return 0;
}
for(int j=i+1;j<=m;j++){
double div=a[j][i]/a[i][i];
for(int k=i;k<=n;k++){
a[j][k]=a[j][k]-a[i][k]*div;
}
b[j]=b[j]-b[i]*div;
}
}
//检查是否无解,存在0=x的情况
for(int i=1;i<=m;i++){ //注意i<=m,之前写成i<=n
int flag=0;
for(int j=1;j<=n;j++){
if(a[i][j]<-fbs||a[i][j]>fbs){
flag=1;
break;
}
}
if(flag==0&&(b[i]>fbs||b[i]<-fbs)){
printf("No solutions\n");
return 0;
}
}
//存在唯一解,
for(int i=n;i>=1;i--){
for(int j=i+1;j<=n;j++){
b[i]=b[i]-a[i][j]*ans[j];
a[i][j]=0;
}
ans[i]=b[i]/a[i][i];
}
return 1;
}
int main()
{
while(scanf("%d%d",&n,&m)!=EOF){
for(int i=1;i<=m;i++){
for(int j=1;j<=n;j++)
scanf("%lf",&a[i][j]);
scanf("%lf",&b[i]);
}
if(solve()){
for(int i=1;i<=n;i++)
printf("%d\n",(int)(ans[i]+0.5));
}
}
return 0;
}
-
C -- 模线性方程组
- Description
https://vjudge.net/problem/HihoCoder-1303
- Input
第1行:1个正整数, N,2≤N≤1,000
第2..N+1行:2个正整数, 第i+1行表示第i组m,r,2≤m≤20,000,000,0≤r<m
计算过程中尽量使用64位整型
- Output
第1行:1个整数,表示满足要求的最小X,若无解输出-1。答案范围在64位整型内
- Sample Input
3
3 2
5 3
7 2
- Sample Output
23
- 题目理解
同上是对着伪代码敲出来的,只是需要增加两个函数以供调用,一个函数、一个
函数,然后就没有然后了
#include<cstdio>
typedef long long ll;
const int maxn=1005;
int n;
ll m[maxn],r[maxn];
ll gcd(ll a,ll b){
if(!b) return a;
else return gcd(b,a%b);
}
ll exgcd(ll a,ll b,ll &x,ll &y){
if(!b){
x=1,y=0;
return a;
}
ll ans=exgcd(b,a%b,x,y);
ll t=x;
x=y;
y=t-(a/b)*y;
return ans;
}
ll solve(){
ll M,R,d,c,k1,k2;
M=m[0],R=r[0];
for(int i=1;i<n;++i){
d=gcd(M,m[i]);
c=r[i]-R;
if(c%d) return -1;
exgcd(M/d,m[i]/d,k1,k2);
k1=(c/d*k1)%(m[i]/d);
R=R+k1*M;
M=M/d*m[i];
R%=M;
}
if(R<0)
R+=M;
return R;
}
int main()
{
while(~scanf("%d",&n)){
for(int i=0;i<n;++i)
scanf("%lld%lld",&m[i],&r[i]);
printf("%lld\n",solve());
}
return 0;
}
-
D -- 中国剩余定理
- Description
https://www.51nod.com/onlineJudge/questionCode.html#!problemId=1079
- Input
第1行:1个数N表示后面输入的质数及模的数量。(2 <= N <= 10) 第2 - N + 1行,每行2个数P和M,中间用空格分隔,P是质数,M是K % P的结果。(2 <= P <= 100, 0 <= K < P)
- Output
输出符合条件的最小的K。数据中所有K均小于10^9
- Sample Input
3
2 1
3 2
5 3
- Sample Output
23
- 题目理解
这道题没有很好的理解这个定理但是C题的代码也可以过
-
F -- Reading comprehension
- Description
Read the program below carefully then answer the question.
#pragma comment(linker, "/STACK:1024000000,1024000000")
#include <cstdio>
#include<iostream>
#include <cstring>
#include <cmath>
#include <algorithm>
#include<vector>
const int MAX=100000*2;
const int INF=1e9;
int main()
{
int n,m,ans,i;
while(scanf("%d%d",&n,&m)!=EOF)
{
ans=0;
for(i=1;i<=n;i++)
{
if(i&1)ans=(ans*2+1)%m;
else ans=ans*2%m;
}
printf("%d\n",ans);
}
return 0;
- Input
Multi test cases,each line will contain two integers n and m. Process to end of file.
[Technical Specification]
1<=n, m <= 1000000000
- Output
For each case,output an integer,represents the output of above program.
- Sample Input
1 10
3 100
- Sample Output
1
5
- 题目理解
通过找规律不难得出、易得
、
将2的幂次方消去得到
,但是这里有一个问题,如果直接进行递推的话数组存不下这么多结果而且时间复杂度也跑不过,于是我们构建矩阵快速幂得到结果,将时间复杂度进行
的降幂
//矩阵快速幂不保留中间过程,也很难保持有些过程并不经过
//指数级降幂
#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
using namespace std;
#define LL __int64
const int maxn=3;
struct Matrix{
LL m[5][5];
}I,A,B,T;
LL a,b,n,mod;
Matrix Mul(Matrix a,Matrix b)
{
Matrix c;
for (int i=1;i<=maxn;i++){
for(int j=1;j<=maxn;j++){
c.m[i][j]=0;
for(int k=1;k<=maxn;k++){
c.m[i][j]+=(a.m[i][k]*b.m[k][j]);
c.m[i][j]%=mod;
}
}
}
return c;
}
Matrix quickpagow(int n)
{
Matrix m=A, b=I;
while(n){
if(n&1)
b=Mul(b,m);
n=n>>1;
m=Mul(m,m);
}
return b;
}
int main()
{
while(~scanf("%I64d%I64d",&n,&mod)){
memset(I.m,0,sizeof(I.m));
memset(A.m,0,sizeof(A.m));
memset(B.m,0,sizeof(B.m));
for(int i=1;i<=maxn;i++){
I.m[i][i]=1;
}
B.m[1][1]=1,B.m[1][2]=2,B.m[1][3]=1;
A.m[1][2]=2;
A.m[2][1]=A.m[2][2]=A.m[3][2]=A.m[3][3]=1;
if(n==1){
printf("%I64d\n",1%mod);
continue;
}
if(n==2){
printf("%I64d\n",2%mod);
continue;
}
T=quickpagow(n-2);
T=Mul(B,T);
printf("%I64d\n",T.m[1][2]%mod);
}
return 0;
}