数论基础训练(2)-解题报告

A-排列

又一次倒在stl-next_permutation函数了，虽然以前用过，但没明白函数运行的原理，原来它是先按字典序排列的，在从后面一次排列，交换顺序的，而且它可以再你给定的一个排列的位置，可以从现在的序列接着从后排列，当排列完的时候，它能回到第一个序列(即为最少字典序)，而这题正是这么来着，无语了，这个函数一定要记住它。

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std ;

int a[1300];

int main()
{
	int t,n,k;
    scanf("%d",&t);
    while(t--){
		scanf("%d %d",&n,&k);
		for(int i=0;i<n;i++)
			scanf("%d",&a[i]);
		for(int i=0;i<k;i++)
			next_permutation(a,a+n);
		for(int i=0;i<n;i++)
			printf(i==0?"%d":" %d",a[i]);
		puts("");
    }
 	return 0 ;
}

C题：这应该不是我们做的！代码也赖得贴了，01背包的大水题，基础的不能再基础了，可惜就是WA几次了，原因，数组开小了。

D题：这题不能暴搞，显然这题是有规律的，取模存在循环节，如果后面取余与前面相同，就没必要判断了。

#include <iostream>
#include <cstdio>
#include <cstring>
using namespace std;
bool flag[100000];
int main()
{
    int n,sum,i;
    while(scanf("%d",&n)!=EOF)
    {
        if((n&1)==0||n==1){
            printf("2^? mod %d = 1\n",n);
            continue;
        }
        i=1,sum=2;
        memset(flag,0,sizeof(flag));
        while(sum%n!=1)
        {
            sum<<=1;
            sum=sum%n;
            i++;
            if(sum==2||flag[sum]) break;
            flag[sum]=1;
        }
        if(sum%n!=1) printf("2^? mod %d = 1\n",n);  
        else printf("2^%d mod %d = 1\n",i,n);
          
    }
    return 0;
}

E-world cup noise

题意：给出有n个0，1的字串长度，当出现两个11在一起，就结束，问最多有多小个排序方案。

很明显的规律题当n=1  ans=2 ,   n=2   ans=3   n=3 ans=5 ,n=4 ans=8;---------斐波拉契数列，就是n=1是2，递推关系还不是一样的。

找规律题就是多枚举几种情况，加上直觉，瞬间可破！

代码没什么可贴的，预先处理到50就行了。

F-Geometry Made Simple

没什么可解释的，完全来增加信心的。

#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
using namespace std ;

int main()
{
	double a[3];
	int t=0,pos;
    while(1){
		t++;
		for(int i=0;i<3;i++){
			scanf("%lf",&a[i]);
			if(a[i]==-1) pos=i;
		}
		if(a[0]==0&&a[1]==0&&a[2]==0) break;
		printf("Triangle #%d\n",t);
		if(pos==0){
			a[0]=a[2]*a[2]-a[1]*a[1];
			if(a[0]<=0) printf("Impossible.\n\n");
			else printf("a = %.3lf\n\n",sqrt(a[0]));
		}
		else if(pos==1){
			a[1]=a[2]*a[2]-a[0]*a[0];
			if(a[1]<=0) printf("Impossible.\n\n");
			else printf("b = %.3lf\n\n",sqrt(a[1]));
		}
		else {
			a[2]=a[0]*a[0]+a[1]*a[1];
			printf("c = %.3lf\n\n",sqrt(a[2]));
		}
    }
 	return 0 ;
}

G题：解高次同余方程。

不会，看看模板吧！

H-C Looooops

题意：

for (variable = A; variable != B; variable += C)

 statement；

题目还说：当数过大时，在计算机存取时会出现循环，而计算机是用二进制表示的，问循环多小次退出循环。

故有:(A+x*C)%(2^k)=B%(2^k)====>CX+y*(2^k)=B-A;其中X代表循环次数，y为一个辅助整数。

看到这个式子很快就想到了扩展欧几里得，求线性同余方程的最小解（X>0），再基础不过。

#include<iostream>
#include<cstdio>
#include<cstring>
typedef __int64 LL;
using namespace std ;

LL exgcd(LL a,LL b,LL &x,LL &y)
{
	if(b==0){
		x=1;
		y=0;
		return a;
	}
	LL d=exgcd(b,a%b,x,y);
	LL t=x-a/b*y;
	x=y;
	y=t;
	return d;	 
}

int main()
{
	LL a,b,c,x,y;
    LL K;
	int k;
     while(scanf("%I64d %I64d %I64d %d",&a,&b,&c,&k),(a||b||c||k)){
		    K=1;
			for(int i=0;i<k;i++){   //这个有点SB了，直接K=1<<k,比赛时急了，果断暴搞。
			     K<<=1;
			}
			LL d=exgcd(c,K,x,y);
			LL C=b-a;
			if(C%d){
				puts("FOREVER");
				continue;
			}
			C/=d,K/=d;
			LL ans=((x*C)%K+K)%K;
			printf("%I64d\n",ans);	
     }
 	return 0 ;
}

I-Goldbach's Conjecture

没看题目，就知道是干什么的，哥德巴赫猜想，你我皆知，求它满足的方案数。

先素数筛选，然后暴力。

#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
using namespace std ;
const int N=100001;
bool vis[N];

void getprime()
{
	for(int i=4;i<=N;i+=2){
		vis[i]=1;
	}
	int m=sqrt(N);
	for(int i=3;i<=m;i+=2){
		if(vis[i]) continue;
		for(int j=i*i;j<=N;j+=2*i)
		     vis[j]=1;
	}
}
int main()
{
	int n,t,ans;
    getprime();
    while(scanf("%d",&n),n){
		t=n/2;
		ans=0;
		for(int i=2;i<=t;i++){
			if(!vis[i]&&!vis[n-i]) ans++;
		}
		printf("%d\n",ans);
    }
 	return 0 ;
}