XDU暑期训练D1训练日志

第一天讲了组合数学和概率期望，懂了有一半。

下午的训练也没写出来多少题，就把我写出来的题做以总结。

A题 整数分解为2的次幂（求出组合总数）

 

这道是PPT的例题，大概思路是

当n为偶数时：

      n的总数是由（n-1）的总数+（n/2）的总数组成；/*（n-1）的每种方案中添加1，

即是所有包含1的方案；剩下的方案同时/2即是（n/2）的情况.*/

当n为奇数时：

      n的总数是（n-1）的总数；//因为只多了1，所以在（n-1）的每种方案中添加1即可.

代码如下:

#include<cstdio>
using namespace std;
const int N=1000000007;
const int M=1100000; 
long long a[M];//记忆化储存
long long f(long long n)
{
	if(n<=1) return 1;
	if(a[n]) return a[n];
	if(n&1)
	{
		a[n]=f(n-1)%N;
		return a[n];
	}
	a[n]=(f(n-1)%N+f(n>>1)%N)%N;
	return a[n];
 } 
int main(){
	long long n;
	scanf("%lld",&n);
	printf("%lld",f(n));
	return 0;
} 

注：来源 51Nod - 1383。

B题 Permutations

 

题意大概是：随意给出一个n的全排列，求出它要经过几次置换才能成为单位元。

大致如上图所示。

本题可以换一个思路来做，去找每一位的数的循环节有几个元素（如1-4-2-1，则1的循环节

有3个元素），最后求出它们的最小公倍数输出即可。

代码如下：

#include<stdio.h>
#include<string.h>
using namespace std;
int gcd(int a,int b)//求最大公约数
{
	if(b==0) return a;
	else return gcd(b,a%b);
}
int lcm(int a,int b)//求最小公倍数
{
	return a/gcd(a,b)*b;//lcm(a,b)=a/gcd(a,b)*b
}
int a[1100],map[1100];
int main()
{
	int n,ans=1;
	scanf("%d",&n);
	memset(map,0,sizeof(map));
	for(int i=1;i<=n;i++)
		scanf("%d",&a[i]);
	for(int i=1;i<=n;i++)
	{
		int temp=a[i];
		int count=1;
		map[i]=1;
		while(temp!=i&&!map[temp])
		{
			map[temp]=1;
			temp=a[temp];
			count++;
		}//寻找循环节
		ans=lcm(ans,count);
	}
	printf("%d\n",ans);
}

注：其实这一道题是关于置换群的基础题，来源poj 2369。

C题 Invoker

题意：c种颜色，n个小球的涂色问题，问有几种染色方案（经过旋转翻折可以相同的算是一种方案）

老实说这道题是一窍不通，抄了模板过的。

大概是分为旋转和翻折两种情况去判断。

代码如下：

#include<stdio.h>
typedef long long ll;
const ll mod=1000000007;
ll c,n;
ll gcd(ll a,ll b)
{
	return b==0?a:gcd(b,a%b);
}
ll pow(ll p,ll n)//快速幂(含取模)
{
	ll ans=1;
	while(n)
	{
		if(n&1) ans=ans*p%mod;
		p=p*p%mod;
		n/=2;
	}
	return ans; 
}
//拓展欧几里得，一直没搞明白，解方程x*a+y*b=1,求逆元
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;
	x=y;
	y=t-a/b*y;
	return d;
}
int main()
{
	int t;
	scanf("%d",&t);
	int cas=1;
	while(t--)
	{
		scanf("%lld%lld",&c,&n);
		int ans=0;
		for(ll i=1;i<=n;i++)//旋转的情况
		{
			ans+=pow(c,gcd(n,i));//sum(c^gcd(n,i))(i=1,2,~~n)
			ans%=mod;
		}
                //翻折的情况
		if(n&1) ans+=(n*pow(c,n/2+1))%mod;
                //n为奇数时, n*c^(n/2+1)(此处除法向下取整)
		else ans+=(n/2*pow(c,n/2+1))%mod+(n/2*pow(c,n/2)%mod)%mod;
                //n为偶数时, n*c^(n/2+1)+n/2*c^(n/2)(此处除法向下取整)
		ans%=mod;
		ll x,y;//x为mod的逆元， 
		exgcd(2*n,mod,x,y);
		x=(x+mod)%mod;
		printf("Case #%d: %lld\n",cas++,ans*x%mod);	
	}
}

 

I题： A * B Problem Plus HDU - 1402 

这题就是模拟乘法，数据比较大，用了JAVA的大数类。

代码如下：

import java.math.BigInteger;
import java.util.Scanner;

public class Main {
	public static void main(String args[])
	{
		String str1,str2;
		Scanner in=new Scanner(System.in);
		while(in.hasNext())
		{
			str1=in.next();
			str2=in.next();
			BigInteger a=new BigInteger(str1);
			BigInteger b=new BigInteger(str2);
			System.out.println(a.multiply(b));
		}
	}
}

注：不得不说Java好强！