多多色-多人伦交性欧美在线观看-多人伦精品一区二区三区视频-多色视频-免费黄色视屏网站-免费黄色在线

國內最全IT社區平臺 聯系我們 | 收藏本站
阿里云優惠2
您當前位置:首頁 > php開源 > 綜合技術 > 2016 ACM/ICPC Asia Regional Shenyang Online 1009 QSC and Master 區間dp

2016 ACM/ICPC Asia Regional Shenyang Online 1009 QSC and Master 區間dp

來源:程序員人生   發布時間:2016-12-03 10:22:56 閱讀次數:7019次

QSC and Master

Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 131072/131072 K (Java/Others)
Total Submission(s): 773    Accepted Submission(s): 291


Problem Description
Every school has some legends, Northeastern University is the same.

Enter from the north gate of Northeastern University,You are facing the main building of Northeastern University.Ninety-nine percent of the students have not been there,It is said that there is a monster in it.

QSCI am a curious NEU_ACMer,This is the story he told us.

It’s a certain period,QSCI am in a dark night, secretly sneaked into the East Building,hope to see the master.After a serious search,He finally saw the little master in a dark corner. The master said:

“You and I, we're interfacing.please solve my little puzzle!

There are N pairs of numbers,Each pair consists of a key and a value,Now you need to move out some of the pairs to get the score.You can move out two continuous pairs,if and only if their keys are non coprime(their gcd is not one).The final score you get is the sum of all pair’s value which be moved out. May I ask how many points you can get the most?

The answer you give is directly related to your final exam results~The young man~”

QSC is very sad when he told the story,He failed his linear algebra that year because he didn't work out the puzzle.

Could you solve this puzzle?

(Data range:1<=N<=300
1<=Ai.key<=1,000,000,000
0<Ai.value<=1,000,000,000)
 

Input
First line contains a integer T,means there are T(1≤T≤10) test case。

  Each test case start with one integer N . Next line contains N integers,means Ai.key.Next line contains N integers,means Ai.value.
 

Output
For each test case,output the max score you could get in a line.
 

Sample Input
3 3 1 2 3 1 1 1 3 1 2 4 1 1 1 4 1 3 4 3 1 1 1 1
 

Sample Output
0 2 0
 

Source
2016 ACM/ICPC Asia Regional Shenyang Online
 


My Solution

     狀態定義: dp[i][j][u]  u == 1 時表示 當端點 i, j 進行合并時(取出 val[i] 、 val[j] 時) 

                                                          或 i < k < k+1 < j, key[i] 和 key[k] 可以合并 且key[k+1] 和 key[j]合并的, 區間 [i, j]內的最大價值;

                                     u == 0 時表示 當端點 i, j 不進行合并時(不取出 val[i] 、 val[j] 時) 的 區間 [i, j]內的最大價值;


     初始化:全部初始化為 0;

     邊界: 當 j - i == 1 時 如果可以合并 則 dp[i][j][1] = val[i] + val[j];

     狀態轉移方程:如果 key[i], key[j] 可以組合, 則 如果 dp[i + 1][j - 1][1] > 0 則 dp[i][j][1] = max(dp[i][j][1], dp[i+1][j⑴][1] + val[i] + val[j]);

                                                           然后對 區間劃分 k = i ~ j 

                                                                    dp[i][j][0] = max(dp[i][j][0], max(dp[i][k][0], dp[i][k][1]) + max(dp[k+1][j][1] , dp[k+1][j][0]));
                                                                    if(dp[i][k][1] > 0 && dp[k+1][j][1] > 0) dp[i][j][1] = max(dp[i][j][1], dp[i][k][1] + dp[k+1][j][1]);


                                                                   接著是對u == 0時的轉移(即keyi keyj 能合并時可以選擇不合并) 

                                                                    dp[i][j][0] = max(dp[i][j][0], max(dp[i+1][j⑴][1], dp[i+1][j⑴][0]));


                            如果 key[i], key[j] 不能合并, 則
                                                                        if(dp[i+1][j⑴][1] > 0){
                                                                                   dp[i][j][0] = max(dp[i][j][0], max(dp[i+1][j⑴][1], dp[i+1][j⑴][0]));
                                                                                   for(k = i; k < j; k++){
                                                                                             dp[i][j][0] = max(dp[i][j][0], max(dp[i][k][0], dp[i][k][1]) + max(dp[k+1][j][1] , dp[k+1][j][0]));
                                                                                             if(dp[i][k][1] > 0 && dp[k+1][j][1] > 0) dp[i][j][1] = max(dp[i][j][1], dp[i][k][1] + dp[k+1][j][1]);
                                                                                   }


                                                                       }
                                                                      else{
                                                                                  dp[i][j][0] = max(dp[i][j][0], dp[i+1][j⑴][0]);
                                                                                  for(k = i; k < j; k++){
                                                                                            dp[i][j][0] = max(dp[i][j][0], max(dp[i][k][0], dp[i][k][1]) + max(dp[k+1][j][1] , dp[k+1][j][0]));
                                                                                            if(dp[i][k][1] > 0 && dp[k+1][j][1] > 0) dp[i][j][1] = max(dp[i][j][1], dp[i][k][1] + dp[k+1][j][1]);
                                                                                  }
                                                                     }


     復雜度 O(n^3)


#include <iostream> #include <cstdio> #include <cmath> #include <cstring> using namespace std; typedef long long LL; const int maxn = 3e2 + 8; LL key[maxn], val[maxn], dp[maxn][maxn][2]; inline LL gcd(LL a, LL b) { return b == 0 ? a : gcd(b, a % b); } inline bool check(LL a, LL b) { return gcd(a, b) != 1; } int main() { #ifdef LOCAL freopen("1009.txt", "r", stdin); //freopen("1009.out", "w", stdout); #endif // LOCAL // ios::sync_with_stdio(false); cin.tie(0); LL T, n; cin >> T; while(T--){ memset(dp, 0, sizeof dp); cin >> n; for(int i = 0; i < n; i++){ cin >> key[i]; } for(int i = 0; i < n; i++){ cin >> val[i]; } LL i, j, k, ans = 0, x, y; for(i = n - 1; i >= 0; i--){ for(j = i + 1; j < n; j++){ if(check(key[i], key[j])){ if(j - i > 2){ if(dp[i+1][j⑴][1] > 0){ dp[i][j][1] = max(dp[i][j][1], dp[i+1][j⑴][1] + val[i] + val[j]); } for(k = i; k < j; k++){ dp[i][j][0] = max(dp[i][j][0], max(dp[i][k][0], dp[i][k][1]) + max(dp[k+1][j][1] , dp[k+1][j][0])); if(dp[i][k][1] > 0 && dp[k+1][j][1] > 0) dp[i][j][1] = max(dp[i][j][1], dp[i][k][1] + dp[k+1][j][1]); } dp[i][j][0] = max(dp[i][j][0], max(dp[i+1][j⑴][1], dp[i+1][j⑴][0])); } else if(j - i == 1) dp[i][j][1] = val[i] + val[j]; } if(j - i > 2){ if(dp[i+1][j⑴][1] > 0){ dp[i][j][0] = max(dp[i][j][0], max(dp[i+1][j⑴][1], dp[i+1][j⑴][0])); for(k = i; k < j; k++){ dp[i][j][0] = max(dp[i][j][0], max(dp[i][k][0], dp[i][k][1]) + max(dp[k+1][j][1] , dp[k+1][j][0])); if(dp[i][k][1] > 0 && dp[k+1][j][1] > 0) dp[i][j][1] = max(dp[i][j][1], dp[i][k][1] + dp[k+1][j][1]); } } else{ dp[i][j][0] = max(dp[i][j][0], dp[i+1][j⑴][0]); for(k = i; k < j; k++){ dp[i][j][0] = max(dp[i][j][0], max(dp[i][k][0], dp[i][k][1]) + max(dp[k+1][j][1] , dp[k+1][j][0])); if(dp[i][k][1] > 0 && dp[k+1][j][1] > 0) dp[i][j][1] = max(dp[i][j][1], dp[i][k][1] + dp[k+1][j][1]); } } } else if(j - i == 1) dp[i][j][0] = 0; //這行可以疏忽 ^_^ //*/ ans = max(ans, max(dp[i][j][0], dp[i][j][1])); } } //cout << dp[0][n⑴][0] << " " << dp[0][n⑴][1]<< endl; if(ans > 0) cout << ans << "\n"; else cout << 0 << "\n"; } return 0; }

  Thank you!

                                                                                                                                               ------from ProLights

生活不易,碼農辛苦
如果您覺得本網站對您的學習有所幫助,可以手機掃描二維碼進行捐贈
程序員人生
------分隔線----------------------------
分享到:
------分隔線----------------------------
關閉
程序員人生
主站蜘蛛池模板: 欧美偷窥自拍 | 日韩免费福利视频 | 日韩精品在线一区 | 爱爱永久免费视频网站 | 精品久久中文网址 | yellow网站在线观看 | 国产精品久久久久久久久免费 | 性欧美video另类hd亚洲人 | 欧美午夜在线观看理论片 | 欧美日韩激情一区二区三区 | 毛片免费观看视频 | 青青草原亚洲视频 | 91极品国产 | 久久久久亚洲精品影视 | 日本www黄| 中文字幕人成乱码在线观看 | 欧美一区二区三区不卡免费观看 | 久久免费久久 | www一区 | 国产在线91区精品 | 黄色avav| 波多野结衣资源在线观看 | 亚洲图片自拍 | 亚洲国产情侣一区二区三区 | 欧美色老汉| 国产又黄又爽又色的免费 | 亚洲综合精品一区 | 成人国产精品一级毛片视频 | 亚洲精品老司机在线观看 | 精品视频一区二区三三区四区 | 国产片欧美片亚洲片久久综合 | 一级片亚洲| 手机看片日韩在线 | 中文字幕一区二区三区久久网站 | 伊人久久大香线焦综合四虎 | 亚洲国产精品久久久久久网站 | 日韩亚洲欧美一区二区三区 | 最近最新中文字幕免费大全3 | 视频免费在线观看 | 日韩精品一区二区三区乱码 | 精品久久综合一区二区 |