2019 ICPC 南昌区域赛 - G Eating Plan（技巧+暴力）

题目链接

---

题目大意

给出大小为$n$的序列，是$n$的排列，对于$m$个询问每次给出一个数$k$，求出从序列中选出最少为多少个数，使得区间内每个数的阶乘的和对998857459取模的结果大于等于$k$。

解题思路

这题有个不容易发现的$trick$，取模的998857459是个合数！对它分解质因子后会发现它最大的那个质因子是2803，这题首先要对初始数组进行阶乘处理，很明显大于等于2803的数，都是998857459的倍数，阶乘后再取模都是0，因此只要求前面2802个阶乘即可

然后暴力写了就是了，用双重循环求每个连续的长度对应的最小和，最后对每个查询线性找答案即可。

代码：

#include <iostream>
#include <cstring>
#include <cstdio>
using namespace std;
typedef long long ll;
const int N=1e5+5;
const int mod=998857459;
const int inf=0x3f3f3f3f;
int cnt;
struct node{
    ll val;
    int id;
}res[2900];
ll sum[2900],ans[N],jc[N],a[N],n,m;
int main()
{
    memset(a,0,sizeof(a));
    memset(ans,0,sizeof(ans));
    jc[0]=1;
    for(ll i=1;i<2803;i++)
    jc[i]=(jc[i-1]*i)%mod;
    scanf("%d%d",&n,&m);
    for(int i=1;i<=n;i++){
        int x;
        scanf("%d",&x);
        a[i]=jc[x];
    }
    sum[0]=0;
    for(int i=1;i<=n;i++){
        if(a[i]){
            res[++cnt].val=a[i];
            res[cnt].id=i;
            sum[cnt]=sum[cnt-1]+res[cnt].val%mod;
        }
    }
    ///暴力求出每个区间长度的最小子段和
    for(int i=1;i<=cnt;i++){
        for(int j=i;j<=cnt;j++){
            ans[res[j].id-res[i].id+1]=max(ans[res[j].id-res[i].id+1],(sum[j]-sum[i-1])%mod);
        }
    }
    while(m--){
        ll x;
        scanf("%lld",&x);
        int len=inf;
        for(int i=1;i<=n;i++)
        if(ans[i]>=x){
            len=i;
            break;
        }
        if(len==inf)
            printf("-1\n");
        else
            printf("%d\n",len);
    }
    return 0;
}