結果

問題 No.695 square1001 and Permutation 4
コンテスト
ユーザー WA_TLE
提出日時 2018-06-08 23:36:06
言語 C++17
(gcc 13.3.0 + boost 1.87.0)
結果
WA  
実行時間 -
コード長 2,794 bytes
コンパイル時間 4,130 ms
コンパイル使用メモリ 127,568 KB
最終ジャッジ日時 2025-01-05 11:00:07
ジャッジサーバーID
（参考情報） 		judge4 / judge2
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 2
other AC * 10 WA * 2
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

#include<deque>
#include<queue>
#include<vector>
#include<algorithm>
#include<iostream>
#include<set>
#include<cmath>
#include<tuple>
#include<string>
#include<chrono>
#include<functional>
#include<iterator>
#include<random>
#include<unordered_set>
#include<array>
#include<map>
#include<iomanip>
#include<assert.h>
using namespace std;
typedef long long int llint;
typedef long double lldo;
#define mp make_pair
#define mt make_tuple
#define pub push_back
#define puf push_front
#define pob pop_back
#define pof pop_front
#define fir first
#define sec second
#define res resize
#define ins insert
#define era erase
/*
cout<<setprecision(20)
cin.tie(0);
ios::sync_with_stdio(false);
*/
const llint mod=100000000000000007LL;
const llint big=2.19e15+1;
const long double pai=3.141592653589793238462643383279502884197;
const long double eps=1e-15;
template <class T,class U>void mineq(T& a,U b){if(a>b){a=b;}}
template <class T,class U>void maxeq(T& a,U b){if(a<b){a=b;}}
llint gcd(llint a,llint b){if(a%b==0){return b;}else return gcd(b,a%b);}
llint lcm(llint a,llint b){return a/gcd(a,b)*b;}
template<class T> void SO(T& ve){sort(ve.begin(),ve.end());}
template<class T> void REV(T& ve){reverse(ve.begin(),ve.end());}
llint LBI(vector<int>&ar,int in){return lower_bound(ar.begin(),ar.end(),in)-ar.begin();}
llint UBI(vector<int>&ar,int in){return upper_bound(ar.begin(),ar.end(),in)-ar.begin();}
llint extgcd(llint a, llint b, llint& x, llint& y) {
    llint d = a;
    if (b != 0) {
        d = extgcd(b, a % b, y, x);
        y -= (a / b) * x;
    } else {
        x = 1; y = 0;
    }
    return d;
}

int main(void){
	cin.tie(0);
	ios::sync_with_stdio(false);
	llint i,j,n,m;cin>>n>>m;
	llint x[5];
	for(i=0;i<m;i++){cin>>x[i];}
	llint modA=168647939LL;
	llint modB=592951213LL;
	//cout<<((llint)modA)*modB<<endl;
	static int ans[10000005]={0};
	ans[0]=1;
	llint ansA=0,ansB=0;
	for(i=0;i<n/2;i++){
		for(j=0;j<m;j++){
			if(i+x[j]<=n/2){ans[i+x[j]]+=ans[i];ans[i+x[j]]%=modA;}
		}
	}
	for(i=0;i<n/2;i++){
		for(j=0;j<m;j++){
			int ter=n-1-(i+x[j]);
			if(ter<0||(n-1)/2<ter){continue;}
			ansA+=((llint)ans[i])*ans[ter];
			ansA%=modA;
		}
	}
	for(i=0;i<=n/2;i++){ans[i]=0;}
	ans[0]=1;
	for(i=0;i<n/2;i++){
		for(j=0;j<m;j++){
			if(i+x[j]<=n/2){ans[i+x[j]]+=ans[i];ans[i+x[j]]%=modB;}
		}
	}
	for(i=0;i<n/2;i++){
		for(j=0;j<m;j++){
			int ter=n-1-(i+x[j]);
			if(ter<0||(n-1)/2<ter){continue;}
			ansB+=((llint)ans[i])*ans[ter];
			ansB%=modB;
		}
	}
	//cout<<ansA<<" "<<ansB<<endl;
	//ansA,ansB
	llint y,z;
	extgcd(modA,modB,y,z);
	llint rrr=0;
	llint itA=ansA*z;itA%=mod;
	llint itB=ansB*y;itB%=mod;
	for(int h=0;h<35;h++){
		if(((llint)modB)&(1LL<<h)){rrr+=itA;}
		if(((llint)modA)&(1LL<<h)){rrr+=itB;}
		rrr%=mod;itA+=itA;itA%=mod;itB+=itB;itB%=mod;
	}
	cout<<rrr<<endl;
	return 0;
}
0