麻豆小视频在线观看_中文黄色一级片_久久久成人精品_成片免费观看视频大全_午夜精品久久久久久久99热浪潮_成人一区二区三区四区

首頁 > 學院 > 開發設計 > 正文

Poj 2976 Dropping tests(01分數規劃 牛頓迭代)

2019-11-14 09:49:40
字體:
來源:轉載
供稿:網友

Dropping tests Time Limit: 1000MS Memory Limit: 65536K Description In a certain course, you take n tests. If you get ai out of bi questions correct on test i, your cumulative average is defined to be Given your test scores and a positive integer k, determine how high you can make your cumulative average if you are allowed to drop any k of your test scores. Suppose you take 3 tests with scores of 5/5, 0/1, and 2/6. Without dropping any tests, your cumulative average is . However, if you drop the third test, your cumulative average becomes . Input The input test file will contain multiple test cases, each containing exactly three lines. The first line contains two integers, 1 ≤ n ≤ 1000 and 0 ≤ k < n. The second line contains n integers indicating ai for all i. The third line contains n positive integers indicating bi for all i. It is guaranteed that 0 ≤ ai ≤ bi ≤ 1, 000, 000, 000. The end-of-file is marked by a test case with n = k = 0 and should not be PRocessed. Output For each test case, write a single line with the highest cumulative average possible after dropping k of the given test scores. The average should be rounded to the nearest integer. Sample Input 3 1 5 0 2 5 1 6 4 2 1 2 7 9 5 6 7 9 0 0 Sample Output 83 100 Hint To avoid ambiguities due to rounding errors, the judge tests have been constructed so that all answers are at least 0.001 away from a decision boundary (i.e., you can assume that the average is never 83.4997). Source Stanford Local 2005

/*裸的01分數規劃問題.令∑a[i]/∑b[i]=ans. 則∑a[i]-∑b[i]*ans=0. 二分一個ans.然后用a[i]-b[i]*ans取前k大檢驗.只能去感性的認識orz...并不會證明.. */#include<iostream>#include<cstdio>#include<algorithm>#define eps 1e-7#define MAXN 1001using namespace std;double ans,a[MAXN],b[MAXN],sum,tmp[MAXN];int n,m,k;bool check(double x){ double tot=0; for(int i=1;i<=n;i++) tmp[i]=a[i]-x*b[i]; sort(tmp+1,tmp+n+1,greater<double>()); for(int i=1;i<=n-k;i++) tot+=tmp[i]; if(tot>=0) return true; else return false;}void slove(){ double l=0,r=1e4,mid; while(l<=r) { mid=(l+r)/2.0; if(check(mid)) l=mid+eps,ans=mid; else r=mid-eps; } printf("%.0f/n",ans*100); return ;}int main(){ while(scanf("%d%d",&n,&k)) { if(!n&&!k) break; sum=ans=0; for(int i=1;i<=n;i++) scanf("%lf",&a[i]); for(int i=1;i<=n;i++) scanf("%lf",&b[i]); slove(); } return 0;}/*發現這題牛頓迭代可做吖.網上的題解都是二分01規劃的.我就寫個牛頓迭代的吧orz(雖然二分的寫過).先選一個估計值s0.我們能保證這個答案是單調的.假設上次迭代的ans為s1,則存在n-k個元素使s1=∑(ai/bi),變形可得到∑ai-s2*∑bi=0,令ans[i]=a[i]-b[i]*s0.取前n-k大統計一個答案.可知必存在n-k個元素使∑ansi=∑ai-s1*∑bi=0,所以當我們按ans排序并取前n-k個元素作為求其∑ans時,∑ansi顯然是>=0的,然后s1=(∑ai-∑ansi)/∑bi)<=(∑ai/∑bi)=s2(i<=n-k).即此迭代過程是收斂的,當等號成立時,s即為答案.有些地方還是有點想不通畢竟弱吖orz. */#include<cstdio>#include<algorithm>#include<cmath>#define MAXN 1001#define eps 1e-7using namespace std;double ans,sum,tmp[MAXN];int n,m,k;struct data{double a,b,ans;}s[MAXN];bool cmp(const data &x,const data &y){ return x.ans>y.ans;}void slove(){ double suma=0,sumb=0,s0=0,s1=0; for(int i=1;i<=k;i++) suma+=s[i].a,sumb+=s[i].b; s0=suma/sumb; while(abs(s0-s1)>eps) { s1=s0;suma=sumb=0; for(int i=1;i<=n;i++) s[i].ans=s[i].a-s[i].b*s0; sort(s+1,s+n+1,cmp); for(int i=1;i<=k;i++) suma+=s[i].a,sumb+=s[i].b; s0=suma/sumb; } printf("%.0f/n",s0*100); return ;}int main(){ while(~scanf("%d%d",&n,&k)) { if(!n&&!k) break; sum=ans=0;k=n-k; for(int i=1;i<=n;i++) scanf("%lf",&s[i].a); for(int i=1;i<=n;i++) scanf("%lf",&s[i].b); slove(); } return 0;}
發表評論 共有條評論
用戶名: 密碼:
驗證碼: 匿名發表
主站蜘蛛池模板: 特级黄aaaaaaaaa毛片 | 一级毛片免费在线 | 久色一区| 久草视频国产在线 | 中文字幕www | 最近国产中文字幕 | 欧美一级色片 | 娇喘在线 | 国产91小视频在线观看 | 一级做a爰片性色毛片2021 | 久久久久久久一区二区三区 | 成人福利在线看 | 91短视频版高清在线观看www | 中文在线免费观看 | 亚洲国产精品久久久久久久久久久 | 成人一区二区在线观看视频 | 免费黄色在线 | 午夜av男人的天堂 | 91短视频版高清在线观看www | 国产一级性生活视频 | 欧洲成人一区 | 亚洲欧美国产高清 | 日本成年网 | 欧洲精品久久久 | 成人免费乱码大片a毛片视频网站 | 色骚综合| 极品美女一级毛片 | 日日草天天干 | 久草在线综合 | 毛片免费大全短视频 | 成人免费看视频 | 国产亚洲欧美日韩在线观看不卡 | 国产精品视频久久久 | 成人小视频免费在线观看 | 欧美精品一区自拍a毛片在线视频 | 免费视频xxxx | 国产精品免费小视频 | 高清视频91 | 视频一区二区三区免费观看 | 黄wwww| 看免费5xxaaa毛片|