2013金山西山居创意游戏程序挑战赛——初赛（2）hdu 4548 4549 4550

A题水题。直接暴力素数表。然后从1开始对美素数进行个数的计算。

c[n]记录从1—n中美素数的个数。

hdu 4548代码

#include <iostream>
#include <stdio.h>
#include <math.h>
#include <stdlib.h>
#include <string.h>
#include <limits.h>
#include <vector>
#include <algorithm>
using namespace std;
bool flag[1000010];
int c[1000010];
void init()
{
    int i,j,k;
    memset(flag,0,sizeof(flag));
    for(i = 2;i < 1000001; i++)
    {
        if(!flag[i])
        {
            for(j = i+i;j < 1000001;j+=i)
                flag[j] = true;
        }
    }
}
bool check(int n)
{
    int tmp = 0;
    while(n)
    {
        tmp += n%10;
        n /= 10;
    }
    if(flag[tmp] == false)    return true;
    else            return false;
}
int main()
{
    init();
    memset(c,0,sizeof(c));
    c[2] = 1;c[3] = 2;
    for(int i = 4;i<=1000000;i++)
    {
        c[i] = c[i-1];
        if(flag[i] == false && check(i))
            c[i]++;
    }
    int n,m;
    int T,cas = 0;;
    scanf("%d",&T);
    while(T--)
    {
        scanf("%d %d",&n,&m);
        printf("Case #%d: %d\n",++cas,c[m] - c[n-1]);
    }
    return 0;
}

B题，容易推算出式子： F[n] = (a^m1)*(b^m2)

设fibonacci数组为fib。则 m1 = fib[n - 1],m2 = fib[n ]，n = 0或1的时候特判定。

费马小定理：

 a^fib[n - 1] = a^(fib[n - 1]%1000000006) (mod 1000000007)

b^fib[n] = b^(fib[n]%1000000006) (mod 1000000007)

剩下的就很简单了。求fib用矩阵快速幂，求答案用整数快速幂。
hdu 4549 代码：

#include <iostream>
#include <stdio.h>
#include <math.h>
#include <stdlib.h>
#include <string.h>
#include <limits.h>
#include <vector>
#include <algorithm>
using namespace std;
#define llg long long
const llg mod = 1000000007;
struct matrix
{
	llg a[2][2];
	void init()
	{
		a[0][0] = a[0][1] = a[1][0] = a[1][1] = 1;
		a[1][1] = 0;
	}
};
matrix mul1(matrix tmp1,matrix tmp2)
{
	matrix ans;
	for(int i = 0;i<2;i++)
		for(int j = 0;j<2;j++)
		{
			ans.a[i][j] = 0;
			for(int k = 0;k<2;k++)
			{
				ans.a[i][j] += (tmp1.a[i][k]*tmp2.a[j][k])%(mod - 1);
				ans.a[i][j] %= (mod - 1);
			}
		}
	return ans;
}
llg mul2(matrix tmp,llg n)
{
	matrix ans;
	ans.init();
	while(n)
	{
		if(n&1) ans = mul1(ans,tmp);
		tmp = mul1(tmp,tmp);
		n>>=1;
	}
	return ans.a[1][1];
}

llg q_mod(llg a,llg n)
{
	llg ans = 1;
	matrix tmp2;
	tmp2.init();
	llg nn = mul2(tmp2,n);
	nn %= (mod - 1);
//cout<<"fib : "<<nn<<endl;	
	llg tmp = a;
	while(nn)
	{
		if(nn&1)	ans *= tmp;
		ans %= mod;
		tmp *= tmp;
		tmp %= mod;
		nn>>=1;
	}
	return ans%mod;
}
int main()
{
	llg a,b,n;
	while(cin >> a >> b >> n)
	{
		if(n == 0)	{cout<<a<<endl;continue;}
		if(n == 1)	{cout<<b<<endl;continue;}
		if(a == 0 || b == 0)	{cout<<0<<endl;continue;}
		a = a%mod;
		b = b%mod;
		n = n%mod;
		llg tmp = q_mod(a,n-1);

		tmp *= q_mod(b,n);
	
		tmp %=mod;
		cout<<tmp<<endl;
	}
	return 0;
}

c题，从网络看别人代码的，很好的贪心，先对最小值进行处理。很精彩的一份代码。

hdu 4550 代码：

#include<iostream>
#include<string.h>
#include<stdio.h>
using namespace std;
#define llg long long
int main()
{
	int T;
	char s[110];
	string a;
	scanf("%d",&T);
	while(T--)
	{
		scanf("%s",s);
		int len = strlen(s);
		a = s[0];
		char minc = '9' + 1;
		int index = 0;
		for(int i = 0;i < len;i++)
		{
			if(s[i] != '0')
			{
				if(minc >= s[i])
				{
					minc = s[i];
					index = i;
				}
			}
		}
		for(int i = 1;i < len;i++)
		{
			if(i == index)	a = s[i] + a;
			else
			{
				if(i > index)	a = a + s[i];
				else
				{
					if(s[i] > a[0])	a = a + s[i];
					else			a = s[i] + a;
				}
			}
		}
		cout<<a<<endl;
	}
	return 0;
}