[luogu7116]微信步数

先判定无解，当且仅当存在一个位置使得移动$n$步后没有结束且仍在原地

暴力枚举移动的步数，记$S_{i}$为移动$i$步（后）未离开范围的点个数，则恰好移动$i$步的人数为$S_{i-1}-S_{i}$（特别的$S_{0}=P$），答案即为$sum_{i=1}^{D}(S_{i-1}-S_{i})i=sum_{i=0}^{D-1}S_{i}$（其中$D=min_{S_{i}=0}i$）

考虑如何求出$S_{i}$，记$len_{j}$为第$j$维的合法范围，由于每一维互不干扰，则有$S_{i}=prod_{j=1}^{m}len_{j}$

由于每$n$步必然重复向一个方向移动，那么$Dle nmax w_{i}$，总复杂度即$o(nmmax w_{i})$

进一步优化，将$S_{i}$按$i mod n$的值分组，对于第$i$组（$iin [0,n)$），$i$单独计算，然后求出$i+n$时的$len_{j}$，记$n$步的总位移为$d_{i}$，则第$kn+i(kge 1)$步时$len'_{j}=max(len_{j}-(k-1)|d_{j}|,0)$

写成式子，即$sum_{k=1}^{D'}prod_{j=1}^{m}(len_{j}+|d_{j}|-k|d_{j}|)$（其中$D'=min_{j=1}^{m}frac{len_{j}}{|d_{j}|}$，因为$max$取0后一定为0，不需要再计算）

将后式$o(m^{2})$暴力展开，那么就是一个关于$k$的一个$m+1$次多项式，多项式的前缀和也可以$o(m^{2})$计算，总复杂度即为$o(nm^{2})$

#include<bits/stdc++.h>
using namespace std;
#define N 500005
#define M 20
#define mod 1000000007
int n,m,ans,w[M],x[N],y[N],d[M],mx[N<<1][M],mn[N<<1][M],g[M],sum[M][M];
int ksm(int n,int m){
    if (!m)return 1;
    int s=ksm(n,m>>1);
    s=1LL*s*s%mod;
    if (m&1)s=1LL*s*n%mod;
    return s;
}
int mul(int x,int y){
    for(int i=m+1;i;i--)g[i]=(1LL*g[i-1]*x+1LL*g[i]*y)%mod;
    g[0]=1LL*g[0]*y%mod;
}
int calc(int k,int x){
    int s=1,ans=0;
    for(int i=0;i<=k+1;i++){
        ans=(ans+1LL*s*sum[k][i])%mod;
        s=1LL*s*x%mod;
    }
    return ans;
}
int main(){
    scanf("%d%d",&n,&m);
    sum[0][1]=1;
    for(int i=1;i<=m;i++){//i+1次多项式， 
        int y=0;
        for(int x=0;x<=i+1;x++){
            y=(y+ksm(x,i))%mod;
            memset(g,0,sizeof(g));
            g[0]=y;
            for(int j=0;j<=i+1;j++)
                if (j!=x){
                    mul(1,mod-j);
                    mul(0,ksm((x-j+mod)%mod,mod-2));
                }
            for(int j=0;j<=i+1;j++)sum[i][j]=(sum[i][j]+g[j])%mod;
        }
    }
    for(int i=1;i<=m;i++)scanf("%d",&w[i]);
    for(int i=1;i<=n;i++){
        scanf("%d%d",&x[i],&y[i]);
        d[x[i]]+=y[i];
        for(int j=1;j<=m;j++){
            mx[i][j]=max(mx[i-1][j],d[j]);
            mn[i][j]=min(mn[i-1][j],d[j]);
        }
    }
    for(int i=n+1;i<=2*n;i++){
        d[x[i-n]]+=y[i-n];
        for(int j=1;j<=m;j++){
            mx[i][j]=max(mx[i-1][j],d[j]);
            mn[i][j]=min(mn[i-1][j],d[j]);
        }
    }
    for(int i=1;i<=m;i++)d[i]/=2;
    bool flag1=0,flag2=0;
    for(int i=1;i<=m;i++)
        if (d[i])flag1=1;
    for(int i=1;i<=m;i++)
        if (mx[n][i]-mn[n][i]>w[i])flag2=1;
    //flag1说明离开了原位置，flag2说明一定是一轮 
    if ((!flag1)&&(!flag2)){
        printf("-1");
        return 0;
    }
    for(int i=0;i<n;i++){
        int s=1;
        for(int j=1;j<=m;j++)s=1LL*s*max(w[j]-(mx[i][j]-mn[i][j]),0)%mod;
        ans=(ans+s)%mod;
    }
    for(int i=n;i<2*n;i++){
        bool flag=0;
        for(int j=1;j<=m;j++)
            if (w[j]-(mx[i][j]-mn[i][j])<0)flag=1;
        if (flag)break;
        memset(g,0,sizeof(g));
        g[0]=1;
        for(int j=1;j<=m;j++)mul(mod-abs(d[j]),w[j]-(mx[i][j]-mn[i][j])+abs(d[j]));
        int D=0x3f3f3f3f;
        for(int j=1;j<=m;j++)
            if (d[j])D=min(D,(w[j]-(mx[i][j]-mn[i][j]))/abs(d[j])+1);
        for(int j=0;j<=m;j++)ans=(ans+1LL*g[j]*calc(j,D))%mod;
    }
    printf("%d",ans);
}