[SGU174 Walls]并差集、排序+二分查找

【题目大意】
给出一些只在端点相交的线段。
问在连了第几条边之后,这些线段中的一部分连成了一个完整的封闭图形。
【算法分析】
封闭图形就是一个环,如果已经成同一连通分量的两个点再被连,就是成环了。
这个可以用并差集判断。
然后我们需要快速查找点。
由于本题都是静态的东西,所以快排+二分足矣。
如果需要动态算法可以考虑hash表和平衡树。
【其它】1A
【CODE】

#include #include #include struct Point{int x,y;}list[400005],a[200005][2];
int n,tot,f[400005];

void input(){
scanf("%d",&n);
for (int i=1;i<=n;i++){
scanf("%d%d%d%d",&a[i][0].x,&a[i][0].y,&a[i][1].x,&a[i][1].y);
list[++tot]=a[i][0];
list[++tot]=a[i][1];
}
}

int Findst(int x){
int l=1,r=tot,mid;
while (l mid=(l+r)/2;
if (x>list[mid].x) l=mid+1;
else r=mid;
}
return l;
}

int Finded(int x){
int l=1,r=tot,mid;
while (l+1 mid=(l+r)/2;
if (x>=list[mid].x) l=mid;
else r=mid-1;
}
if (x==list[r].x) return r;
return l;
}

int Find(Point tmp){
int l=Findst(tmp.x),r=Finded(tmp.x),mid;
while (l mid=(l+r)/2;
if (tmp.y>list[mid].y) l=mid+1;
else r=mid;
}
return l;
}

int cmp(const void *x,const void *y){
if (((Point*)x)->x==((Point*)y)->x) return ((Point*)x)->y-((Point*)y)->y;
return ((Point*)x)->x-((Point*)y)->x;
}

int GF(int p){
if (f[p]==p) return p;
return f[p]=GF(f[p]);
}

void solve(){
for (int i=1;i<=tot;i++) f[i]=i;
for (int x,y,i=1;i<=n;i++){
x=Find(a[i][0]);
y=Find(a[i][1]);
if (GF(x)==GF(y)){
printf("%dn",i);
return;
}
f[f[x]]=f[y];
}
printf("0n");
}

int main(){
input();
qsort(list+1,tot,sizeof(Point),cmp);
solve();
}

[SGU206 Roads]KM算法的妙用

【题目大意】

给定N个城市和M条道路,其中前N-1条道路是大路,其它都是小路,并保证大 路可构成N个城市 的生成树。

第i条路费用为ci,要求你谎报费用(设第i条路上报的费用为di),使得对费用数组d来说,N-1条大路恰为最小生成树,并且sum|ci-di|尽可能小。

(摘自:http://crfish.blogbus.com/logs/62846336.html

【算法分析】

本题的解题报告来自:http://wenku.baidu.com/view/7bfaa802de80d4d8d15a4faf.html

大概是这样:

令w[i]表示第i条路需要改变的量,且w[i]>=0

假设该路是大路,显然d[i]=c[i]-w[i];

假设该路是小路,显然d[i]=c[i]+w[i];

然后对于树环有限制条件:c[i]-w[i]<=c[j]+w[j]  (其中i为大路,j为小路)

变形成:w[i]+w[j]>=c[i]-c[j]。

于是KM里有不等式w[i]+w[j]>=P[i,j]。

然后由于最后KM的顶标之和等于最佳匹配的总权值,而对于该题而言,这个总权又是必须的。

于是两个问题等价了。

【CODE】

#include #include #include struct gtp{int x,y,w,op,next;}g[810];
int n,m,p,dis;
int w[405][405],fa[66],depth[66],lx[405],ly[405],link[405],ls[66];
bool cx[405],cy[405];

inline void swap(int &x,int &y){x=x^y; y=x^y; x=x^y;}
inline int max(int x,int y){return x>y?x:y;}
inline int min(int x,int y){return x scanf("%d%d",&n,&m);
for (int i=1;i<=m;i++){
scanf("%d%d%d",&g[i].x,&g[i].y,&g[i].w);
g[i].next=ls[g[i].x];
ls[g[i].x]=i;
g[i].op=i+m;
}
for (int i=m+1;i<=2*m;i++){
g[i].x=g[i-m].y;
g[i].y=g[i-m].x;
g[i].w=g[i-m].w;
g[i].next=ls[g[i].x];
ls[g[i].x]=i;
g[i].op=i-m;
}
}

void buildtree(int p){
for (int t=ls[p];t;t=g[t].next)
if ((t fa[g[t].y]=t;
depth[g[t].y]=depth[p]+1;
buildtree(g[t].y);
}
}

void buildgraph(){
memset(w,0,sizeof(w));
int i,x,y,Fa;
for (i=n;i<=m;i++){
x=g[i].x; y=g[i].y;
while (y!=x){
if (depth[x]>depth[y]) swap(x,y);
Fa=fa[y];
if (Fa>=n) Fa=g[Fa].op;
w[Fa][i-n+1]=max(0,g[Fa].w-g[i].w);
y=g[fa[y]].x;
}
}
p=max(n-1,m-n+1);
}

bool find(int i){
cx[i]=true;
int q;
for (int j=1;j<=p;j++)
if (lx[i]+ly[j]==w[i][j] && !cy[j]){
cy[j]=true;
q=link[j];
link[j]=i;
if (!q || find(q)) return true;
link[j]=q;
}
else if (!cy[j]) dis=min(dis,lx[i]+ly[j]-w[i][j]);
return false;
}

void KM(){
int i,j;
for (i=1;i<=p;i++)
for (j=1;j<=p;j++)
if (w[i][j]>lx[i])
lx[i]=w[i][j];
for (i=1;i<=p;i++){
for (;;){
for (j=1;j<=p;j++) cx[j]=false;
for (j=1;j<=p;j++) cy[j]=false;
dis=0x7FFFFFFF;
if (find(i)) break;
for (j=1;j<=p;j++) if (cx[j]) lx[j]-=dis;
for (j=1;j<=p;j++) if (cy[j]) ly[j]+=dis;
}
}
}

void output(){
for (int i=1;i<=p;i++) g[i].w-=lx[i];
for (int i=1;i<=p;i++) g[i+n-1].w+=ly[i];
for (int i=1;i<=m;i++)
printf("%dn",g[i].w);
}

int main(){
input();
buildtree(1);
buildgraph();
KM();
output();
}

[HDOJ3402 Ants run!]排序

【算法分析】
排序,然后判断就可以。
对于每只蚂蚁,我们只考虑他什么时候能追上前面那只蚂蚁就可以,因为被追上与之是对应的。
所以我们要尽量让后面那只蚂蚁难追上前面那只蚂蚁。
所以我们排序即可。
【CODE】
#include #include #include #include #include using namespace std;
const double pi=3.141592653589793238462643383279502884197169399375;
int n;
int a[100005];
double r;

void input(){
scanf("%d %lf",&n,&r);
for (int i=0;i scanf("%d",&a[i]);
r=2*r*pi/n;
}

void solve(){
double ans=1e50;
for (int i=1;i if (a[i]!=a[i-1] && r/(a[i]-a[i-1]) ans=r/(a[i]-a[i-1]);
if (ans==1e50) printf("Infn");
else printf("%.3lfn",ans);
}

int main(){
int Tc;
scanf("%d",&Tc);
for (int i=1;i<=Tc;i++){
input();
sort(a,a+n);
solve();
}
}

[HDOJ3403 Palindrome day]暴力

【算法分析】
直接预处理出所有答案。
枚举回文数的前半截。
这样就可以确保升序。
注意判断闰年。
比赛时因为这个闰年的某个判断没AC。= =
【其它】
我真是沙茶。。。
【CODE】
#include #include #include typedef __int64 lld;
const int N=4000000;
int tot,Tc,len,odd;
lld list[4000005];
int L1[13]={0,31,29,31,30,31,30,31,31,30,31,30,31};
int now,mid,dd,mm,Ed[10];
lld yy;
char ch[100];

inline lld filp(lld x){return x/10+x%10*10;}

void Add(int LL){
if (tot==N) return;
int i,j;
if (len%2==0) j=LL;
else j=LL-1;
while (j>0){
ch[++LL]=ch[j];
j–;
}
yy=dd*100+mm;
for (i=1;i<=LL;i++)
yy=yy*10+ch[i]-‘0’;
if (!((yy%4==0 && yy%100!=0) || yy%400==0))
if (filp(mm)==2 && filp(dd)==29) return;
yy=0;
for (i=1;i<=LL;i++)
yy=yy*10+ch[i]-‘0’;
list[++tot]=yy*10000+mm*100+dd;
}

void Try(){
int LL,i,j;
LL=(len+1)/2;
for (i=1;i<=LL;i++) ch[i]='0';
Add(LL);
if (tot==N) return;
while (1){
i=LL;
while (i>0 && ch[i]==’9′) i–;
if (i==0) return;
for (j=i+1;j<=LL;j++) ch[j]='0';
ch[i]++;
Add(LL);
if (tot==N) return;
}
}

bool can(){
int td=filp(dd),tm=filp(mm);
if (tm<1 || tm>12) return false;
if (td<1 || td>L1[tm]) return false;
return true;
}

void Get(){
for (dd=10;dd<=99;dd++)
for (mm=1;mm<=99;mm++)
if (can())
Try();
}

void init(){
tot=0;
for (len=0;tot Get();
Ed[len]=tot;
}
}

void PUT(lld x,int LL){
if (LL==0) return;
lld tmp=1;
for (int i=1;i while (tmp>0){
printf("%d",(int)(x/tmp%10));
tmp/=10;
}
}

void output(lld x,int LL){
PUT(x%100,2);
PUT(x/100%100,2);
PUT(x/10000,LL);
PUT(filp(x/100%100),2);
PUT(filp(x%100),2);
printf("n");
}

void solve(){
int Tc,x;
scanf("%d",&Tc);
for (int i=1;i<=Tc;i++){
scanf("%d",&x);
int j=0;
while (x>Ed[j]) j++;
output(list[x],j);
}
}

int main(){
init();
solve();
}

HDU月赛——2010.5.1

本次月赛出了6题。Rank 2。
1002果断TLE。然后再用拓展欧几里得算法,再次WA到死。
然后最后rejudge了。
回文数的正解居然是要排序的,我写的是按顺序弄出来的。。。一路狂WA。。。
最后就这样悲剧了。
悲剧+沙茶。。。
没啥说的。。。BS我这沙茶吧
1001:http://hi.baidu.com/edwardmj/blog/item/7973c631706b6af41a4cff36.html
1002:http://hi.baidu.com/edwardmj/blog/item/432eec95b9c7bc6155fb9628.html
1003:http://hi.baidu.com/edwardmj/blog/item/bdd0b333ad4f3215ebc4af59.html
1004:http://hi.baidu.com/edwardmj/blog/item/b40d62d19d4bc038970a1612.html
1005:http://hi.baidu.com/edwardmj/blog/item/82fcbcc49ecc4f179d163d21.html
1006:http://hi.baidu.com/edwardmj/blog/item/46f4fcd30c4a34349b5027e2.html
1007:http://hi.baidu.com/edwardmj/blog/item/ddefd41529d7340d4a90a71c.html
另外1003在SCOI测试中,匹配是比并查集快的,但是这里匹配会TLE,只能用并差集。