結果
問題 | No.3015 アンチローリングハッシュ |
ユーザー | Lay_ec |
提出日時 | 2015-11-04 22:20:40 |
言語 | C++11 (gcc 11.4.0) |
結果 |
WA
|
実行時間 | - |
コード長 | 3,909 bytes |
コンパイル時間 | 1,224 ms |
コンパイル使用メモリ | 107,744 KB |
実行使用メモリ | 6,948 KB |
最終ジャッジ日時 | 2024-09-13 12:27:24 |
合計ジャッジ時間 | 4,822 ms |
ジャッジサーバーID (参考情報) |
judge5 / judge2 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | WA | - |
testcase_01 | WA | - |
testcase_02 | WA | - |
testcase_03 | WA | - |
testcase_04 | WA | - |
testcase_05 | WA | - |
testcase_06 | WA | - |
testcase_07 | WA | - |
testcase_08 | WA | - |
testcase_09 | WA | - |
testcase_10 | WA | - |
testcase_11 | WA | - |
testcase_12 | WA | - |
testcase_13 | WA | - |
testcase_14 | WA | - |
testcase_15 | WA | - |
testcase_16 | WA | - |
testcase_17 | WA | - |
testcase_18 | WA | - |
testcase_19 | WA | - |
testcase_20 | WA | - |
testcase_21 | WA | - |
testcase_22 | WA | - |
ソースコード
#include <iostream> #include <string> #include <vector> #include <cmath> #include <algorithm> #include <cstdlib> #include <ctime> #include <cstdio> #include <functional> #include <set> #include <sstream> #include <map> #include <queue> #include <stack> using namespace std; typedef long long ll; typedef long double ld; template <typename T1,typename T2> ld IP(vector<vector<T1> > &B1, int x, vector<vector<T2> > &B2, int y){ ld res=0; for(int i=0;i<B1.size();i++) res+=(ld)B1[i][x]*(ld)B2[i][y]; return res; } //https://en.wikipedia.org/wiki/Gram-Schmidt_process vector<vector<ld> > GS(vector<vector<ll> > &Basis){ //vector<vector<ld> > Basis_=Basis; const int n=Basis.size(); vector<vector<ld> > Basis_(n,vector<ld>(n)); for(int i=0;i<n;i++){ for(int j=0;j<n;j++) Basis_[i][j]=Basis[i][j]; } for(int i=1;i<n;i++){ for(int j=0;j<i;j++){ ld ip=IP(Basis,i,Basis_,j); ld ip2=IP(Basis_,j,Basis_,j); for(int k=0;k<n;k++) Basis_[k][i]-=ip/ip2*Basis_[k][j]; } } return Basis_; } ld MU(vector<vector<ll> > &Basis,int x,vector<vector<ld> > &Basis_,int y){ return IP(Basis,x,Basis_,y)/IP(Basis_,y,Basis_,y); } ld PI2(vector<vector<ll> > &Basis,vector<vector<ld> > &Basis_,int i,int x){ const int n=Basis.size(); vector<ld> pi(n,0); for(int j=i;j<n;j++){ ld ip1=IP(Basis,x,Basis_,j); ld ip2=IP(Basis_,j,Basis_,j); for(int k=0;k<n;k++){ pi[k]+=ip1/ip2*Basis_[k][j]; } } ld res=0; for(int j=0;j<n;j++) res+=pi[j]*pi[j]; return res; } vector<ll> LLL(vector<vector<ll> > Basis){ const int n=Basis.size(); const ld delta=(1.0/4.0)+powl(3.0/4.0,(ld)n/(ld)(n-1)); vector<ll> svp(n); vector<vector<ld> > Basis_=GS(Basis); while(true){ /* for(int i=0;i<n;i++){ for(int j=0;j<n;j++) cout<<Basis[i][j]<<" "; cout<<endl; } cout<<endl; */ for(int i=0;i<n;i++){ for(int j=i-1;j>=0;j--){ ll c_ij=floor(IP(Basis,i,Basis_,j)/IP(Basis_,j,Basis_,j)+0.5); for(int k=0;k<n;k++) Basis[k][i]-=c_ij*Basis[k][j]; Basis_=GS(Basis); } } bool cont=false; for(int i=0;i<n-1;i++){ ld p1=PI2(Basis,Basis_,i,i); ld p2=PI2(Basis,Basis_,i,i+1); if(delta*p1>p2){ for(int j=0;j<n;j++) swap(Basis[j][i],Basis[j][i+1]); Basis_=GS(Basis); cont=true; break; } } if(cont) continue; else{ /* cout<<"result"<<endl; for(int i=0;i<n;i++){ for(int j=0;j<n;j++){ cout<<Basis[i][j]<<" "; } cout<<endl; } */ for(int i=0;i<n;i++) svp[i]=Basis[i][0]; return svp; } } } unsigned long long get_hash(string s, unsigned long long a, unsigned long long b){ unsigned long long hash = 0; for(int i = 0; i < s.size(); i++){ hash = (hash * a + s[i]) % b; } return hash; } int main(){ /* vector<vector<ll> > mat(4,vector<ll>(4)); mat[0][0]=105; mat[1][0]=821; mat[2][0]=404; mat[3][0]=328; mat[0][1]=881; mat[1][1]=667; mat[2][1]=644; mat[3][1]=927; mat[0][2]=181; mat[1][2]=483; mat[2][2]=87; mat[3][2]=500; mat[0][3]=893; mat[1][3]=834; mat[2][3]=732; mat[3][3]=441; LLL(mat); */ ll p,b; cin>>p>>b; swap(p,b); //No.3015 vector<ll> coe; int n; for(n=4;;n+=4){ cout<<n<<endl; vector<vector<ll> > Basis(n,vector<ll>(n,0)); for(int i=0;i<n-1;i++){ Basis[i][i]=-b; Basis[i+1][i]=1; } Basis[0][n-1]=p; coe=LLL(Basis); reverse(coe.begin(),coe.end()); bool ok=true; for(int i=0;i<n;i++) ok&=(abs(coe[i])<=25); if(ok){ /* for(int i=0;i<n;i++) cout<<coe[i]<<" "; cout<<endl; */ break; } } map<int,pair<char,char> > st; for(int i=-25;i<=25;i++){ bool ok=false; for(char s='a';s<='z' && !ok;s++){ for(char t='a';t<='z' && !ok;t++){ if(s-t==i){ st[i]=make_pair(s,t); ok=true; } } } } string s,t; for(int i=0;i<n;i++){ s+=st[coe[i]].first; t+=st[coe[i]].second; } cout<<s+"\n"<<t<<endl; //cout<<get_hash(s,b,p)<<endl; //cout<<get_hash(t,b,p)<<endl; return 0; }