2015暑假第一次积分赛

题目地址

A Little Keng（快速求幂）

题目大意：计算1n + 2n + 3n +
... + mn 的值的末尾有多少个0，其中(1 <= m <= 100 , 1 <= n <= 1000000)
解题思路：看到mn翻书找数论的知识，发现和快速求幂很像，想了一下那个mod应该是不能超过int，就用了平时最喜欢用的1000000000，具体原因不知道，直接套用了模板

#include <iostream>
#include <cstdio>
#include <cmath>
#include <algorithm>
#include <cstring>
#include <queue>
#include <string>

using namespace std;

const int mod = 1000000000;

long long f(int m,int n)
{
if(n == 1)  return m;
long long ret = f(m,n/2);
ret = (ret*ret) % mod;
if(n&1)
ret = (ret*m) % mod;
return ret;
}

int main()
{
int n, m;
while(scanf("%d%d",&m,&n) != EOF)
{
long long sum = 0;
for(int i = 1; i <= m; i++)
{
sum = sum + f(i,n);
}
int ans = 0;
while(sum%10 == 0 && sum != 0)
{
ans++;
sum /= 10;
}
printf("%d\n", ans);
}
return 0;
}<span style="font-family:SimSun;font-size:10px;"><strong>
</strong></span>

B Big Keng
题目大意：给定一个立体的图形，上面是圆柱，下面是圆台，圆柱的底面半径和圆台的上半径相等，然后体积的V时，问这个图形的表面积最小可以是多少。（不算上表面）。据说是三分，没想出怎么做

C  Bloodsucker（概率DP）

题目大意：第0天有n-1个人，1个吸血鬼，每天只能有2个生命相遇，若为同类则没反应，若为1人1鬼，则人变成鬼的概率为p，给出n和p，求所有人变成吸血鬼的期望

解题思路：逆推，dp[i]代表i个人到所有人变成吸血鬼的概率，往上逆推，则：dp[i]=(dp[i+1]+1)*p1+(dp[i]+1)*（1-p1）

移项后化简得： p1*dp[i]=dp[i+1]*p1+1

#include <iostream>
#include <cstdio>
#include <cmath>
#include <algorithm>
#include <cstring>
#include <queue>
#include <string>

using namespace std;

double dp[100005];

int main()
{
int t,n;
double p;
scanf("%d",&t);
while(t--)
{
scanf("%d%lf",&n,&p);
dp
= 0;
for(int i = n-1; i >= 1; i--)
{
double s1,s2,p1;
s1 = (double)n*(n-1)/2;//从n个人中选两个人出来的选法
s2 = (double)i*(n-i);//为从i个人中选出一人的方法，n-i为从吸血鬼中选出一个吸血鬼的总数，共有这么多种
p1 = s2/s1*p;//人与吸血鬼相遇的概率
dp[i] = (dp[i+1]*p1+1)/p1;
}
printf("%.3f\n",dp[1]);
}
return 0;
}

H How Many Sets I（容斥原理+快速求幂）

题目大意：一个有n个元素的集合S，其子集为Si， 求有多少个有序的(S1, S2, ..., Sk)
组合满足S1 ∩ S2 ∩ ... ∩ Sk =
∅，题目给出n和k
解题思路：总情况数减去有交集的子集数。
2^n^k - C(n, 1)*2^(n-1)^k + C(n, 2)*2^(n-2)^k ... C(n, n)*2^(n-n)^k 化简为(2^k - 1)^n，再 两次快速幂

#include <iostream>
#include <cstdio>
#include <cmath>
#include <algorithm>
#include <cstring>
#include <queue>
#include <string>

using namespace std;

const int mod = 1000000007;

int powMod(long long a,int b)
{
long long res=1;
while(b)
{
if(b&1)
res = (res*a)%mod;
b >>= 1;
a = (a*a)%mod;
}
return res;
}

int main()
{
int n,k;
while(scanf("%d%d",&n,&k) != EOF)
printf("%d\n", powMod((powMod(2, k)-1+mod)%mod, n));
return 0;
}