[POI2014]Beads

题目大意：
　　有$n(n\leq10^6)$种颜色，第$i$种颜色有$c_i(\sum c_i\leq10^6)$个，指定第一个颜色为$a$，最后一个颜色为$b$，问对于一个长度为$m=\sum c_i$的序列，是否能构造出一个染色方案满足相邻的颜色不相同。如果能，试构造出一种方案。

思路：
　　贪心。如果序列中有多个元素，开头结尾颜色相同而这种颜色只有一个，则显然不存在合法方案。每次选取当前数量最多的不同于前一个颜色的颜色，如果有一样多的，就尽量取和最后一个位置的颜色相同的颜色。如果没有可以取的元素或者倒数第二个元素颜色只能选和最后一个元素相同的颜色就说明不存在合法方案。洛谷和BZOJ随便AC，然而SZKOpuł上怎么卡常都是T。

正常代码：

 #include<queue>
 #include<cstdio>
 #include<cctype>
 #include<algorithm>
 inline int getint() {
     register char ch;
     while(!isdigit(ch=getchar()));
     register int x=ch^'';
     while(isdigit(ch=getchar())) x=(((x<<)+x)<<)+(ch^'');
     return x;
 }
 const int N=1e6+;
 int a[N],c[N];
 std::priority_queue<std::pair<int,int> > q;
 int main() {
     const int n=getint();
     a[]=getint(),a[]=getint();
     int m=;
     for(register int i=;i<=n;i++) {
         m+=c[i]=getint();
     }
     std::swap(a[],a[m]);
     c[a[]]--;
     if(m!=) c[a[m]]--;
     if(c[a[]]<) {
         putchar('');
         return ;
     }
     for(register int i=;i<=n;i++) {
         if(c[i]) q.push(std::make_pair(c[i],i==a[m]?i+n:i));
     }
     for(register int i=;i<m;i++) {
         int x=q.top().second%n;q.pop();
         if(!x) x+=n;
         if(x==a[i-]) {
             if(q.empty()) {
                 putchar('');
                 return ;
             }
             int y=q.top().second%n;q.pop();
             if(!y) y+=n;
             q.push(std::make_pair(c[x],x==a[m]?x+n:x));
             a[i]=y;
             if(--c[y]) q.push(std::make_pair(c[y],y==a[m]?y+n:y));
         } else {
             a[i]=x;
             if(--c[x]) q.push(std::make_pair(c[x],x==a[m]?x+n:x));
         }
     }
     if(a[m-]==a[m]) {
         putchar('');
         return ;
     }
     for(register int i=;i<=m;i++) {
         printf("%d%c",a[i]," \n"[i==m]);
     }
     return ;
 }

卡常代码：

 #pragma GCC optimize(3)
 #pragma GCC optimize("unroll-loops")
 #pragma GCC target("sse,sse2,sse3,ssse3,sse4,popcnt,abm,mmx,avx,tune=native")
 #include<cstdio>
 #include<climits>
 #include<algorithm>
 #include<sys/mman.h>
 #include<sys/stat.h>
 typedef long long int64;
 class BufferedInputStream {
     private:
         char *buf,*p;
         int size;
     public:
         BufferedInputStream() {
             register int fd=fileno(stdin);
             struct stat sb;
             fstat(fd,&sb);
             size=sb.st_size;
             p=buf=reinterpret_cast<char*>(mmap(,size,PROT_READ,MAP_PRIVATE,fileno(stdin),));
         }
         char getchar() {
             return (p==buf+size||*p==EOF)?EOF:*p++;
         }
 };
 BufferedInputStream in;
 inline int getint() {
     register char ch;
     while(!__builtin_isdigit(ch=in.getchar()));
     register int x=ch^'';
     while(__builtin_isdigit(ch=in.getchar())) x=(((x<<)+x)<<)+(ch^'');
     return x;
 }
 const int N=1e6+;
 int a[N],c[N];
 class priority_queue {
     #define _left <<1
     #define _right <<1|1
     #define _par >>1
     private:
         int size;
         int64 val[N];
     public:
         priority_queue() {
             val[]=LLONG_MAX;
         }
         void push(const int64 &x) {
             register int p=++size;
             for(;val[p _par]<x;p=p _par) {
                 val[p]=val[p _par];
             }
             val[p]=x;
         }
         int64 top() const {
             return val[];
         }
         void pop() {
             register int p=;
             const int64 x=val[size--];
             for(;p _left<=size&&(x<val[p _left]||x<val[p _right]);p=val[p _left]>val[p _right]?p _left:p _right) {
                 val[p]=std::max(val[p _left],val[p _right]);
             }
             val[p]=x;
         }
         bool empty() const {
             return !size;
         }
     #undef _left
     #undef _right
     #undef _par
 };
 inline int64 make_pair(const int &x,const int &y) {
     return (int64)x<<|y;
 }
 inline int second(const int &x) {
     return x&0x7fffffff;
 }
 priority_queue q;
 int main() {
     const int n=getint();
     a[]=getint(),a[]=getint();
     int m=;
     for(register int i=;i<=n;i++) {
         m+=c[i]=getint();
     }
     std::swap(a[],a[m]);
     c[a[]]--;
     if(__builtin_expect(m!=,)) c[a[m]]--;
     if(__builtin_expect(c[a[]]<,)) {
         __builtin_puts("");
         __builtin_exit();
     }
     for(register int i=;i<=n;i++) {
         if(__builtin_expect(c[i],)) q.push(make_pair(c[i],i==a[m]?i+n:i));
     }
     for(register int i=;i<m;i++) {
         int x=second(q.top());q.pop();
         if(__builtin_expect(x>n,)) x-=n;
         if(__builtin_expect(x==a[i-],)) {
             if(__builtin_expect(q.empty(),)) {
                 __builtin_puts("");
                 __builtin_exit();
             }
             int y=second(q.top());q.pop();
             if(__builtin_expect(y>n,)) y-=n;
             q.push(make_pair(c[x],x==a[m]?x+n:x));
             if(__builtin_expect(--c[a[i]=y],)) q.push(make_pair(c[y],y==a[m]?y+n:y));
         } else {
             if(__builtin_expect(--c[a[i]=x],)) q.push(make_pair(c[x],x==a[m]?x+n:x));
         }
     }
     if(__builtin_expect(a[m-]==a[m],)) {
         __builtin_puts("");
         __builtin_exit();
     }
     for(register int i=;i<=m;i++) {
         __builtin_printf("%d%c",a[i]," \n"[i==m]);
     }
     __builtin_exit();
 }