再加个数学专题

先是矩快吧
题先是两道模拟赛的矩阵快速幂优化dp
先开个坑等补
目前写法最好的：

点击查看代码

#include<bits/stdc++.h>
#define int long long
using namespace std;
const int N=2005,M=105,inf=1e18,mod=1092515507;
int n,m,w[N],v[N],p;
struct T{int fi,se;inline T(int f=-inf,int s=0){fi=f,se=s;}inline T operator*(const T &b){return T{fi+b.fi,se*b.se%mod};}inline void operator+=(const T &b){if(b.fi>fi)*this=b;else if(b.fi==fi) {se+=b.se;if(se>mod)se-=mod;}}
};
struct TT{T a[M][M];int r,c;inline T *operator[](int id){return a[id];}inline const T *operator[](int id)const{return a[id];}inline void init(int x,int y){r=x,c=y;}inline TT operator*(const TT &b){TT res;res.init(r,b.c);for(int k=1;k<=c;k++)for(int i=1;i<=r;i++)for(int j=1;j<=b.c;j++)res[i][j]+=a[i][k]*b[k][j];return res;}
}op[40];
T pow1(int y){TT x;x.init(100,1);x[1][1]=T{0,1};int now=0;while(y){if(y&1)x=op[now]*x;//到底谁乘谁要注意好，不然时间会爆炸的哦～y>>=1,now++;}return x[1][1];
}
signed main(){freopen("pack.in","r",stdin);freopen("pack.out","w",stdout);ios::sync_with_stdio(false);cin.tie(0),cout.tie(0);cin>>n>>m;op[0].init(100,100);for(int i=1;i<=n;i++)cin>>w[i];for(int i=1;i<=n;i++)cin>>v[i],op[0][1][w[i]]+=T{v[i],1};for(int i=2;i<=100;i++)op[0][i][i-1]+=T{0,1};for(int i=1;i<=35;i++)op[i]=op[i-1]*op[i-1];while(m--){cin>>p;T f=pow1(p);// f[1]+={0,1};// for(int k=0;k<=35;k++){//     if(!((p>>k)&1))continue;//     for(int i=1;i<=100;i++){//         for(int j=1;j<=100;j++){//             g[i]+=f[j]*op[k][i][j];//         }//     }//     for(int i=1;i<=100;i++)swap(g[i],f[i]),g[i]=T{};// }if(f.fi<0)cout<<-1<<" "<<-1<<endl;else cout<<f.fi<<" "<<f.se<<endl;}
}