2019 ICPC 沈阳网络赛 K. Guanguan's Happy water

2019 ICPC 沈阳网络赛 K. Guanguan's Happy water


题目大意:已知一个数列f(n):
f[x]=a[x] (1<=x<=k)
f[x]=f[x-1]*p[1]+f[x-2]*p[2]………+f[x-k]*p[k] (x>k)

给你所有的a[i],再给你接下来k个f[i],求f这个数列的前n项和。

解法:第一名10分钟以内AC,是人性的扭曲还是道德的沦丧?高斯消元加等比矩阵求和,我们组写了半天,WA,TLE。直接放弃。
赛后同校队伍:如果n<=k,直接求和输出,如果n>k,先求出前k个的和,接下来每k个为一个循环节,答案直接为sum[k]+循环节和*(n/k-1)+f[k+1]+…+f[n%k+k]。
nmd,wsm??
我打表打了半天,公式写了这么多,你和我说就这么过了????!!!!
我随便找一组数据就把你hack了,这种代码都敢提交,怕不是没读题
然后我交了一发,这是我的表情:
在这里插入图片描述
下面是错误的AC代码:

#include <bits/stdc++.h>
using namespace std;
using namespace chrono;
const int N = 1005;
const int M = 1000000007;
const int INF = 0x3f3f3f3f;
const double PI = acos(-1);
const double eps = 1e-8;
#define ms(x, y) memset((x), (y), sizeof(x))
#define mc(x, y) memcpy((x), (y), sizeof(y))
typedef long long ll;
typedef unsigned long long ull;
#define fi first
#define se second
#define mp make_pair
typedef pair<int, int> pii;
typedef pair<ll, int> pli;
#define bg begin
#define ed end
#define pb push_back
#define al(x) (x).bg(), (x).ed()
#define st(x) sort(al(x))
#define un(x) (x).erase(unique(al(x)), (x).ed())
#define fd(x, y) (lower_bound(al(x), (y)) - (x).bg() + 1)
#define ls(x) ((x) << 1)
#define rs(x) (ls(x) | 1)
template <class T>
bool read(T & x) {
    char c;
    while (!isdigit(c = getchar()) && c != '-' && c != EOF);
    if (c == EOF) return false;
    T flag = 1;
    if (c == '-') { flag = -1; x = 0; } else x = c - '0';
    while (isdigit(c = getchar())) x = x * 10 + c - '0';
    x *= flag;
    return true;
}
template <class T, class ...R>
bool read(T & a, R & ...b) {
    if (!read(a)) return false;
    return read(b...);
}
mt19937 gen(steady_clock::now().time_since_epoch().count());
int a[N];
int main()
{
    time_point<steady_clock> start = steady_clock::now();

    int T, k;
    ll n;
    read(T);
    while (T--) {
        read(k, n);
        for (int i = 1; i <= 2 * k; i++) read(a[i]);
        ll sum = 0;
        for (int i = 1; i <= min(1LL * k, n); i++) (sum += a[i]) %= M;
        if (n / k > 1) {
            ll tmp = 0;
            for (int i = k + 1; i <= 2 * k; i++) (tmp += a[i]) %= M;
            ll kk=((n/k)-1)%M;
            sum = (sum + tmp * kk%M) % M;
        }
        if (n / k > 0) {
            for (int i = 1; i <= n % k; i++)
                (sum += a[i + k]) %= M;
        }
        printf("%lld\n", sum);
    }

   // cerr << endl << "------------------------------" << endl << "Time: "
    //     << duration<double, milli>(steady_clock::now() - start).count()
     //    << " ms." << endl;

    return 0;
}

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值