标签：style blog http color os io ar for art
POJ 3253 Fence Repair STL堆操作
我想说，STL里堆是我目前见到最蛋疼的操作。
#include <cstdio> #include <cstring> #include <iostream> #include <algorithm> #include <cstdlib> #include <cmath> #include <utility> #include <vector> #include <queue> #include <map> #include <set> #define max(x,y) ((x)>(y)?(x):(y)) #define min(x,y) ((x)>(y)?(y):(x)) using namespace std; int n,a[20005]; int main() { scanf("%d",&n); for(int i=0; i<n; i++) scanf("%d",&a[i]); vector< int > v(a,a+n); make_heap(v.begin(),v.end(),greater< int >()); long long ans=0; for(int i=1; i<=n-1; i++) { long long tmp=v.front(); pop_heap(v.begin(),v.end(),greater< int >()); v.pop_back(); tmp+=v.front(); pop_heap(v.begin(), v.end(),greater< int >()); v.pop_back(); v.push_back(tmp); push_heap(v.begin(),v.end(),greater< int >()); ans+=tmp; printf("%d\n", tmp); } printf("%lld\n",ans); return 0; }
POJ 3254 Corn Fields 状态压缩递推
uc[]表示相邻两位不同时为1的情况，377种，因此n*377^2即可
注意，位运算优先级是低于==的！
#include <cstdio> #include <cstring> #define MOD 100000000 int n,m,a[20],k; int uc[500],all; int f[20][5000]; int main() { scanf("%d%d",&n,&m); for(int i=0; i<n; i++) for(int j=0; j<m; j++) { scanf("%d",&k); a[i]=(a[i]<<1)+k; } for(int i=0; i<(1<<m); i++) { int tmi=i,flag=0; while(tmi) { if((tmi&3)==3) flag=1; tmi>>=1; } if(!flag) { uc[all++]=i; if((i|a[0])==a[0]) f[0][i]=1; } } for(int i=1; i<n; i++) for(int j=0; j<all; j++) if((a[i]|uc[j])==a[i]) { for(int k=0; k<all; k++) if(!(uc[j]&uc[k])) { f[i][uc[j]]=(f[i][uc[j]]+f[i-1][uc[k]])%MOD; } } int ans=0; for(int j=0; j<all; j++) ans=(ans+f[n-1][uc[j]])%MOD; printf("%d\n", ans); return 0; }
POJ 3255 A*算法求第k短路
第一次接触A*，了解了下。
—————————————————维基百科————————————————————————————
如果以 g(n)表示从起点到任意顶点n的实际距离，h(n)表示任意顶点n到目标顶点的估算距离，那么 A*算法的公式为：f(n)=g(n)+h(n)。 这个公式遵循以下特性：
—————————————————————————————————————————————————
#include <cstdio> #include <cstring> #include <iostream> #include <algorithm> #include <cstdlib> #include <cmath> #include <utility> #include <vector> #include <queue> #include <map> #include <set> #define max(x,y) ((x)>(y)?(x):(y)) #define min(x,y) ((x)>(y)?(y):(x)) #define INF 0x3f3f3f3f using namespace std; struct Edge { int y,w,ne; }e[200005]; int x,y,w,n,m; int be[5005],all; int h[5005]; bool vis[5005]; struct Point { int x,g; bool operator < (const Point T) const { return g+h[x]>T.g+h[T.x]; } }; void add(int x, int y, int w) { e[all].y=y; e[all].w=w; e[all].ne=be[x]; be[x]=all++; } void SPFA(int s) { queue< int > q; for(int i=0; i<=n; i++) { h[i]=INF; vis[i]=0; } h[s]=0; vis[s]=1; while(!q.empty()) q.pop(); q.push(s); while(!q.empty()) { int u=q.front(); q.pop(); vis[u]=0; for(int i=be[u]; i!=-1; i=e[i].ne) { int v=e[i].y; if(h[v]>h[u]+e[i].w) { h[v]=h[u]+e[i].w; if(!vis[v]) { vis[v]=1; q.push(v); } } } } } int Astar(int s, int t, int k) { SPFA(t); priority_queue< Point > q; while(!q.empty()) q.pop(); memset(vis,0,sizeof(vis)); Point cur,nxt; cur.x=s; cur.g=0; q.push(cur); while(!q.empty()) { cur=q.top(); q.pop(); if(cur.x==t && !--k) return cur.g; for(int i=be[cur.x]; i!=-1; i=e[i].ne) { nxt.x=e[i].y; nxt.g=cur.g+e[i].w; q.push(nxt); } } return -1; } void init() { all=0; memset(be,-1,sizeof(be)); } int main() { scanf("%d%d",&n,&m); init(); for(int i=0; i<m; i++) { scanf("%d%d%d",&x,&y,&w); add(x,y,w); add(y,x,w); } printf("%d\n",Astar(1,n,2)); return 0; }
标签：style blog http color os io ar for art
原文地址：http://www.cnblogs.com/Mathics/p/3955210.html
