2019ICPC南昌-G - Eating Plan题解

2019ICPC南昌-G - Eating Plan题解

题目链接：
https://nanti.jisuanke.com/t/42582

题面：

题目大意：
给定n(1e5)个数,m(1e4)次询问，对于n个数（1到n的排列）,每个数的权值是它阶乘值%998857459，每次询问一个数k，找到区间 使得这个区间内的所有值的权值和%998857459是大于等于k的，输出最小的区间的长度。

我们可以发现998857459这个神奇的数并不是素数，它的最大因子除了自己是2803，也就是说只有小于2803的数才会有权值，其他的权值都是0（因为要取模998857459），所以我们可以用一个vector存所有的小于2803的数，然后对这些数跑一个n^2来得到一个lim数组，lim [ i ] 表示区间长度位i的最大可以满足的k，最后对lim数组进行从1到n的一个从小到大的覆盖，即:
if(lim[i-1]>lim[i])：lim[i] = lim[i-1];

这样子我们就能二分得到答案了。

#include<bits/stdc++.h>
#define ll long long
#define ull unsigned long long
#define mem(a,b) memset(a,b,sizeof(a))
#define inf 0x3f3f3f3f
#define IO std::ios::sync_with_stdio(false);cin.tie(0)
#define endl '\n'
#define lb(x) x&(-x)
#define pii pair<int,int>
#define pll pair<ll,ll>
using namespace std;
const double pi = 3.141592653589793238462;
const int maxn = 1e5+5;
const double eps = 1e-18;
const ull base  = 13331;
const ll mod = 998857459;


ll bk[maxn];
int a[maxn];
int ne[maxn];
int n,m;

ll lim[maxn];

vector<int> pos;

ll ff[maxn];

void init(){
    bk[0] = 1;
    for(int i = 1;i<=2900;i++){
        bk[i] = bk[i-1] * i;
        bk[i]%=mod;
    }
}

int solve(ll x){
    int l = 1,r = n;
    int ans = -1;
    while(l<=r){
        int mid = (l+r)>>1;
        if( lim[mid] >= x ){
            ans = mid;
            r--;
        }
        else l++;
    }
    return ans;
}

int main(){
    IO;
    cin>>n>>m;
    init();
    for(int i = 1;i<=n;i++)cin>>a[i];
    for(int i = 1;i<=m;i++){
        cin>>ff[i];
    }
    for(int i = 1;i<=n;i++){
        if(a[i]<2803){
            pos.push_back(i);
        }
    }

    for(int i = 0;i<pos.size();i++){
        int l = pos[i];
        ll sum = bk[a[l]];
        lim[1] = max(lim[1],sum);
        for(int j = i+1;j<pos.size();j++){
            int r = pos[j];
            sum+=bk[a[r]];
            sum%=mod;
            int cb = r-l+1;
            lim[cb] = max(lim[cb],sum);
        }
    }
    for(int i = 2;i<=n;i++){
        lim[i] = max(lim[i],lim[i-1]);
    }
    for(int i = 1;i<=m;i++){
        ll v = ff[i];
        ne[i] = solve(v);
    }
    for(int i = 1;i<=m;i++){
        cout<<ne[i]<<endl;
    }
}