Time Limit:5s Memory Limit:256MByte
Submissions:1597Solved:301
今天HHHH考完了期末考试,他在教学楼里闲逛,他看着教学楼里一间间的教室,于是开始思考:
如果从一个坐标为 (x1,y1,z1)(x1,y1,z1)的教室走到(x2,y2,z2)(x2,y2,z2)的距离为|x1−x2|+|y1−y2|+|z1−z2||x1−x2|+|y1−y2|+|z1−z2|
那么有多少对教室之间的距离是不超过RR的呢?
#include <iostream>
#include <cstdio>
#include <algorithm>
#include <cstring>
using namespace std;
typedef long long ll;
int cnt[15][15][15], book[15][15][15];
ll dp[50];
int main()
{
int t, n, q, x, y, z;
cin >> t;
while(t--)
{
memset(dp, 0, sizeof(dp));
memset(cnt, 0, sizeof(cnt));
scanf("%d%d", &n,&q);
for(int i = 1; i <= n; i++)
scanf("%d%d%d", &x, &y, &z), cnt[x][y][z]++;
for(int i = 0; i <= 10; i++)
for(int j = 0; j <= 10; j++)
for(int k = 0; k <= 10; k++)
for(int x = 0; x <= 10; x++)
for(int y = 0; y <= 10; y++)
for(int z = 0; z <= 10; z++)
{
int dis = abs(i-x) + abs(j-y) + abs(k-z);
if(i == x && j == y && k == z)
{
dp[dis] = dp[dis] + (ll)(cnt[i][j][k]-1)*(cnt[i][j][k])/2;
}
else dp[dis] = dp[dis] + (ll)(cnt[i][j][k]*cnt[x][y][z]);
}
for(int i = 1; i <= 30; i++)
dp[i] = dp[i]/2 + dp[i-1];
while(q--)
{
scanf("%d", &x);
if(x >= 30) x = 30;
printf("%lld\n", dp[x]);
}
}
return 0;
}
Time Limit:5s Memory Limit:256MByte
Submissions:660Solved:160
今天HHHH在操场上跑步,HHHH作为一个爱运动的人,肯定会想方设法把跑步所消耗的能量减到最少.
现在这个操场上有nn个可以休息的点,他们的坐标分别为x1,x2...xn(xi≤xi+1)x1,x2...xn(xi≤xi+1),HHHH刚开始在x1x1 ,并且他只能在这些点休息,在中途不能停下来,否则会因为旁边的音浪太强而被晃到.
如果HHHH连续跑一段长度为ll的距离,那么他将会消耗2l+a2l+a的能量(aa为HHHH的可爱值)
现在给你这些点的坐标,请帮HHHH计算他跑到xnxn点所需要消耗的能量最少是多少.
2x很快就会比a高到不知道哪里去,所以当这个距离超过一定值的时候肯定是分开
跑比较划算了,所以我们只需要枚举当前点前面的几十个点再转移一下就行了
#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;
typedef long long ll;
const ll INF = 1e18;
const int maxn = 1e5 + 5;
int t, n, a;
ll x[maxn], dp[maxn], cnt[50];
int main()
{
for(int i = 0; i <= 30; i++) //预处理
cnt[i] = (1<<i);
scanf("%d", &t);
while(t--)
{
scanf("%d%d", &n, &a);
for(int i = 0; i < maxn; i++)
dp[i] = INF;
for(int i = 1; i <= n; i++)
scanf("%d", &x[i]);
dp[1] = 0;
for(int i = 2; i <= n; i++)
{
for(int j = 1; j < i; j++)
{
if(x[i]-x[i-j] <= 30) //三十个点就大于a的最大值了
dp[i] = min(dp[i], dp[i-j] + cnt[x[i]-x[i-j]]+a);
else break;
}
}
printf("%lld\n", dp[n]);
}
return 0;
}
2l+a 2l+a
Time Limit:5s Memory Limit:256MByte
Submissions:405Solved:85
今天HHHH在学数论,他看到一个很优美的式子:
∑ni=1ik mod p∑i=1nik mod p
一向热衷于抱队友大腿的HHHH便问队友ZZZZ怎么做
ZZZZ:"n,kn,k多大?"
HHHH:"105,105105,105"
ZZZZ:"快速幂嘛"
HHHH:"109,105109,105"
ZZZZ:"拉格朗日插值嘛"
HHHH:"1018,10181018,1018"
队友:"让我想想.."
官方题解:
看到n比p大了不到100个,我们可以思考一下假如我们只算到p或者p−1会是怎么样的
即考虑函数∑p−1i=1ik mod p
我们先打个表能发现当k是p−1的倍数的时候这个值是p−1,别的时候都是0
然而这是否正确的呢? 我们可以简单的推导一下:
假设g是p的原根,那么根据原根的性质,g的1次方到p−1次方在模p的时候
对应了1到p−1的某个全排列,即原式等价于
gk+g2k+g3k+...+g(p−1)k
然后我们用等比数列求和一下则得到
(gk−gpk)/(1−gk)=(gk−gpk)×inv(1−gk)
再根据费马小定理我们便得到了gp=g,那么这个式子就变成了0
不过这是建立在这是等比数列的情况下的,所以我们对几种不是等比数列的情况需要特判:
(1) p为2的时候原根是1
(2) k为p-1的时候由费马小定理可以得到所有项均为1
所以我们现在只需要算剩下的100个左右的数就可以了,快速幂搞一下就行.
不过要小心这个范围乘法会爆long long 所以需要快速乘/O(1)快速乘
①首先我们知道,i^k%p和(i+p)^k%p的价值是一样的。
那么肯定以长度p为循环节,是有循环内容的。
②当k不是p-1的倍数的时候,我们发现,每一段长度为p的内容加和%p都是0.
②当k是p-1的倍数的时候,我们发现,每一段长度为p的内容加和都是p-1.
贯彻“遇事不决,就打表”
代码:
#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;
typedef long long LL;
LL p, k;
LL quick_mul(LL x,LL y,LL MOD)
{
x = x%MOD,y = y%MOD;
return ((x*y-(LL)(((long double)x*y+0.5)/MOD)*MOD)%MOD+MOD)%MOD;
}
LL quick_Mod(LL x, LL n)
{
LL ans = 1;
while(n)
{
if(n&1)
{
ans = quick_mul(ans, x, p);
}
x = quick_mul(x, x, p);
n >>= 1;
}
return ans;
}
int main()
{
LL n;
int t;
scanf("%d", &t);
while(t--)
{
scanf("%lld%lld%lld", &n, &k, &p);
if(k==0)
{
printf("%lld\n",n%p);
continue;
}
LL ans = 0;
if(k%(p-1) == 0)
{
ans = (ans + quick_mul(p-1, n/p, p))%p;
}
n %= p;
for(LL i = 1; i <= n; i++)
{
ans = (ans + quick_Mod(i, k))%p;
}
printf("%lld\n", ans);
}
return 0;
}