麻豆小视频在线观看_中文黄色一级片_久久久成人精品_成片免费观看视频大全_午夜精品久久久久久久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;}
發表評論 共有條評論
用戶名: 密碼:
驗證碼: 匿名發表
主站蜘蛛池模板: 欧美18—19sex性护士中国 | 日本在线不卡一区二区三区 | 国产91av视频 | 成人羞羞视频在线观看免费 | 国产精品美女久久久免费 | 国产精品久久久久久久久久久久久久久久 | 亚洲一区国产二区 | 国产免费久久久久 | 日韩电影一区二区 | 麻豆视频免费网站 | 久久久久久久久久久久久久久久久久 | 久操福利视频 | 国产精品久久久久网站 | 成人短视频在线观看免费 | 久久人人97超碰国产公开结果 | 国产免费久久久久 | 精品国产精品久久 | 日本黄色免费片 | 国产精品久久久久久久久粉嫩 | 久久国产一 | 日本一区二区高清不卡 | 国产精品av久久久久久网址 | 成人黄色小视频在线观看 | 国产亚洲精品yxsp | 国产91九色在线播放 | 成人福利在线免费观看 | 国产精选91 | 午夜视频在线 | 成人在线视频黄色 | 国内精品伊人久久久久网站 | 日韩欧美激情视频 | 久久精品网站视频 | 免费欧美精品 | 成人福利在线视频 | 成人毛片100部免费观看 | 九九热精品视频在线 | 国产午夜精品在线 | 黄色一级片免费在线观看 | 成人毛片在线免费看 | 黄视频网站免费观看 | av电影免费在线看 |