題目大意:給出一個正整數n(n為合數),求n的一個劃分(a1,a2,...,ak,...)(k>=2)。使得其在存在最大的最大公約數之下,存在最大的最小公倍數。
Time Limit:500MS Memory Limit:65536KB 64bit IO Format:%I64d & %I64u
數據規模:200<=n<=10^9。
理論基礎:關於LCM的一個定理:LCM(a1,a2,a3,...,ai)為各個數分解質因數以後,所有存在的素數的最高冪次之積。具體參見鏈接1。
題目分析:
首先,要想有最大公公約數,那麼假設它們的最大公約數為G,那麼:G*(a1/G,a2/G,a3,...)=n,所以要想讓G最大,那麼只需要a1/G+a2/G+a3/G+...最小記為M,那麼取M為n的最小素數因子即可。這樣一來:M<=sqrt(10^9)<31625。問題轉換為求素數M的一個劃分,使劃分中的數的嘴小公倍數最大。
根據定理,我們可以斷言:將拆分成M=若干個1+p1^k1+p2^k2+p3^k3+...(其中p1,p2,p3,...等數兩兩互質且為素數)時,最小公倍數最大。
證明:首先,我們證明兩兩互質時最大。假設LCM最大時,仍有a1<=a2存在最大公約數k,那麼a1=k*p,a2=k*q,這樣一來a1與a2的最小公倍數為:k*p*q。那麼我們進行如下操作。因為:a1,a2不互質,所以k>=2,那麼我們將:a2拆成:a2=q+(k-1)*q,這樣q,q(k-1),k*p,三個數兩兩互質且最小公倍數大於原來的最小公倍數,所以必然兩兩互質。
其次,我們證明,每個數都可以表示為:若干個1+p^k(p為素數),假設:存在合數(k*q)^i(k,p均為質數,且>=2),那麼此數可以證明:k^i*p^i>k^i+p^i。(利用二元函數f(x,y)=(xy)^i-x^i-y^i可以證明,這裡不多展開了)。所以,最後變為對各個素數的指數的動態規劃問題。那麼指標該怎麼定呢?因為每兩個數都是互質的,所以其最小公倍數即為已經存在的各個數之積,因為不需要求解LCM,所以大可用double存儲。如果害怕double過大時,存在前十六位相同,導致狀態轉移錯無,那我們可以用log來表示,因為:log是單增函數,且log(a*b)=loga+logb,便於狀態轉移。
double dp[N][M];
int pre[i][j];
我們用dp[i][j]表示用前i個數表示j時最優解,即log的最大值。
用pre[i][j]表示取dp[i][j]為最優解表示j時,j剩余的數。
狀態轉移方程:dp[i][j]=max(dp[i-1][j-p[i]]+logp[i],dp[i][j]);這樣我們最後的最優解即為:dp[N][j]裡最大的那一個然後:如果:j-pre[N][j]不為零的話就要輸出(j-pre[N][j])*G,下來搜索dp[N-11][j-prime[N]]。最後dp[1][j]剩下來的值即為要輸出的1*G的個數了。
代碼如下:
#include<iostream> #include<cstring> #include<string> #include<cstdlib> #include<cstdio> #include<cmath> #include<algorithm> #include<queue> #include<ctime> #include<vector> using namespace std; typedef double db; #define DBG 0 #define maa (1<<31) #define mii ((1<<31)-1) #define sl(c) ((int)(c).size()) //取字符串長度; #define forl(i, a, b) for(int i = (a); i < (b); ++i) //不帶邊界值循環,正序 #define forle(i, a, b) for(int i = (a); i <= (b); ++i) //帶邊界值循環,正序 #define forh(i, a, b) for(int i = (a); i > (b); --i) //不帶邊界值,逆序 #define forhe(i, a, b) for(int i = (a); i >= (b); --i) //帶邊界值,逆序 #define forlc(i, a, b) for(int i##_b = (b), i = (a); i < i##_b; ++i) //帶別名的循環,用於不可改變值 #define forlec(i, a, b) for(int i##_b = (b), i = (a); i <= i##_b; ++i) #define forgc(i, a, b) for(int i##_b = (b), i = (a); i > i##_b; --i) #define forgec(i, a, b) for(int i##_b = (b), i = (a); i >= i##_b; --i) #define forall(i, v ) forl(i, 0, sz(v)) //循環所有 #define forallc(i, v ) forlc(i, 0, sz(v)) #define forlla(i, v ) forhe(i, sz(v)-1, 0) #define forls(i, n, a, b) for(int i = a; i != b; i = n[i]) //搜表用 #define rep(n) for(int repp = 0; repp < (n); ++repp) #define repc(n) for(int repp_b = (n), repp = 0; repp < repp_b; ++repp) #define rst(a, v) memset(a, v, sizeof a) //把字符v填充到a reset 重置 #define cpy(a, b) memcpy(a, b, sizeof a) //copy b 的sizeof(a)個字符到a #define rstn(a, v, n) memset(a, v, (n)*sizeof((a)[0])) //把字符v填充到a[n]之前的字節 #define cpyn(a, b, n) memcpy(a, b, (n)*sizeof((a)[0])) //copy b 的 n 個字符到a #define ast(b) if(DBG && !(b)) { printf("%d!!|\n", __LINE__); while(1) getchar(); } //調試 #define dout DBG && cout << __LINE__ << ">>| " #define pr(x) #x"=" << (x) << " | " #define mk(x) DBG && cout << __LINE__ << "**| "#x << endl #define pra(arr, a, b) if(DBG) {\ dout<<#arr"[] |" <<endl; \ forlec(i, a, b) cout<<"["<<i<<"]="<<arr[i]<<" |"<<((i-(a)+1)%8?" ":"\n"); \ if((b-a+1)%8) puts("");\ } //數列查看 #define rd(type, x) type x; cin >> x //讀數 inline int rdi() { int d; scanf("%d", &d); return d; } inline char rdc() { scanf(" "); return getchar(); } inline string rds() { rd(string, s); return s; } inline db rddb() { db d; scanf("%lf", &d); return d; } template<class T> inline bool updateMin(T& a, T b) { return a>b? a=b, true: false; } template<class T> inline bool updateMax(T& a, T b) { return a<b? a=b, true: false; } typedef long long LL; typedef long unsigned LU; const int N=130; const int M=40000; double dp[N][M]; bool nprime[M]={true,true}; vector<int> prim; int pre[N][M],n,m; void init() { int bor=sqrt((float)M),temp; forl(i,2,bor) { if(!nprime[i]) { prim.push_back(i); for(int j=i;(temp=j*i)<=M;j++) nprime[temp]=true; } } forlec(j,bor+1,M) { if(!nprime[j])prim.push_back(j); } pra(prim,0,prim.size()-1); } void solve(int n) { rst(dp,0); int gcd=0,div; for(int i=2;prim[i]*prim[i]<=n;i++) { if(n%prim[i]) continue; gcd=n/prim[i]; break; } div=n/gcd; forl(i,1,N) { int k=prim[i-1]; forlec(j,1,div) { dp[i][j]=dp[i-1][j]; pre[i][j]=j; } while(k<=div) { forlec(j,k,div) { if(dp[i][j]<dp[i-1][j-k]+log((double)k)) { dp[i][j]=dp[i-1][j-k]+log((double)k); pre[i][j]=j-k; } } k*=prim[i-1]; } } int res[500],cnt=0,state=div,now=N-1; int lcm=0; forle(j,0,div) { if(lcm<dp[now][j]) { lcm=dp[now][j]; state=j; } } while(now) { int temp=pre[now][state]; if((state-temp)!=0)res[cnt++]=state-temp; state=temp; now--; } rep(state)res[cnt++]=1; for(int i=0;i<cnt;i++) printf("%d%c",res[i]*gcd,i==cnt-1?'\n':' '); } int main() { init(); while(scanf("%d",&n)!=EOF) { if(n%2==0) { printf("%d %d\n",n/2,n/2); continue; } if(n%3==0) { printf("%d %d\n",n/3*2,n/3); continue; } solve(n); } return 0; } #include<iostream> #include<cstring> #include<string> #include<cstdlib> #include<cstdio> #include<cmath> #include<algorithm> #include<queue> #include<ctime> #include<vector> using namespace std; typedef double db; #define DBG 0 #define maa (1<<31) #define mii ((1<<31)-1) #define sl(c) ((int)(c).size()) //取字符串長度; #define forl(i, a, b) for(int i = (a); i < (b); ++i) //不帶邊界值循環,正序 #define forle(i, a, b) for(int i = (a); i <= (b); ++i) //帶邊界值循環,正序 #define forh(i, a, b) for(int i = (a); i > (b); --i) //不帶邊界值,逆序 #define forhe(i, a, b) for(int i = (a); i >= (b); --i) //帶邊界值,逆序 #define forlc(i, a, b) for(int i##_b = (b), i = (a); i < i##_b; ++i) //帶別名的循環,用於不可改變值 #define forlec(i, a, b) for(int i##_b = (b), i = (a); i <= i##_b; ++i) #define forgc(i, a, b) for(int i##_b = (b), i = (a); i > i##_b; --i) #define forgec(i, a, b) for(int i##_b = (b), i = (a); i >= i##_b; --i) #define forall(i, v ) forl(i, 0, sz(v)) //循環所有 #define forallc(i, v ) forlc(i, 0, sz(v)) #define forlla(i, v ) forhe(i, sz(v)-1, 0) #define forls(i, n, a, b) for(int i = a; i != b; i = n[i]) //搜表用 #define rep(n) for(int repp = 0; repp < (n); ++repp) #define repc(n) for(int repp_b = (n), repp = 0; repp < repp_b; ++repp) #define rst(a, v) memset(a, v, sizeof a) //把字符v填充到a reset 重置 #define cpy(a, b) memcpy(a, b, sizeof a) //copy b 的sizeof(a)個字符到a #define rstn(a, v, n) memset(a, v, (n)*sizeof((a)[0])) //把字符v填充到a[n]之前的字節 #define cpyn(a, b, n) memcpy(a, b, (n)*sizeof((a)[0])) //copy b 的 n 個字符到a #define ast(b) if(DBG && !(b)) { printf("%d!!|\n", __LINE__); while(1) getchar(); } //調試 #define dout DBG && cout << __LINE__ << ">>| " #define pr(x) #x"=" << (x) << " | " #define mk(x) DBG && cout << __LINE__ << "**| "#x << endl #define pra(arr, a, b) if(DBG) {\ dout<<#arr"[] |" <<endl; \ forlec(i, a, b) cout<<"["<<i<<"]="<<arr[i]<<" |"<<((i-(a)+1)%8?" ":"\n"); \ if((b-a+1)%8) puts("");\ } //數列查看 #define rd(type, x) type x; cin >> x //讀數 inline int rdi() { int d; scanf("%d", &d); return d; } inline char rdc() { scanf(" "); return getchar(); } inline string rds() { rd(string, s); return s; } inline db rddb() { db d; scanf("%lf", &d); return d; } template<class T> inline bool updateMin(T& a, T b) { return a>b? a=b, true: false; } template<class T> inline bool updateMax(T& a, T b) { return a<b? a=b, true: false; } typedef long long LL; typedef long unsigned LU; const int N=130; const int M=40000; double dp[N][M]; bool nprime[M]={true,true}; vector<int> prim; int pre[N][M],n,m; void init() { int bor=sqrt((float)M),temp; forl(i,2,bor) { if(!nprime[i]) { prim.push_back(i); for(int j=i;(temp=j*i)<=M;j++) nprime[temp]=true; } } forlec(j,bor+1,M) { if(!nprime[j])prim.push_back(j); } pra(prim,0,prim.size()-1); } void solve(int n) { rst(dp,0); int gcd=0,div; for(int i=2;prim[i]*prim[i]<=n;i++) { if(n%prim[i]) continue; gcd=n/prim[i]; break; } div=n/gcd; forl(i,1,N) { int k=prim[i-1]; forlec(j,1,div) { dp[i][j]=dp[i-1][j]; pre[i][j]=j; } while(k<=div) { forlec(j,k,div) { if(dp[i][j]<dp[i-1][j-k]+log((double)k)) { dp[i][j]=dp[i-1][j-k]+log((double)k); pre[i][j]=j-k; } } k*=prim[i-1]; } } int res[500],cnt=0,state=div,now=N-1; int lcm=0; forle(j,0,div) { if(lcm<dp[now][j]) { lcm=dp[now][j]; state=j; } } while(now) { int temp=pre[now][state]; if((state-temp)!=0)res[cnt++]=state-temp; state=temp; now--; } rep(state)res[cnt++]=1; for(int i=0;i<cnt;i++) printf("%d%c",res[i]*gcd,i==cnt-1?'\n':' '); } int main() { init(); while(scanf("%d",&n)!=EOF) { if(n%2==0) { printf("%d %d\n",n/2,n/2); continue; } if(n%3==0) { printf("%d %d\n",n/3*2,n/3); continue; } solve(n); } return 0; }