結果
問題 | No.695 square1001 and Permutation 4 |
ユーザー |
![]() |
提出日時 | 2019-02-15 00:32:23 |
言語 | C++11 (gcc 13.3.0) |
結果 |
AC
|
実行時間 | 2,352 ms / 7,000 ms |
コード長 | 13,890 bytes |
コンパイル時間 | 1,918 ms |
コンパイル使用メモリ | 173,084 KB |
実行使用メモリ | 42,368 KB |
最終ジャッジ日時 | 2024-07-22 19:35:55 |
合計ジャッジ時間 | 12,757 ms |
ジャッジサーバーID (参考情報) |
judge4 / judge2 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 2 |
other | AC * 12 |
ソースコード
#include <bits/stdc++.h>// #include <boost/multiprecision/cpp_int.hpp>// #define int long long#define inf 1000000007#define pa pair<int,int>#define ll long long#define pal pair<double,double>#define ppap pair<pa,int>#define PI 3.14159265358979323846#define paa pair<int,char>#define mp make_pair#define pb push_back#define EPS (1e-10)int dx[8]={0,1,0,-1,1,1,-1,-1};int dy[8]={1,0,-1,0,-1,1,1,-1};using namespace std;class pa3{public:int x;int y,z;pa3(int x=0,int y=0,int z=0):x(x),y(y),z(z) {}bool operator < (const pa3 &p) const{if(x!=p.x) return x<p.x;if(y!=p.y) return y<p.y;return z<p.z;//return x != p.x ? x<p.x: y<p.y;}bool operator > (const pa3 &p) const{if(x!=p.x) return x>p.x;if(y!=p.y) return y>p.y;return z>p.z;//return x != p.x ? x<p.x: y<p.y;}bool operator == (const pa3 &p) const{return x==p.x && y==p.y && z==p.z;}bool operator != (const pa3 &p) const{return !( x==p.x && y==p.y && z==p.z);}};class pa4{public:int x;int y,z,w;pa4(int x=0,int y=0,int z=0,int w=0):x(x),y(y),z(z),w(w) {}bool operator < (const pa4 &p) const{if(x!=p.x) return x<p.x;if(y!=p.y) return y<p.y;if(z!=p.z)return z<p.z;return w<p.w;//return x != p.x ? x<p.x: y<p.y;}bool operator > (const pa4 &p) const{if(x!=p.x) return x>p.x;if(y!=p.y) return y>p.y;if(z!=p.z)return z>p.z;return w>p.w;//return x != p.x ? x<p.x: y<p.y;}bool operator == (const pa4 &p) const{return x==p.x && y==p.y && z==p.z &&w==p.w;}};class pa2{public:int x,y;pa2(int x=0,int y=0):x(x),y(y) {}pa2 operator + (pa2 p) {return pa2(x+p.x,y+p.y);}pa2 operator - (pa2 p) {return pa2(x-p.x,y-p.y);}bool operator < (const pa2 &p) const{return y != p.y ? y<p.y: x<p.x;}bool operator > (const pa2 &p) const{return x != p.x ? x<p.x: y<p.y;}bool operator == (const pa2 &p) const{return abs(x-p.x)==0 && abs(y-p.y)==0;}bool operator != (const pa2 &p) const{return !(abs(x-p.x)==0 && abs(y-p.y)==0);}};class Point{public:double x,y;Point(double x=0,double y=0):x(x),y(y) {}Point operator + (Point p) {return Point(x+p.x,y+p.y);}Point operator - (Point p) {return Point(x-p.x,y-p.y);}Point operator * (double a) {return Point(x*a,y*a);}Point operator / (double a) {return Point(x/a,y/a);}double absv() {return sqrt(norm());}double norm() {return x*x+y*y;}bool operator < (const Point &p) const{return x != p.x ? x<p.x: y<p.y;}bool operator == (const Point &p) const{return fabs(x-p.x)<EPS && fabs(y-p.y)<EPS;}};typedef Point Vector;#define pl pair<int,pas>struct Segment{Point p1,p2;};double dot(Vector a,Vector b){return a.x*b.x+a.y*b.y;}double cross(Vector a,Vector b){return a.x*b.y-a.y*b.x;}bool parareru(Point a,Point b,Point c,Point d){// if(abs(cross(a-b,d-c))<EPS)cout<<"dd "<<cross(a-b,d-c)<<endl;return abs(cross(a-b,d-c))<EPS;}double distance_ls_p(Point a, Point b, Point c) {if ( dot(b-a, c-a) < EPS ) return (c-a).absv();if ( dot(a-b, c-b) < EPS ) return (c-b).absv();return abs(cross(b-a, c-a)) / (b-a).absv();}bool is_intersected_ls(Segment a,Segment b) {if(a.p1==b.p1||a.p2==b.p1||a.p1==b.p2||a.p2==b.p2) return false;if(parareru((a.p2),(a.p1),(a.p1),(b.p2))&¶reru((a.p2),(a.p1),(a.p1),(b.p1))){// cout<<"sss"<<endl;if(dot(a.p1-b.p1,a.p1-b.p2)<EPS) return true;if(dot(a.p2-b.p1,a.p2-b.p2)<EPS) return true;if(dot(a.p1-b.p1,a.p2-b.p1)<EPS) return true;if(dot(a.p1-b.p2,a.p2-b.p2)<EPS) return true;return false;}else return ( cross(a.p2-a.p1, b.p1-a.p1) * cross(a.p2-a.p1, b.p2-a.p1) < EPS ) && ( cross(b.p2-b.p1, a.p1-b.p1) * cross(b.p2-b.p1,a.p2-b.p1) < EPS );}double segment_dis(Segment a,Segment b){if(is_intersected_ls(a,b))return 0;double r=distance_ls_p(a.p1, a.p2, b.p1);r=min(r,distance_ls_p(a.p1, a.p2, b.p2));r=min(r,distance_ls_p(b.p1, b.p2, a.p2));r=min(r,distance_ls_p(b.p1, b.p2, a.p1));return r;}Point intersection_ls(Segment a, Segment b) {Point ba = b.p2-b.p1;double d1 = abs(cross(ba, a.p1-b.p1));double d2 = abs(cross(ba, a.p2-b.p1));double t = d1 / (d1 + d2);return a.p1 + (a.p2-a.p1) * t;}string itos( int i ) {ostringstream s ;s << i ;return s.str() ;}int gcd(int v,int b){if(v>b) return gcd(b,v);if(v==b) return b;if(b%v==0) return v;return gcd(v,b%v);}double distans(double x1,double y1,double x2,double y2){double rr=(x1-x2)*(x1-x2)+(y1-y2)*(y1-y2);return sqrt(rr);}int mod;ll extgcd(ll a, ll b, ll &x, ll &y) {if (b == 0ll) {x = 1ll;y = 0ll;return a;}ll d = extgcd(b, a%b, y, x);y -= a/b * x;return d;}pa operator+(const pa & l,const pa & r) {return {l.first+r.first,l.second+r.second};}int pr[1200010];int inv[1200010];int beki(int wa,int rr,int warukazu){if(rr==0) return 1%warukazu;if(rr==1) return wa%warukazu;wa%=warukazu;if(rr%2==1) return ((ll)beki(wa,rr-1,warukazu)*(ll)wa)%warukazu;ll zx=beki(wa,rr/2,warukazu);return (zx*zx)%warukazu;}double bekid(double w,int r){if(r==0) return 1.0;if(r==1) return w;if(r%2) return bekid(w,r-1)*w;double f=bekid(w,r/2);return f*f;}int comb(int nn,int rr){int r=pr[nn]*inv[rr];r%=mod;r*=inv[nn-rr];r%=mod;return r;}void gya(int ert){pr[0]=1;for(int i=1;i<=ert;i++){pr[i]=(pr[i-1]*i)%mod;}for(int i=0;i<=ert;i++) inv[i]=beki(pr[i],mod-2,mod);}// cin.tie(0);// ios::sync_with_stdio(false);//priority_queue<pa3,vector<pa3>,greater<pa3>> pq;//sort(ve.begin(),ve.end(),greater<int>());//mt19937(clock_per_sec);// mt19937_64 rng(chrono::steady_clock::now().time_since_epoch().count()) ;//----------------kokomade tenpure------------ll modinv(ll r,ll mo){// mod mo での r の逆数 ないときは-1// r*x +mo*r = 1ll aa,bb;ll gc=extgcd(r,mo,aa,bb);if(gc!=-1){// cout<<-1<<endl;}if(gc!=1) return -1;aa%=mo;if(aa<0)aa+=mo;return aa;}ll garner(vector<pair<ll,ll>> vec,ll MOD){// vec[i]: x % vec[i].first == vec[i].second// 0<= vec[i].second<vec[i].first// vec[i].first はどの2つも互いに素!!!!// MOD はどうでもいいvec.pb({MOD,0ll});int siz=vec.size();vector<ll> ruiseki(siz,1ll);vector<ll> teisuu(siz,0ll);ll val;for(int i=0;i<siz-1;i++){// (teisuu[i] + val*ruiseki[i]) % vec[i].first = vec[i].secondval= (ll)vec[i].second-teisuu[i];while(val<0) val+=(ll)vec[i].first;val%=(ll)vec[i].first;val*= modinv(ruiseki[i],(ll)vec[i].first);val%=(ll)vec[i].first;for(int j=i+1;j<siz;j++){teisuu[j]+=(ll)ruiseki[j]*val%(ll)vec[j].first;teisuu[j]%=(ll)vec[j].first;// cout<<teisuu[j]<<endl;ruiseki[j]*=(ll)vec[i].first;ruiseki[j]%=(ll)vec[j].first;}}return teisuu.back();}int mo[3]={17,9920467,592951213};int a[10000020];signed main(){cin.tie(0);ios::sync_with_stdio(false);int n,m;cin>>n>>m;vector<int> ve;for(int i=0;i<m;i++){int y;cin>>y;ve.pb(y);}vector<pair<ll,ll>> w;for(int k=0;k<3;k++){a[1]=1;int m=(n+1)/2;for(int i=2;i<=m;i++){int e=0;for(auto v:ve)if(i-v>0){e+=a[i-v];e%=mo[k];}a[i]=e%mo[k];}int m2=n-m;ll q=0;for(auto v:ve){for(int i=m;i>max(m-v,0);i--){int pos=m2-(v-(m-i)-1);if(pos>0){q+=(ll)a[i]*a[pos]%mo[k];q%=mo[k];}}}w.pb({(ll)mo[k],(ll)q});}cout<<garner(w,100000000000000007ll)<<endl;return 0;}