bzoj1059题解

【解题思路】

　　因为只要验证可行性，所以考虑行和考虑列是等价的，故我们只考虑行的交换操作。

　　这样，拆一波点，把每一行拆成两个点，左边为原交换行，右边为目标交换行，原问题等价于能否对这个二分图进行完全匹配，题中给定邻接矩阵即为二分图的邻接矩阵。

　　于是直接匈牙利一波，复杂度O(n³)。

【参考代码】

 #include <cctype>
 #include <cstdio>
 #define REP(I,low,high) for(register int I=(low);I<=(high);I++)
 #define PER(I,high,low) for(register int I=(high);I>=(low);I--)
 inline int getint()
 {
     char ch=getchar();
     for(;!isdigit(ch)&&ch!='+'&&ch!='-';ch=getchar());
     bool impositive=ch=='-';
     if(impositive)
         ch=getchar();
     int result=;
     for(;isdigit(ch);ch=getchar())
         result=(result<<)+(result<<)+ch-'';
     return impositive?-result:result;
 }
 template<typename integer> inline int write(integer n)
 {
     integer now=n;
     bool impositive=now<;
     if(impositive)
     {
         putchar('-');
         now=-now;
     }
     char sav[];
     sav[]=now%+'';
     int result=;
     for(;now/=;sav[result++]=now%+'');
     PER(i,result-,)
         putchar(sav[i]);
     return result+impositive;
 }
 //Header Template
 #include <cstring>
 bool used[],map[][];
 int n,result[];
 bool find(int S)
 {
     REP(T,,n)
         if(map[S][T]&&!used[T])
         {
             used[T]=true;
             if(!result[T]||find(result[T]))
             {
                 result[T]=S;
                 return true;
             }
         }
     return false;
 }
 int main()
 {
     for(int T=getint();T--;)
     {
         n=getint();
         REP(i,,n)
             REP(j,,n)
                 map[i][j]=getint();
         memset(result,,sizeof(result));
         int ans=;
         REP(i,,n)
         {
             memset(used,,sizeof(used));
             ans+=find(i);
         }
         puts(ans==n?"Yes":"No");
     }
     return ;
 }